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Abstract 

We demonstrate that the dipole-hadron cross-section computed from the non-linear evo- 
lution equation for the Colour Glass Condensate saturates the Froissart bound in the case 
of a fixed coupling and for a small dipole (Q 2 ^> ^qcd)- That is, the cross-section increases 
as the logarithm squared of the energy, with a proportionality coefficient involving the pion 
mass and the BFKL intercept (a s -/V c /7r)4 In 2. The pion mass enters via the non-perturbative 
initial conditions at low energy. The BFKL equation emerges as a limit of the non-linear 
evolution equation valid in the tail of the hadron wavefunction. We provide a physical pic- 
ture for the transverse expansion of the hadron with increasing energy, and emphasize the 
importance of the colour correlations among the saturated gluons in suppressing non-unitary 
contributions due to long-range Coulomb tails. We present the first calculation of the satu- 
ration scale including the impact parameter dependence. We show that the cross-section at 
high energy exhibits geometric scaling with a different scaling variable as compared to the 
intermediate energy regime. 
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1 Introduction 



One of the striking features of the physics of strong interaction is that at high energies, cross 
sections are slowly, but monotonously, increasing with s (= the total center-of-mass energy 
squared). For instance, the data for the pp and pp total cross sections at high energy can be 
reasonably well fitted by both s s (with 6 ~ 0.09) and In 2 s, although the In 2 s growth appears 
to be favored by recent investigations Q. 

At a theoretical level, it was proven many years ago that in the limit s — > oo, the hadronic 
total cross sections must rise no faster thanQ (7r/m 2 ) In 2 s |||, ||, |J]. This is the Froissart bound, 
which is a consequence of very general principles, such as unitarity, crossing and analiticity, 
but does not rely on any detailed dynamical information. So, in reality, this bound may very 
well be not saturated, although the measured cross-sections seem to do so (at least, in so far 
as they are consistent with a squared log dependence upon the energy). In fact, a simple, 
qualitative mechanism realising this In 2 s growth has been proposed by Heisenberg j|] already 
before Froissart bound has been rigorously proven (a similar argument is briefly discussed by 
Froissart Q). However, we are not aware of any field-theoretical implementation of this, or 
other, argument. (See also Ref. || for a recent review and more references, and Ref. [jjj for a 
recent analysis of the data.) 

The goal of this paper is to see how the Froissart bound can be consistent with modern 
pictures of high energy strong interactions. To keep the discussion as simple as possible while 
still encompassing the interesting physics, we shall consider the scattering of a colour dipole off a 
hadronic target at very high energy. The "colour dipole" may be thought of as a quark-antiquark 
pair in a colourless state, like a quarkonium, or a fluctuation of the virtual photon in deep 
inelastic scattering. More generally, arbitrary hadronic probes can be considered as collections 
of "colour dipoles" at least in the approximation in which the number of colours N c is large (so 
that a gluon excitation can be effectively replaced by a qq pair). For our approximations to be 
justified, we shall assume that the dipole is "small", in the sense that it has a small transverse 
size, rj_ <C 1/Aqcd, or a large transverse resolution: Q 2 = 1/r 2 3> A-qcd- 

On the particular example of the dipole-hadron scattering, we shall develop an argument 
involving both perturbative and non-perturbative features of QCD, like non-linear quantum 
evolution, parton saturation, and confinement, which will lead us to the conclusion that the 
total cross-section respects, and even saturates, the Froissart bound (at least, in the case of a 
fixed coupling, that we shall exclusively consider in this paper). In its essence, our argument 
may be viewed as a modern version of Heisenberg's original mechanism. But our main point is 
to place this mechanism within the context of our present theoretical understanding, and clarify 
to which extent perturbation theory may play a role in determining this result. In brief, we shall 
demonstrate that this mechanism is naturally realized by the solutions to non-linear evolution 



equations derived within perturbation theory, in Refs. [7], ||, |J, 1C], but with non-perturbative 
initial conditions at low energy. In particular, our analysis will confirm, clarify and extend 
the conclusions reached in an early study ]ll]] based on the GLR equation [|l2|, 13 1, which is a 
simplified version of the non-linear evolution equations that we shall use in this paper. 



1 The coefficient C in front of in 2 s (atot(s) < C in 2 s) has been ieft unspecified in the original analysis by Frois- 
sart The upper bound C = 7r/m 2 has been first derived by Lukaszuk and Martin fl], from general assumptions 
(unitarity, crossing and analiticity) on the pion-pion scattering amplitude. This value is consistent with an early 
argument by Heisenberg in which the Froissart bound is actually saturated : a to t(s) = (7r/m 2 ) In 2 s. 
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The fact that perturbation theory is the appropriate tool to describe the high energy be- 
haviour of hadronic cross-sections is by itself non-trivial, and deserves some comment. The first, 
and most certainly true, objection to the use of weak coupling methods comes from the fact 
that the Froissart bound involves the pion mass, and this surely must arise from confinement. 
And, indeed, this is how the pion mass enters also our calculations: via the non-perturbative 
initial condition that we shall assume, and which specifies the impact parameter dependence of 
the dipole-hadron scattering amplitude at low energy. But starting with this initial condition, 
we shall then study its evolution with increasing energy within perturbation theory, and show 
that the resulting cross-section saturates the Froissart bound at high energy. 

The adequacy of perturbation theory even for such a limited purpose — the study of the 
quantum evolution of the cross-section with s — is still non-trivial ||, and has been in fact 
disputed in recent pape rs |i~4| ] . It is first of all clear that ordinary perturbation theory, in the form 
of the BFKL equation |l5| — which resums the dominant radiative corrections at high energy, 
but neglects the non-linear effects associated with high parton densities — , fails to describe the 
asymptotic behaviour at high energy: The BFKL equation predicts a power-law growth for the 
total cross section: a ~ s WQs , with u> = 41n2 and a s = a s N c /ir, which clearly violates the 
Froissart bound. Besides, with increasing energy, the solution to the BFKL equation "diffuses" 
towards smaller and smaller transverse momenta, thus making the applicability of perturbation 
theory questionable. 

But the objections in Ref. [14] actually apply to the more recently derived non-linear evo- 
lution equations j?], |8| pi, in which non-linear effects cure the obvious pathologies of the 
BFKL equation. (Alternative derivations of some of these equations, at least in specific limits, 
have been given in Refs. |12| , |i~3| , 16, [D]]. See also Refs. [18, 19] for recent reviews and more 
references.) Specifically, the non-linear effects associated with the high density of gluons in the 
hadron light-cone wavefunction (the "Colour Glass Condensate" [|l(J) lead to gluon saturation 
]l2| , 2C, 21, |2^[, with important consequences for the high-energy scattering: First, this intro- 
duces a hard intrinsic momentum scale, the "saturation scale" , which is a measure of the density 
of the saturated gluons in the impact parameter space, and grows like a power of the energy. 
This limits the infrared diffusion [^] and thus provides a better justification for using weak 
coupling methods at high energy. Second, this ensures the unitarization of the dipole-hadron 



scattering amplitude at fixed impact parameter [|24j, |25|, [2£|, 27, |8j. That is, with increasing 
energy, the hadron eventually turns "black" (i.e., the dipole is completely absorbed, or the scat- 
tering amplitude reaches the unitarity limit) at any given point b± in the impact parameter 
space. 

However, by itself, the unitarization at given b± is not enough to guarantee the Froissart 
bound for the total cross-section, which involves an integration over all the impact parameters. 
Indeed, as a quantum mechanical object, the hadron has not a sharp edge, but rather a dif- 
fuse tail, so, with increasing energy, the "black disk" can extend to larger and larger impact 
parameters. The radial expansion of the black disk is controlled by the scattering in the sur- 
rounding "grey area", where the gluon density is relatively low and the BFKL equation still 
applies. Given the problems of the BFKL equation alluded to before, it is not a priori clear 
whether this expansion is slow enough to ensure that the Froissart bound is respected. This 
would require a black disk radius which increases at most logarithmically with s. In Heisenberg's 
original argument (but using the current terminology), such a logarithmic increase was ensured 
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by a compensation between the power-law increase of the scattering amplitude with s and its 
exponential decrease with b = \b±\. Such an exponential fall-off at large impact parameters is, of 
course, a true property of full QCD, and also the crucial assumption about our initial condition, 
but it is not clear whether this property can be preserved by the perturbative quantum evolution, 
which involves massless gluons and therefore long range interactions. 

In fact, it was the main point of Ref. gj to argue that the Coulomb tails associated with the 
saturated gluons should replace the exponential fall-off with b = \b±\ of the initial distribution 
by just a power-law fall-off, which would be then too slow to ensure the Froissart bound: the 
corresponding black disk radius would increase as a power of s. 



However, as we shall explain in this paper, the argument in Ref. [14| is irrelevant for the 
problem at hand, and also incorrect in its original formulation] 2 ]. Part of the confusion in Ref. 
EM comes from non-recognizing that the saturated gluons are actually colour neutral over a 



relatively short distance (of the order of the inverse saturation scale) [22, 18], and thus cannot 
produce Coulomb tails at large distances. Rather, the colour field produced by the saturated 
gluons is merely a dipolar field, whose fall-off with b is sufficiently fast to respect the Froissart 
bound for the scattering of an external dipole. (In Ref. [14], this was masked by the fact that the 
authors were truly computing the scattering of a coloured external probe, although their verbal 
arguments were formally developed for a "dipole".) As we shall see, for the relevant impact 
parameters, this long-range dipole-dipole scattering is subleading at high energy as compared 
to the short-range scattering off the local sources. 

More precisely, we shall demonstrate that the region which controls the evolution of the 
cross section is the "grey area" outside, but close to the black disk. In this area, we expect 
perturbation theory to apply, since the local saturation scale is much larger than Aqqd- For an 
incoming dipole in this area, the dominant interactions are those with the non-saturated colour 
sources within a "saturation disk" (i.e., a disk with radius equal to the inverse saturation scale) 
around the impact parameter b±. This is a consequence of two physical facts: a ) the non-linear 
effects limit the contribution of the distant colour sources, and b) being colourless, the dipole 
couples only to the local electric field (as opposed to the long-range gauge potentials), so it is 
less influenced by colour sources which are far away. 

Since, moreover, the local saturation length is much shorter than the typical scale ~ 1/Aqcd 
for transverse inhomogeneity in the hadron (this is where the condition that the dipole is "small" , 
i.e., Q 2 3> Aq CD , is essential), it follows that the quantum evolution proceeds quasi-locally in the 
impact parameter space. This in turn implies that, within the grey area, the b ± -dependence 
of the scattering amplitude factorizes out, and is therefore determined by the initial condition 
at low energy. On general physical grounds, we shall assume this initial condition to have an 
exponential fall-off as e~ 2m ^ b at large distances (indeed, pion pairs must control the long distance 
tail of the hadron wavefunction; see, e.g., [^] and Refs. therein). Then an argument similar 
to the original one by Heisenberg can be used to conclude that the total cross-section increases 
like In 2 s. 

The coefficient in front of In 2 s in our final result is also interesting, as it reflects the subtle 
interplay between the perturbative and non-perturbative physics contributing to this result. 

2 See however the new preprint 0] where a modified version of this argument has been presented; we shall 
comment on this new argument in the Note added at the end of this paper. 
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Specifically, we shall find that, for any hadronic target, 
vr (uja s ) 2 2 

a w ^ — In s as s — > oo, f 1.11 

2 m~ 

where the pion mass in the denominator enters via the exponential fall-off of initial condition, 
while the factor ua s = 4(ln2)a s N c /ir in the numerator is recognized as the "BFKL intercept". 
This latter comes up because, in deriving this result, we will have to consider the solution to 
the BFKL equation at large energy for fixed Q 2 . This is of course the limit for which the 
BFKL equation has been originally proposed |fL5|| , but not also the limit used in more recent 
applications of this equation within the context of saturation (e.g., in studies of the saturation 



scale [21, |3(J, or of the "geometric scaling" ]29fl) for a homogeneous hadron. 

The difference with Refs. pl| , |29| , pOfl occurs because we consider here a different physical 
problem: Rather than studying the quantum evolution at a fixed impact parameter — which 
would then limit the applicability of the BFKL equation to not so high energies, such that the 
local saturation scale remains below Q 2 — , we rather follow the expansion of the black disk 
with increasing s, and use the BFKL equation only in the outer grey area at sufficiently large 
b, where the gluon density remains small even when the energy is large. In other terms, by 
increasing the energy at fixed Q 2 and simultaneously moving towards larger impact parameters, 
one always finds a corona where the local saturation scale is much smaller than Q, but much 
larger than Aqcd- For points b± in this corona, we can prove the factorization of the scattering 
amplitude into a b± -dependent "profile function" which is determined by the initial condition 
at low energy, and an energy- and Q 2 -dependent factor which can be computed by solving the 
homogeneous (i.e., no b± -dependence) BFKL equation. In the high energy limit at fixed Q 2 , 
this calculation yields the cross-section in eq. fll.l|) . 



In addition to eq. (LI), we shall find a variety of new features within our approach. For 
example, we shall compute the impact parameter dependence of the saturation scale, and dis- 
cover an entirely nontrivial structure for the radial distribution of matter inside the hadron, and 
its evolution with increasing energy. Related to that, we shall find that the geometric scaling 
arguments which have been used to characterize deep inelastic scattering at high energies be- 
come modified by the appearence of two different scales (associated both with gluon saturation) 
which govern the geometric scaling of the total cross-section in different ranges of the energy. 
As a result of our analysis, we shall derive an intuitive picture for the expansion of the hadron 
in the transverse plane. 

The outline of the paper is as follows: 

In the second section, we qualitatively and semi-quantitatively discuss how the Froissart 
bound becomes saturated within the framework of our knowledge of gluon saturation, non- 
linear evolution, and its linearized, BFKL, approximation. This discussion introduces the main 
arguments to be demonstrated by the technical developments in the nextcoming sections. 

In the third section, we study the properties of a non-linear evolution equation for the 
scattering amplitude, the Balitsky-Kovchegov (BK) equation [0, We argue that the dominant 
contribution to the scattering in the grey area comes from virtual dipoles whose size is much 
smaller then the saturation length. This leads us to conclude that the corresponding scattering 
amplitude factorizes in the way alluded to before. In the rest of the third section, we explore 
the physical significance of this factorization by using the efective theory for the Color Glass 
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Condensate [[ji], 18|. We show that an essential ingredient for factorization and Froissart bound 
is the colour-neutrality of the saturated gluons within the black disk. 

In the fourth section, we compute the radius of the black disk using the solution to the 
BFKL equation. Then, we extend our results to compute the impact parameter dependence of 
the saturation momentum. We identify two range of values of the impact parameter at which 
also the saturation scale factorizes, i.e., it is the product of an exponentially decreasing function 
of b± times a factor increasing like a power of s. These two ranges correspond to the two limits 
of the BFKL solution alluded to before: For b± sufficiently close to the center of the hadron, 
the increase with the energy is the same as for the saturation scale of a homogeneous hadron, 
previously studied in Refs. pl| , [2^ , 3C]. This comes from a solution to the BFKL equation in 
the intermediate regime where In Q 2 ~ a s In s. On the other hand, near the edge of the hadron, 
the increase with s is rather controlled by the high energy solution at InQ 2 <C a s lns. 

In the fifth section, we discuss the implication of the different behaviours found at small and 
large impact parameters for geometrical scaling in deep inelastic collisions. 

In the last section, we summarize our results and present our conclusions. 



2 Saturation scale and the Froissart bound 

The dipole-hadron collision will be considered in a special frame, the "dipole frame" [24, [TsR , 
in which the physical interpretation of our results becomes most transparent: This is the frame 
in which the effects of the quantum evolution are put solely in the wavefunction of the hadron, 
which carries most of the total energy. This being said, it should be stressed that our final 
results are independent of this choice of the frame — although their interpretation may look 
different in other frames — since at a mathematical level they are based on boost invariant 
equations, namely, the non- linear evolution equation for the scattering amplitude 0, |sj, [llj and 
its linearized, BFKL |fH[ ], approximation. In fact, the relevant non-linear equation — namely, 
the Balitsky-Kovchegov (BK) equation to be presented in Sect. 3 — has been independently 
derived in the hadron rest frame [0, ||, where the evolution refers to the incoming dipole wave- 
function, and in the infinite momentum frame (which for the present purposes is equivalent to 
the "dipole frame" alluded to before) , from the evolution of the "Colour Glass Condensate" [[h]] . 

So, let us consider an incoming dipole of transverse size r±, with Q 2 = l/r\ 3> Aq CD , which 
scatters off the hadron target at large invariant energy squared s, or rapidity gap r = ln(s/Q 2 ). 
In the dipole frame, most of the total energy is carried by the hadron, which moves nearly at 
the speed of light in the positive z direction. Moreover, any further increase in the total energy 
is achieved by boosting the hadron alone. Thus, the dipole rapidity Tdipole is constant, and 
chosen such as a s Tdi po ie "C 1> so that we can neglect higher Fock space components in dipole 
wavefunction: the dipole is just a quark-antiquark pair, without additional gluons. On the other 
hand, the hadron wavefunction has a large density of small-x gluons, which increases rapidly 
with r. Here, x is the longitudinal momentum fraction of the gluons which participate in the 
collision, and is related to the rapidity r as r = ln(l/x). 

In this special frame, the unitarization effects in the dipole-hadron collision can be assimilated 
to the saturation effects in the target wavefunction. This is so since it is the same momentum 
scale, namely the saturation scale Q 2 , which sets the border for both types of effects. 

A priori, Q s is an intrinsic scale of the hadron, proportional to the gluon density in the 
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transverse plane at saturation [12, 13, pH, FL8| : 



Qlij,b x ) ~ ^-xG(x,Q2, 6±) . (2.1) 
In this equation, 

is the number of gluons with longitudinal momentum fraction x and transverse size Axj_ ~ 1/Q 
per unit rapidity and per unit transverse area, at transverse location b±. The more standard 
gluon distribution xG(x, Q 2 ) is obtained by integrating xG(x, Q 2 , b±) over all the points b± in 
the transverse plane. Note that the distribution in b± is indeed a meaningful quantity since 
we consider gluons with relatively large transverse momenta, Q 2 3> Aq CD , which are therefore 
localized over distances Ax± much smaller than the typical scale for transverse inhomogeneity 
in the hadron, namely, Ab± ~ I/Aqcd- 

The saturation scale separates between two physical regimes: At high transverse momenta 
k± 3> Qs( T , b±), we are in the standard, perturbative regime: The gluon density is low, but it in- 
creases very fast, (quasi) exponentially with r, according to linear evolution equations like BFKL 



[15] or DGLAP |pj] |. At low momenta k± 5i Q s ( T ,b±), the non-linear effects are strong even if 
the coupling is weak, and lead to saturation: The gluon phase-space density is parametrically 
large, dN/drd 2 k±d 2 b± ~ l/a s , but increases only linearly with r [[H|, |2^, |l8|| . 

From eq. (p.lfl, one expects the saturation scale to increase rapidly with r, so like the gluon 



distribution at high momenta. For a hadron which is homogeneous in the transverse plane (no 



dependence upon b±), the r-dependence of Q s is by now well understood |2l], 2£, [3C], [3^, [33|, p3{ ], 
and will be reviewed in Sect. 4 below. Understanding the r- and 6-dependences of the saturation 
scale in the general inhomogeneous case is intimately related to the problem of the high energy 
behaviour of the total cross-section, so this will be a main focus for us in this paper. 

To appreciate the relevance of the saturation scale for the dipole scattering, note that a small 
dipole, i.e., a dipole with transverse size r± <C 1/Q S (t, b± ) , where b± is the impact parameter, 
couples to the electric field created by the colour sources in the target. Thus, its scattering 
amplitude is proportional to the correlator of two electric fields, which is the same as the gluon 
distribution xG(x, Q 2 ,b±) with Q 2 ~ \jr\ @ : 



TT 2 OL C 

Mr (r± ,b ± ) ~ r± j^ZY xG(x, \jr\ , b± ) (single scattering). (2.3) 

This equation, together with the solution xG(x, 1/tj_, b±) to the BFKL equation, predicts an 
exponential increase of the scattering amplitude with r. If extrapolated at high energy, this 
behaviour would violate the unitarity requirement J\f T < 1. However, eq. ( |2.3D assumes single 
scattering, so it is valid only as long as the gluon density is low enough for the condition 
■A/" T (r_L, <C 1 to be satisfied. At high energies, where the gluon density is large, multiple 
scattering becomes important, and leads to unitarization. Assuming that the successive collisions 
are independent, one obtains [24] : 

M T (r±,b±) ~ 1 — exp | — r\ x G(x, l/rj_, (multiple scattering), (2.4) 
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which clearly respects unitarity. 

So far, this only demonstrates the role of multiple scattering in restoring unitarity. The deep 
connection to saturation follows after observing that the condition for multiple scattering to be 
important — that is, that the exponent in eq. (|2.4j) is of order one — is the same as the condition 



(2.1) for gluon saturation in the hadron wavefunction. This is natural since the dipole is a direct 
probe of the gluon distribution in the hadron, so the non-linear effects in the dipole-hadron 
scattering and in the gluon distribution become important at the same scale. But this also 
shows that, in the nondinear regime at Q 2 J5 Ql{ T i b±) one cannot assume independent multiple 
scatterings, as in eq. fl2.4|): Rather, the dipole scatters coherently off the saturated gluons, with 



a scattering amplitude which, unlike eq. (2.4), cannot be related to the gluon distribution (a 



2-point function) alone, but involves also higher n-point functions. This amplitude satisfies a 
nondinear evolution equation to be discussed in Sect. 3. 

But it is nevertheless true that, as suggested by eq. fl2,4|) , the scattering amplitude becomes 
of order one in the saturation regime at Q 2 ^ Qs( T > (This has been verified via both analytic 



[32, 22, 29] and numerical [32, 33, |23J investigations of the nondinear evolution equation.) Thus, 
when increasing the energy at fixed Q 2 = l/r\, the unitarity limit M T (r±, b±) = 1 is eventually 
reached at any given impact parameter b± : for the incoming dipole, the hadron looks locally 
"black". 

However, the unitarization of the local scattering amplitude is not enough to guarantee the 
Froissart bound for the total cross-section. The latter is obtained by integrating the scattering 
amplitude over all impact parameters: 

a(r,r ± ) = 2 [ d 2 b ± M T (r ± ,b ± ). (2.5) 



It is easy to see that difficulties with the Froissart bound can arise only because the hadron does 
not have a sharp edge. Indeed, if in transverse projection the hadron was a disk of finite radius 
i?0) then for sufficiently large energy it would become black at all the points within that disk, 
and the total cross-section would saturate at the geometrical value 27ri?g. 

But in reality a hadron is a quantum bound state of the strong interactions, so its wave- 
function has necessarily an exponential tail, with the scale set by the lowest mass gap in QCD, 
that is, the pion mass. Specifically, in the rest frame of the hadron, the distribution of matter 
is typically of the Woods-Saxon type: 

P(r) = TTi (2-6) 



1 + exp (r=^j 



where Rq is the typical radial size of the hadron under consideration (this increases as A 1 / 3 for 
a nucleus with atomic number A), while the thickness a = l/2m n is universal (i.e., the same 
for all hadrons). This latter involves twice the pion mass because of isospin conservation: At 
high energy, one probes the gluons in the hadron wavefunction, and gluons have zero isospin, so 
they couple to the external probe via the exchange of (at least) two pions. It is therefore 2m n 
which controls the exponential fall off of the scattering amplitude, or of p(r), at large distances: 
p(r) oc e -^( r -Ro) f or r _ Rq > i/2m n . 

Since, moreover, small-x gluons have large longitudinal wavelengths, a high energy scattering 
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is sensitive only to the distribution integrated over z, that is, to the transverse profile function : 



/ dzp(Jb\ + : 

S ^ = ~ rTT ( \ . ( 2 - 7 ) 

Jdzp{z) 

(b = \b±\) which is normalized at the center of the hadron: S(b = 0) = 1. Note that, independent 
of the detailed form of p(r) in the central domain at r < Rq, the function S(b) decreases 
exponential^, 5(6) ~ e -*n*Q>-Ro) i f or b-R > l/2m*. 

Based on these considerations, we shall assume that the scattering amplitude for the (rel- 
atively low energy) dipole-hadron scattering in the target rest frame — which is our initial 
condition for the quantum evolution with r — has the following factorized structure: 

K (r ±1 b ± ) =K (r±)S(b), (2.8) 

with To = Tdipoie and the profile function S(b) introduced above. At low energy and high Q 2 , the 
factorization of the b ± -dependence is natural, since consistent with the DGLAP equation (see, 



e.g., p^l ). For what follows, the crucial feature of the initial amplitude (2.8) is its exponential 
fall off ~ e -^(b-R ) at large distances b > R . 

Starting with this initial condition, we increase the energy by boosting the hadron to higher 
and higher rapidities. Clearly, the gluon distribution xG(x, 1/tj_, b±) at any b± will increase with 
r, and the BFKL approximation suggests that this increase should be exponential. Thus, even 
points b± which were originally far away in the tail of the hadron wavefunction (b — Rq 1/277^), 
and did not contribute to scattering at the initial rapidity To, will eventually give a significant 
contribution, and even become black when the energy is high enough. That is, with increasing 
r, the black disk may extend to arbitrarily large impact parameters. 

At this point, it is useful to introduce some more terminology. By the "black disk" we mean 
the locus of the points b± in the transverse plane at which the unitarity limit A^-(rj_,6j_) = 1 
has been reached in a dipole-hadron collision at rapidity r and transverse resolution Q 2 = 
Equivalently, the condition Q 2 < Q 2 (r,b±) is satisfied — i.e., the gluons with transverse 
momenta ~ Q are saturated — at all the points in the black disk. Given the shape of the initial 
matter distribution Q2.G| ) — which is isotropic and decreases from the center of the hadron 
towards its edge — it is clear that the "black disk" is truly a disk, with center at b = 
and a radius R(t,Q 2 ) which increases with r and decreases with Q 2 . The black disk radius 
is determined by any of the two following conditions (for more clarity, we shall often rewrite 
J\f T (Q 2 ,b±) =AT T (r± = 1/Q,b±) in what follows): 

Af T (Q 2 , b±) = k for b = R(t,Q 2 ), (2.9) 



or 



Q2M±) = Q 2 for b = R(t,Q 2 ), (2.10) 

which are equivalent since, in turn, the saturation scale is defined by: 

K(Q 2 , b±) = k for Q 2 = Q 2 (r,b ± ). (2.11) 
3 But power law corrections to this exponential decrease may be numerically important |28|. 
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In these equations, k is a number smaller than one, but not much smaller (e.g., k = 1/2), whose 
precise value is a matter of convention^. For qualitative arguments, and also for quantitative 
estimates at the level of the approximations to be developed below, one can take k = 1. 

We shall also need below the "edge of the hadron" at rapidity r, by which we mean the radial 
distance Rh{t) at which the saturation scale becomes of order Aqcd (this would correspond to 
the black disk seen by a large dipole with resolution Q 2 ~ Aq CD , e.g., a pion) : 

Q 2 s (r,b ± ) = A 2 QCD for b = R H (r). (2.12) 

With these definitions at hand, we now return to the discussion of the Froissart bound. At 
sufficiently high energy, the total cross-section is dominated by the contribution of the black 
disk (this will be verified in Sect. 5) : a ~ obd with 

°bd(t,Q 2 ) = 2 J d 2 bj_ 0(i?(r, Q 2 ) - b)M T {Q 2 , b ± ) » 2ttR 2 (t,Q 2 ). (2.13) 

Thus, the question about the Froissart bound becomes a question about the expansion of the 
black disk with r : To respect this bound, R(t, Q 2 ) must grow at most linearly with r. 

One can easily construct a "naive" argument giving such a linear increase (this is similar in 
spirit to the old argument by Heisenberg ||): Starting with an initial distribution like (]2.6| ), 
assume that, with increasing r, the gluon density increases in the same way at all the points b± 
(outside the black disk), so that the b± -dependence of the scattering amplitude factorizes out, 
and is fixed by the initial condition: 

K(Q 2 ,b ± ) « S(b ± )Af T (Q 2 ). (2.14) 

Let us furthermore assume that the function M T (Q 2 ) at large r is given by standard perturbation 
theory, that is, by the solution to the BFKL equation at high energy |15| : J\f T {Q 2 ) oc e ujasT , 
where u = 4 In 2 and a s = N c a s /n. Under such (admittedly crude) assumptions, the scattering 
amplitude at large r and b 3> Rq is given by: 

A/V(Q 2 , &]_) « e^ T e~ M , (2.15) 

V Q 

where we have also included the leading Q 2 -dependence of the asymptotic BFKL solution at 
high energy [ |15[ . (A 2 is some arbitrary reference scale, of order Aq CD .) This expression together 



with the saturation condition (2.9) imply: 



and the resulting cross-section saturates the Froissart bound indeed: 

~ / I . (,'- \ n / .rc\ x \ 2 2 












) 2 ~l 






Km n J 



a( > T,Q ^ ~ 2^2 \^ sT - o ln T2 I ~o( — ) T as r^oo. (2.17) 



4 The separation between the saturation regime at small Q 2 and the low-density regime at high Q 2 being not 
a sharp one, ther e is some ambiguity in defining the borderline Q 2 = Q 2 (r,bj_). This is fixed by choosing the 



number k in eq. (2.11) 
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Our main objective in this paper will be to show that the above, seemingly naive, argument 
is essentially correct, and the results in eqs. ( 2.16| ) and ( 2.171) are truly the predictions of the 



non-linear evolution equation for N T (Q 2 , b±) at sufficiently large r. This is non-trivial since the 
naive argument might go wrong for, at least, two reasons: 

i) At very high energies, the non-linear effects become important, and the use of the BFKL 
equation becomes questionable. For instance, the unitarization of the local scattering amplitude 
J\f T (Q 2 ,b±) is precisely the result of such non-linear effects, which are taken into account by 
replacing the BFKL equation with the BK equation. 

ii) Although non-linear, the quantum evolution described by the BK equation remains per- 
turbative, so it involves massless gluons and long-range effects which could not only invalidate the 
factorization property ( [2.14D , but also replace the exponential fall-off of the initial distribution 
by just a power-law fall-off (an eventuality in which the Froissart bound would be violated). 

Nevertheless, as we explain now (and will demonstrate in Sect. 3 below), none of these two 
objections apply to the problem of interest. Indeed: 

i) However large is r, there exists an outer corona at R(t, Q 2 ) < b < Rh(t) where the hadron 
looks still "grey", i.e., where J\f T (Q 2 , b±) <C 1 and the BFKL approximation applies. It is this 
"grey area" which controls the expansion of the black disk, and therefore the evolution of the 
total cross-section at high energy. 

ii) The quantum evolution within the grey area is quasi- local in b± , because of the non- linear 
effects which limit the range of the relevant interactions to Ab± <C 1/Q s (t, b±). 

To be more specific, note that, in order to study the expansion of the black disk with r, 
one needs to consider the evolution of the scattering amplitude J\f T (Q 2 ,b±) at points b± which 
lie outside the black disk, but relatively close to it. Indeed, when r — ► r + dr with a s dr ~ 1 
(which is the typical increment in the high energy regime of interest: a s r 3> 1), the black disk 
expands by incorporating the points b± within the range R < b < R + dR with R = R(r, Q 2 ) 



and dR ~ cf. eq. (2.16). Such points are sufficiently far away from the black disk for 

the local saturation scale to be small compared to Q 2 — that is, they are in the "grey area" 
— , but also sufficiently far away from the edge of the hadron (cf. eq. ( |2,12| )) for Q s (T,b±) 
to be a "hard" scale. That is, the following conditions are satisfied for any b± of interest: 
A-q CD -C Q 2 (r,b) -C Q 2 . Both inequalities are important for our argument, as we explain now: 
The fact that Q 2 s (r,b) <C Q 2 ensures that the dominant contribution to the evolution 
dJ\f T (Q 2 , b±) of the scattering amplitude comes from nearby colour sources, i.e., from the sources 
which are located within a saturation disk around b± (see also Fig. [lj) : 

\z±-b±\ « - 1 . , (2.18) 
Qs{T, b±) 

and therefore lie themselves inside the grey area. This is so because sources which lie further 
away are shielded by the non-linear effects. Besides, being a colour singlet, the dipole is not 
sensitive to the long-range gauge potentials. 

The fact that 1/Q s (t, b±) <C 1/A.qcd implies that the transverse inhomogeneity in the 
hadron can be neglected when computing the contribution of such nearby sources to cbV r (Q 2 , b±). 
That is, all the relevant sources act as being effectively at the same impact parameter, equal to 
b±. This explains the factorized expression ( 2.14 ) f° r t ne scattering amplitude. 



Since, moreover, J\f T (Q 1 , z±) <C 1 for any z± satisfying ( 2.18] ), it follows that the function 



N T (Q 2 ) in eq. ( 2.14| ) can be computed by solving the linearized (and homogeneous) version of 
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Figure 1: A pictorial representation of the dipole-hadron scattering in transverse projection 
(only half of the hadron disk is shown). The 6-dependence of the saturation scale illustrated by 
the lower plot is the one to be found in Sect. 4.3. 
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the BK equation, namely the BFKL equation without b± dependence. 

Clearly, it was essential for the previous arguments that the dipole is "perturbative" : Q 2 S> 
A.q CD . This ensures that the separation Rh(t) — R(t, Q 2 ) between the black disk and the hadron 
edge (i.e., the width of the "grey area") is sufficiently large for the condition Q 2 (r,b) 3> A.q CD 
to apply at all the points b± of interest. Besides, one can argue that Q 2 is the scale at which 
the QCD coupling should be evaluated (see the discussion in the Conclusions). 

A factorization assumption similar to eq. ( |2.14j ) has been already used in the literature, in 
particular, in relation with the Froissart bound |35|, and also as an Ansatz in the search for 
approximate analytical [32| or numerical solutions to the BK equation. But in previous work, 
this assumption was always based on experience with the (homogeneous) DGLAP equation, and 
not duly justified in the small-x regime. 

That such a factorization is highly non-trivial in the presence of long-range gauge interac- 
tions is also emphasized by a recent controversy about this point, put forward in Ref. [[14[|. 
Specifically, in Ref. [14| it has been shown that, as far as the scattering of a coloured probe 
off the hadron is concerned, the long-range fields created by the saturated gluons provide a 
non-unitarizing contribution to the respective cross-section. On the basis of this example, the 
authors of Ref. fT4[| have concluded that the non-linear BK equation provides "saturation with- 
out unitarization" . In Sect. 3.2 below, we shall carefully and critically examine the arguments in 
Ref. []l4|], and demonstrate that, for the physically interesting case where the external probe is a 
(colourless) dipole, there is no problem with unitarity at all. The long-range interactions between 
the incoming dipole and the saturated gluons give only a small contribution to the scattering 
amplitude in the grey area, because the saturated gluons form themselves a dipole (i.e., they are 
globally colour neutral), and the dipole-dipole interaction falls off sufficiently fast with the sep- 
aration between the two dipoles. The dominant contribution comes rather from the short-range 
scattering within the grey area (cf. eq. ( 2.18 )), for which the factorization assumption ( |2.14 ) is 
indeed justified. 



3 Quantum evolution and black disk radius 

In the effective theory for the Colour Glass [fill], the dipole-hadron scattering is described as 
scattering of the qq pair off a stochastic classical colour field which represents the small-x gluons 
in the hadron wavefunction. At high energy, one can use the eikonal approximation to obtain: 

M T (r ± ,b ± ) = l-S T (x ± ,y ± ), S T (x ± ,y ± ) = ^(tv(VHx ± )V(y ± ))) T , (3.1) 

with rj_ = x± — y± the size of the dipole and b± = (x± + y±)/2 the impact parameter (the quark 
is at xj_, and the antiquark at y±). The S"-matrix element S T ( x± , y± ) involves the Wilson lines 
(path ordered exponentials along the straightline trajectories of the quark and the antiquark) 
V* and V built with the colour field of the target hadron. For instance, 

V\x x ) = Pexplig f dna a Jx ± )t a ) , (3.2) 



o 
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where a rj (x±) is the stochastic "Coulomb field" created by color sources (mostly gluons) at 
rapidities r' < r, and has longitudinal support at (space-time) rapidity^ rj < r. Thus, the 
integral over rj in eq. ( |3.2D is in fact an integral over the longitudinal extent of the hadron (in 
units of space-time rapidity) seen by the external probe in a scattering at rapidity r. That is, 
the actual width of the hadron depends upon the energy of the collision. This is so since, with 
increasing energy, gluon modes with larger and larger longitudinal wavelengths participate in 
the collision, so that the hadron looks effectively thicker and thicker fllQf . 

The brackets in the definition (3.1) of the S'-matrix element refer to the average over all the 



configurations of the classical field with some appropriate probability distribution W T [a] : 

MV\x ± )V(y ± ))) T = J [da] tv(V^(x ± )V(y ± ))W T [a\. (3.3) 

This probability distribution is not known directly, but its variation corresponding to integrating 
out gluons in the rapidity window (r, r + dr) can be computed [34, 10]. This leads to a functional 



evolution equation for W r [a] whose precise form is not needed here (see Ref. [10] for details). 



Suffices it to say that, via equations like (3.3), the functional equation for W r [a] can be translated 
into an hierarchy of ordinary evolution equations for the n-point functions of the Wilson lines [pX(] . 
This procedure yields the same equations as obtained by Balitsky within a different approach, 
which focuses directly on the evolution of Wilson line operators 0. (The fact that the infinite 
hierarchy of coupled equations by Balitsky can be reformulated as a single functional equation 
has been first recognized by Weigert ||.) In the limit where the number of colours N c is large, 
a closed equation can be written for the 2-point function ( |3.l[) (with a s = N c a s /ir) : 

d _ , , _ f d 2 z ± (x ± - y ± ) 2 

-S T {x±,y±) = -a s 



dr J 2-k (x± - z±) 2 (y± - z±) 2 

x(s T (x±,y±) - S T (x±, z±)S T (z±,y±)y (3.4) 

The same equation has been derived independently by Kovchegov H within the Mueller's dipolc 



model |25] . We shall refer to eq. ( pA ) as the Balitsky-Kovchegov (BK) equation. 



3.1 Scattering in the grey area 

In this subsection, we shall study the scattering amplitude in the grey area, and prove the 
factorization property ( |2.14| ). For more clarity, we shall formulate our arguments at the level of 
the BK equation. But one should keep in mind that our final conclusions are not specific to the 
large N c limit: the same results would have been obtained starting with the general non- linear 
evolution equations in Refs. @, [9|, ^Q] . 

Specifically, we shall use eq. (|3.4| ) to demonstrate that the dominant contribution to dS T /dr 
in the grey area comes from short-range scattering, i.e. from points z± such that 

\z± -b ± \< 1 « -r— . (3.5) 
Q s (T,b±) Aqcd 



5 The space-time rapidity is defined as rj = ln(a; _ P + ), where x~~ — (t — z)/\/2 is the light-cone longitudinal 
coordinate, and P + is the light-cone momentum of the hadron. Light-cone vector notations are defined in the 
standard way, that is, v = (l/v / 2)(v° ± v 3 ). With the present conventions, the hadron is a right mover, while 
the dipole is a left mover. 
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As explained in Sect. 2, the impact parameters of interest are such that the following in- 
equalities are satisfied (with Q 2 = l/rj_): Q 2 3> Q 2 (r,b) S> Aq CD . That is, the dipole is small 
not only as compared to the typical scale for non-perturbative physics and transverse inhomo- 
geneity in the hadron, namely 1/Aqcd, but also as compared to the shorter scale 1/Q s (t, b±), 
which is the local saturation length. This implies that the dipole is only weakly interacting with 
the hadron: S T (x±,y±) ~ 1, or M T (r±,b±) <C 1. But this does not mean that we are a priori 



allowed to linearize eq. (3^) with respect to N T . Indeed, the r.h.s. of this equation involves an 
integral over all z±, so the virtual dipoles with transverse coordinates (x±, z±) or (z±,y±) can be 
arbitrarily large. In fact, we shall see below that the dominant contribution comes nevertheless 
from z± which is relatively close to in the sense of eq. Q3,5| ), but the upper limit 1/Q s (r, b±) 
in this equation is a consequence of the non-linear effects. 

To see this, it is convenient to divide the integral over z± in eq. ( |3.4[ ) into two domains 
("short-range" and "long-range"): 



(A) « 1/Q s (r,b ± ), (B) 1/Q s (r,b ± ) <. \z± 



(3-6) 



It is straightforward to compute the contribution of domain (B) to the r.h.s. of eq. (3^): In this 
range, — z±\ ~ \z± — y±\ ~ \z± — b±\, so the virtual dipoles are both relatively large, and 
therefore strongly absorbed. Thus, to estimate their contribution, one can set S r (x±,z±) « 
and S T (z±,y±) 



0, and approximate the ("dipole" |24|, 25], [T 



(x± - y±) 



(x ± -z ± ) 2 (y ± -z ± ) 2 (z ± -b ± )±- 
This gives (with u\ = (z ± - b ± ) 2 > 1/Q 2 s (r,b±)) : 

-5 T (x_L,y_L) 

(B) 



kernel in the BK equation as: 

(3.7) 



0t 



~ -a s r 2 ± S T (x ± ,y ± ) J 



du. 



IQ1 u ± 



Y(rlQ 2 s (r,b ± ))S T (x ± ,y ± ). 



(3.8) 



Since S T (x±,y±) is of order one for the small dipole of interest, we deduce the following order- 
of-magnitude estimate (which we write for dM T /dT, for further convenience) : 



d_ 

dr 



Mr(r±,h 



(B) 



Us rj_Q 2 s (T, b±). 



(3.9) 



Note that, even for this "long range" contribution, the integral in eq. ( |3.8[) is dominated by 
points u± which are relatively close to the lower limit 1/Q s (t, b±) ; this is so because the dipole 
kernel (|3.?1 ) is rapidly decreasing at large distances | z± — b± \ 3> r± . 

To evaluate the corresponding contribution of domain (A), we note first that, in this domain, 
all the dipoles are small, so the scattering amplitude is small, N T ^ 1, for any of them. It is 
therefore appropriate to linearize the r.h.s. of eq. (|3.4|) with respect to M T (below, u± = x± — z±): 



-a. 



(A) 



2n uj_(r± — u±Y z 



(3.10) 



x lM T (r ± ,b ± ) -J\f T (u±,b± 



u± - r ± 



I'L ~U±,b±- -y) 
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This is recognized as the BFKL equation in the coordinate representation. Since both u± and r± 
are small as compared to 1/Q s (t, b±), and therefore much smaller than 1/Aqcd, it is appropriate 
to neglect the hadron inhomogeneity when evaluating the r.h.s. of this equation. That is, all 
the functions N T in the r.h.s can be evaluated at the same impact parameter, namely b±. 

To obtain an order-of-magnitude estimate for the r.h.s. of eq. Q3.10), we need an estimate 



for the function M T (r±,b±) in the regime where r± <C 1/Q s (t, b±). An approximate solution 
valid in this regime will be constructed in Sect. ||. But for the present purposes, we do not need 
all the details of this solution. Rather, it is enough to use the following "scaling approximation" 

K(r±,b ± ) ~ (rlQ 2 s (r,b ± )) X , (3.11) 

with A < 1. In Ref. |29f| , this approximation has been justified for a homogeneous hadron (no 
dependence upon b±). In Sect. 5 below, we shall find that, in the regime of interest, geometric 
scaling remains true also in the presence of inhomogeneity. 

The highest value A = 1 corresponds to the "double logarithmic regime"^] in which the dipole 
is extremely small, ln(Q 2 /Q 2 (t, b±)) S> 1, or, equivalently, its impact parameter b± is far outside 
the black disk, b R(t,Q 2 ). Here, however, we are mostly interested in points b± which are 
not so far away from the black disk, since we would like to study how the latter expands by 
incorporating points from the grey area. In this regime, i.e., for b > R(t,Q 2 ) but such that 
(6 — R(t,Q 2 ))/R(t,Q 2 ) <C 1, the scattering amplitude is given by eq. ( |3.11 ) with a power A 
which is strictly smaller than one (see Sect. 5). 



To simplify the evaluation of eq. ( 3.10| ), we shall divide domain (A) in two subdomains, in 



which further approximations are possible: (A.I) one of the two virtual dipoles, say u±, is much 
smaller than the other oneQ : u± -C r±; (A. II) both virtual dipoles are larger than the original 
one, although still smaller than one saturation length: r± <C u± ~ \u± — r±\ <C 1/Q S (t, b± ) . 

In domain (A.I), the first and third "dipoles" in the r.h.s. of eq. ( |3.10| ) cancel each other 
(since r\ ~ (r± — u±) 2 ), and we are left with 



■^M T (r ± ,bj 



^/^^r^ («MM±)) A ~ a s (rlQ 2 (r,b ± )) X , (3.12) 



(A.I) 



where the factor of 2 takes into account that one can choose any of the two virtual dipoles as 
the small one. 

In domain (A. II), one can neglect M T (r±, b±) <C M T (u±, b±) ~ N T (u± —r±,b±), and obtain: 



■^M T (r ± ,b A 



- 2a s 

(A.II) 

= a. 



2 



l/Qs d 2 u±_ r_ 
2vr 



(u\Q 2 s (t, b ± ) 



(rlQ^)Y ^^^ , (3 , 3) 



Strictly speaking, there is no geometric sca ling in this regime, but for power counting purposes one can assume 



that M T (fj_, 6_l) oc r ± xG(x, l/r ± , b±) (cf. eq. (2.4)) is linear in r ± . Indeed, at very high Q , the gluon distribution 
xG(x, Q 2 ,b±_) is only weakly depen dent upon Q . 



7 That is, in the notations of eq. (3.4), the point z± is within the area occupied by the original dipole (x±,y±), 
and much closer to x± than to y±. 
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which is of the same order as the (A.I)-contribution ( |3.12| ) when A < 1, but is logarithmically 
enhanced over it, and also over the long-range contribution (3.9), when A — * 1 : 



d 

-Q-N T {r±,b±) 



~ a, 



(A.II) 



(rlQ 2 s (r,b ± ))ln 



_Ql(r,b ± ) 



when A = 1. 



(3.14) 



By comparing eqs. (|3.9[) and fl3.12|) -( [3.14 ), it should be clear by now that, for any A < 1, 
the short-range contribution, domain (A), dominates over the long-range one, domain (B). In 
other words, from the analysis of the non-linear BK equation, we have found that, for a "small" 
incoming dipole, the dominant contribution to the quantum evolution of M T (r± , b± ) comes 
from still "small" virtual dipoles {(x±, z±), (z±, y±)} (see Figure [J). This has two important 
consequences, (i) One can linearize the BK equation with respect to M T , as we did already in 



eq. (3.10). This gives the BFKL equation, (ii) One can ignore the transverse inhomogeneity in 
the BFKL equation. That is, one can replace eq. (|3.10| ) by 



d 

— N T {r±,b ± ) 



Q'.s 



d 2 z ± 



2vr (z± - x±) 2 (z ± - y ± ) 



:{j\f T (z± -x±,b±) +N T (z± -y±,b±) -N T (r ± ,b ± )}, 



(3.15) 



m 



in which all the amplitudes N T are evaluated at the same impact parameter, namely at b±. 

Note that, as compared to eq. ( 3. 13] ), there is no need to insert an upper cutoff ~ \/Q 
the integral in eq. ( 3.15j ) . This is so since, to the accuracy of interest, the solution Af T (r±, b±) 
to eq. ( 3. 15| ) is actually insensitive to such a cutoff. One can understand this on the basis of 
eq. fl3.13|) : For large r± (with rj_Ql(r, b±) <C 1 though), the solution has the "scaling" behaviour 
in eq. (l3.1l] ) with a power A which is strictly smaller than one (see the discussion in Sects. 4 and 
5). Then the integral in eq. (3.13) is dominated by points u± which are close to the lower limit 
rj_, i.e., by virtual dipoles which are not much larger than the incoming dipole. The dependence 
upon the upper cutoff ~ 1/Q S is therefore a sub leading effect, which can be safely ignored. 

We finally come to the last step in our argument: Since the dependence of eq. ( 3.15| ) upon 
b± is only "parametric", it is clear that the impact parameter dependence of the solution is 
entirely fixed by the initial condition. This, together with eq. ( |2.8| ), implies that M T (r±, b±) has 
the factorized structure: 



M T (r ± ,b ± ) = S(b ± W T (r ± ) 



(3.16) 



where S(b±) is the transverse profile of the initial condition, while M T (r±) satisfies the homo- 
geneous BFKL equation and will be discussed in Sect. ||. A brief inspection of the previous 
arguments reveals that the terms neglected in our approximations are suppressed by either 
powers of r±Q 2 (r,b±) (e.g., the long-range contribution (|3.9|) , or the cutoff-dependent term in 
eq. ( [3.13 )), or by powers of Aq CD /Q 2 (t, b±) (the inhomogeneous effects in eq. Q3.10 )). This 
specifies the accuracy of the factorized approximation in eq. (3.16). 



3.2 More on the saturated gluons 

In the previous subsection, we have seen that colour sources located far away from the impact 
parameter of the dipole, such as \z± — b±\ > 1/Q s (r, 6j_), do not significantly contribute to 
the scattering amplitude in the grey area. In what follows, we shall examine more carefully a 
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particular contribution of this type, namely, that associated with the saturated gluons within 



the black disk: \z±\ < R(t,Q 2 ). Indeed, it has been recently argued 14] that, by itself, this 
contribution would lead to unitarity violations. To clarify this point, we shall compute this 
contribution within the effective theory for the Colour Glass Condensate, where the physical 
interpretation of the result is transparent. The same result will be then reobtained from the 
BK equation. Our analysis will confirm that, for impact parameters b± within the grey area, 
this long-range contribution is indeed subleading, and can be safely neglected at high energies. 



As we shall see, non-unitary contributions of the type discussed in Ref. [14| appear only in 
the physically uninteresting case where the exernal probe carries a non-zero colour charge (as 
opposed to the colourless dipole). 

According to eqs. ( |3.lD and (|3.2| ), the scattering amplitude at rapidity r depends upon the 
Coulomb field a'} l (x±) at all the space-time rapidities r/ < r. In general, this is related to the 
colour sources in the hadron via the two-dimensional Poisson equation — V\a^{x±) = p^(x±), 
with the solution: 

<4(x±) = J ^ d 2 z ± {x 1 _\^ T \z i _) p«(z ± ). (3.17) 

In this equation, p^(x±) is the colour charge density (per unit transverse area per unit space- 
time rapidity) of the colour sources at space-time rapidity r]. Here, we are only interested in 
such colour sources which are saturated. To isolate their contribution, it is important to remark 
that these sources have been generated by the quantum evolution up to a "time" equal to r/, 
so the corresponding saturation scale is Q s (rj, b±) (and not Q s (t, b±)). Thus, the integration in 
eq. (|3 . 1 7| ) must be restricted to \z±\ < R(r], Q 2 ), with R(rj, Q 2 ) the black disk radius at rapidity 
7] (cf. eq. Q2.10D ), and Q 2 the typical momentum carried by the Fourier modes of a^(x±) (as 
usual, this is fixed by the transverse size of the incoming dipole). 

There is a similar restriction on the values of r/ : For given Q 2 , there is a minimum rapidity 
fi(Q 2 ) below which there is no black disk at all: R(rj,Q 2 ) = for r/ < fi(Q 2 )- This is the 
rapidity at which the black disk first emerges at the center of the hadron, namely, at which 
Q 2 (fi,b = 0) = Q 2 (see Sect. 4.4 below). Thus, in order to count saturated sources only, the 
integral over rj in Wilson lines like ( |3.2| ) must be restricted to the interval f\{Q 2 ) < i] < r. In 
Fig. the saturated sources (with momentum Q 2 ) occupy the lower right corner, below the 
dashed line which represents the profile of the black disk as a function of r/. 

The external point x± in eq. ( |3.17 ) is at the impact parameter of the quark (or the antiquark) 
in the dipole, so it satisfies \x±\ 3> R(t, Q 2 ) > R(r), Q 2 ) for any r] < r. It is therefore appropriate 
to evaluate the field ( |3.17| ) in a multipolar expansion: 

a a v (x ± ) = ( x± \^\0 ± }J R d 2 z ±P a v (z ± )-^(x ± \-^\0 ± )J R d 2 z ± z*p a v (z ± ) +••• 



(x ± \-^ T \0 ± }Q« + ^- T Vl+-.., ( o. I S ) 



where R = R{r] : Q 2 ), Q a is the total colour charge within the black disk, T> l a is the corresponding 
dipolar moment, etc. 

To compute the scattering amplitude ( |3.lD , one has to construct the Wilson lines and V 
with the field ( |3.18| ) and then average over a (or, equivalently, over p) as in eq. ( |3.3| ). In what 
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1\ T\ X 

Figure 2: The longitudinal profile of the hadron as it appears in a scattering at given r and Q 2 . 
The longitudinal coordinate is on the horizontal axis, and is measured in units of space-time 
rapidity. The tranverse coordinate is on the vertical axis. A longitudinal layer at rapidity rj is 
delimited for more clarity. The wavy line represents the colour field a v created at point x± by 
the (saturated) source p v at z±. The enclined dashed line represents the limit of the black disk, 
which increases linearly with 77, as we shall see in Sect. 4. 
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follows, it is more convenient to work with the probability distribution for p, i.e., W r [p]. In 
general, this distribution is determined by a complicated functional evolution equation, which 
is very non- linear g |n]]. 



However, as observed in Ref. [22], this equation simplifies drastically in the saturation 
regime, where the corresponding solution W r [/o] is essentially a Gaussian in p. This can be 
understood as follows: The non-linear effects in the quantum evolution enter via Wilson lines 
like eq. (|3.2| ). At saturation, the field a in the exponential carries typical momenta k± ~ Q s and 
has a large amplitude a ~ 1/g. Thus, the Wilson lines are strongly varying over a transverse 
distance 1/Q s (t, b±). When observed by a probe with transverse resolution Q 2 <C Q 2 {t, b±), 
these Wilson lines are rapidly oscillating and average to zero. Thus, at saturation, one can 
drop out the Wilson lines, and all the associated non-local and non-linear effects. Then, the 
probability distribution W T [p] becomes indeed a Gaussian, which by the same argument is local 
in colour and space-time rapidity, and also homogeneous in all the (longitudinal and transverse) 
coordinates. The only remaining correlations are those in the transverse plane, and, importantly, 



these are such as to ensure colour neutrality [18]. 



Specifically, the only non-trivial correlation function of the saturated sources is the two-point 
function, which reads [^] (see also Sect. 5.4 in Ref. flql ) : 

(p a v (z ± )p b v ,(u ± )) T = 8 ab 8( V - V ')X(z ± -u ± ), 

X(k ± ) = -k\. (3.19) 

For given rj and k"j_ ~ Q 2 , eq. ( |3.19| ) holds for points z± and u± within the black disk of radius 
R(r],Q 2 ). The crucial property of the 2-point function ( |3.19| ) is that it vanishes as k\ — ► 0. 
Physically, this means that, globally, the saturated gluons are colour neutral^], as anticipated: 
{Q a Q a ) = 0, where Q a is the total colour charge (at given rj) in the transverse plane. In fact, 
since the Wilson lines average to zero over distances Ab± <; 1/Q s (rj,b±), it follows that colour 
neutrality is achieved already over a transverse scale of the order of the saturation length: 

f d 2 x ± [ d 2 y ± (p$(x ± )f£(y ± )) = 0, (3.20) 

JAS_l JAS ± 

where AS± is, e.g., a disk of radius R > 1/Q s (r/, b±) centered at b±. 

This immediately implies that, as soon as the black disk is large enough, the overall charge of 
the saturated gluons vanishes, Q a = 0, so we can ignore the monopole field in eq. ( p,18| ). Here, 
"large enough" means, e.g., R{t,Q 2 ) > 1/Q, which guarantees that the radius of the black disk 
is larger than the saturation length 1/Q s ( r ],b±) at any b± within the disk and at any r] in the 
interval f\{Q 2 ) < t) < r (since then Q s (^,6j_) > Q). In these conditions, the dominant field of 
the saturated gluons at large distances is the dipolar field in eq. ( |3.18| ). 

Below, we shall need the two-point function of this field: 

(a$(x ± )afy(y ± )) T = 5 ab 5{r, - rj) lv (x ± - y ± ). (3.21) 
From eq. (|3.18j ), we deduce (with colour indices omitted, since trivial): 



8 Since the distribution of p is a Gaussian, the fact that (Q a Q a ) vanishes is equivalent to Q a = 0, which means 
colour neutrality indeed. 
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with (cf. eq. (pPD): 



d 2 u± 



Q° d 2 k, , ^ <9 2 A- 2 

(2vr) 2 SJb^fcj vr 



S»2R 2 (rj,Q 2 



(3.23) 



where i? = R(rj,Q 2 ), Q s = Q s (v^±) (with 6j_ = (z± + u±)/2 < R), and formal manipulations 
like integrations by parts or the use of the Fourier representation of the 5-function were permitted 
since RQ S 3> 1. Thus, finally, 



7tj(z_l,1/±) = ^ ^ 2 (^>Q 2 



(3.24) 



which, we recall, is valid only as long as x±, y± 3> R(r), Q 2 ). 

We are now in a position to compute the scattering amplitude ( |3.1| ) for the scattering between 
the incoming dipole and the dipolar colour charge distribution within the black disk. To this 
aim, we have to average the product tr(W (x±)V(y±)) over the Gaussian random variable a v 
with two-point function ( |3.21 ). The result of this calculation is well-known (see, e.g., [pjjfl ): 



S T (x±,y±) = exp 



9 2 C F 



dr] -f v (x±,x±) + 7r?(y±,y±) - ^ ri (x±,y±) 



(3.25) 



where C F = t a t a = (N 2 - 1)/2N C and the lower limit r in the integral is a shorthand for 
f\{Q 2 ) with Q 2 = \jr 2 _. The "tilde" symbol on S is to remind that this is not the total S- 
matrix element, but just the particular contribution to it coming from the saturated gluons. An 
immediate calculation using eqs. ( 3.24 ) and (|3.25|) yields (with a s = 2a s Ci?/7r ~ a s N c /ir in the 
large _/V c limit) : 



Sr(x±,y±) 



exp < -a s 



(x± - yi_f 
2x\y\ 



dr 1 R 2 (7 1 ,Q 2 ) 



~ exp <^ -a s I dr] R 2 (r], Q 2 ) \ , 



(3.26) 



where in the second line we have replaced in the denominator x\ ~ y\ ~ b\ (which is appropriate 
since b± 3> R(t,Q 2 ) 3> r±). The exponent in eq. ( 3.26| ) vanishes when the dipole shrinks to a 
point, r± — * 0. This is the expected dipole cancellation, manifest already on eq. ( |3.25 ). 



The previous derivation makes the physical interpretation of eq. (3.26) very clear: The 
exponent there is the square of the potential ~ gt a r±^R(r]) for the interaction between two 
dipoles — the "external dipole" of size r± and the dipole made of the saturated gluons (at a 
given space-time rapidity rj), with size R(rj) — separated by a large distance b. There exists one 



layer of saturated gluons at any r\ within the interval fi(Q 2 ) < 77 < r, so eq. ( 3.26 ) involves an 
integral over this interval. 

Of course, eq. ( |3.26 ) can be also obtained directly from the BK equation, although, in that 
context, its physical interpretation in terms of dipole-dipole scattering may not be so obvious. 
In fact, this is just a particular piece of what we have called "the contribution (B)" in the 
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previous subsection, i.e., the contribution of the points z± satisfying \z± — b±\ 3> 1/Q s (t, b±). 
If z± now refers to the saturated gluons within the black disk, then it is further restricted by 
\z±\ < R(t,Q 2 ), which implies \z± — b±\ ~ b for the same reasons as above. Then, a simple 
calculation similar to eq. fl3.8|) immediately yields 



T 2 

-^S T {x ± ,y ± ) ~ -a s ^R\r,Q 2 )S T (x ± ,y ± ), (3.27) 

which after integration over r is indeed equivalent to eq. fl3-26| ) PI 

For comparison with eq. ( |3.26| ), it is interesting to compute also the S'-matrix element for 
a coloured external probe, e.g., a quark, which scatters off the saturated gluons in the eikonal 
approximation. A calculation entirely similar to that leading to eq. Q3.26 ) yields (b± is the 
transverse location of the quark) : 

1 /. T rt„ A f 9 2 C F 



Wc (trV\b ± )) r = exp l-i±*L £ dr, 7fl (b ± ,b ± ) 



exp j-a s ±- / d n R-U,.Q-)\ (3.28) 



where the second line follows after using eq. ( |3.24 ) 



To summarize, the amplitude for the scattering off the saturated gluons decreases like 1/6* 
for an external dipole, but only as 1/6^ for a coloured probe. This difference turns out to be 
essential: because of it, this long-range scattering plays only a marginal role for the dipole, while 
it leads to unitarity violations in the case of the coloured probe (although the very question of 
unitarization makes little physical sense for a "probe" which is not a colour singlet). 

To see this, assume the long-range contributions shown above to be the only contributions, 
or, in any case, those which give the dominant contribution to the cross-section. Then, one 
can rely on the previous formulae to estimate the rate of expansion of the black disk. Namely, 
assume that, for the purpose of getting an order-of-magnitude estimate, one can extrapolate 
eqs. ( ^26|) and (|3|) up to energies where the black disk approaches the incidence point b± of 
the external probe. Then, we expect the exponents in these equations to become of order one 
for b ~ R(t,Q 2 ). For the external dipole, this condition implies: 

This gives (recall that Q 2 = 1/rjJ: 

Q 4 R\t, Q 2 ) = a s J T _ d<n Q 2 R 2 (v, Q 2 ), (3.30) 
or, after taking a derivative w.r.t. r, 

2^(Q 2 R 2 (t,Q 2 )) = a s , (3.31) 



9 The mismatch by a factor of two betw een eqs. (3.26) and (3.27) is inherent to the mean field approximation 
used in Refs. |l8j to derive eq. (3.25), and is completely irrelevant for the kind of estimates that we are 
currently interested in. 
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whose solution R 2 (t,Q 2 ) increases linearly with r. 

By contrast, for a coloured probe, the same condition yields: 



a s j_ d V R 2 (7],Q 2 ) = R 2 (t,Q 2 ) (3.32) 
or after taking a derivative w.r.t. r: 

^-R 2 (t,Q 2 ) = a s R 2 (r,Q 2 ), (3.33) 



which gives an exponential increase with r, as found in Ref. [14]. 



We thus see that the violation of unitarity by long-range Coulomb scattering reported by 



the authors of Ref. [14| is related to their use of an external probe which carries a non-zero 
colour charge. This case is physically ill defined, and therefore uninteresting (note, indeed, 
that (trW) is not a gauge-invariant quantity); in particular, its relevance for the problem of 
gluon saturation in the target wavefunction remains unclear to us (since the relation between 
"blackness" and saturation holds only for dipole probes; cf. the discussion prior to eq. ( |2.9| )). 



On the other hand, for the physically interesting case of an external dipole, the contribution 



(3.31) to the expansion of the black disk not only is consistent with unitarity — if this was the 
only contribution, the cross-section oc R 2 (t, Q 2 ) would increase linearly with r — , but at large 
r, is even negligible as compared to the corresponding contribution of the short-range scattering 
(which gives a cross-section increasing like r 2 , cf. eqs. ( |2.16|) -( ^17 )). 



These considerations are conveniently summarized in the following, schematic, approxima- 
tion to the scattering amplitude in the grey area, which follows from the previous analysis in this 
section: Af T (r± , bj_ ) is the sum of two contributions, a short-range contribution, cf. eq. ( 2.15| ), 
and a long-range one, cf. eq. ( 3.26| ) : 



M T (r ± ,b ± ) « v /^A2>^e- 2m - 6 + a s ^ f dq i? 2 (r?,Q 2 ), (3.34) 

with the short-range contribution determined by the solution to the homogeneous BFKL equa- 
tion ( |3.15| ) together with the assumed exponential fall-off of the initial condition (see Sects. 4 and 
5 below for more details), and the long-range contribution obtained by keeping only the lowest- 
order term in eq. ( 3.26|) (which is enough since we are in a regime where S T (x±,y±) ~ 1). At 



the initial rapidity To, the long-range contribution vanishes, while the short-range contribution 
reduces to e -^(b-R{r ,Q 2 )) with #( r0) Q2) given by eq (|1|) ) as i t should^. 

For given r and r± = 1/Q, eq. ( |3.34 ) applies at impact parameters in the grey area, 



R(t,Q ) < b < Rh(t), but it can be extrapolated to estimate the boundaries of this area, 
according to eqs. (|2.9|)-( |2"?L2| ). It is easy to check that, at high energy, both these boundaries are 



determined by the short-range contribution, i.e., the first term in the r.h.s. of eq. ( 3.34 ). Thus, 
this contribution dominates the scattering amplitude at any b± in the grey area. By comparison, 
the long-range contribution is suppressed by one power of l/a s r. 

For instance, the black disk radius is obtained by requiring M T (r±, b±) ~ 1 for b± ~ R(t, Q 2 ) 
(cf. eq. (|2.9| )). If one assumes the short-range contribution to dominate in this regime, one 



That is, the initial scattering amplitude M T n (r± , b±) is equal to one within the black disk (b±_ < R(ro,Q 2 )), 



and decreases exponentially outside it. Eq. (3.34) applies, of course, only at impact parameters outside the black 
disk. 
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obtains the estimate ( |2.16| ) for the black disk radius: R(t,Q 2 ) ~ (w/2m ir )a s r. By using this 
result, one can evaluate the corresponding long-range contribution, and thus check that this is 
comparatively small, as it should for consistency with the original assumption: 



drjR 2 (r,,Q 2 ) 



1 1 

7, rat 



< 1. 



(3.35) 



By contrast, if one starts by assuming that the long-range contribution dominates, then one 
is running into a contradiction, since in this case R(t,Q 2 ) oc y/a s r, cf. eq. ( |3,31 ), and the 
short-range contribution increases exponentially along the "trajectory" b = R{r,Q 2 ). 

A similar conclusion holds for b = Rjj{t) [since Rh{t) = R(t,Q 2 = Aq CD ), cf. eq. ([2. 121) ], 
and therefore for any point b± within the grey area. (The hadron radius Rh{t) will be evaluated 
in Sect. 4.2 below.) Thus, the short-range contribution is indeed the dominant one in the grey 
area. This contribution preserves the exponential fall-off of the initial condition, and therefore 
saturates the Froissart bound, as explained in Sect. 2. 

It is essential for the consistency of the previous arguments — which combine perturbative 
quantum evolution with non-perturbative initial conditions — that perturbation theory has 
been applied only in the regime where it is expected to be valid, namely, in the central region 
at b < Rh{t), where the gluon density is high and the local saturation scale Q s (r,b) is much 
larger than Aqcd- It has been enough to consider this region for the present purposes since this 
includes both the black disk and the grey area which controls its expansion. Within this region, 
the perturbative evolution equations of Refs. @, [8|, |9|, [n| can be trusted, and the additional 
approximations that we have performed on these equations are under control as well. When 
supplemented with appropriate boundary conditions — which are truly non-perturbative, since 
reflecting the physics of confinement — , these equations allow one to compute the rate for the 
expansion of the black disk, and, more generally, to follow this expansion as long as the black 
disk remains confined within the region of applicability of perturbation theory (which includes, 
at least, the central area at the initial rapidity tq : b < Rh(to)). 

But if one attempts to follow this expansion up to much higher rapidities, where R(t, Q 2 ) S> 
Rh(tq), then a strict application of the perturbative evolution with initial conditions at r = To 
would run into difficulties! 1 " 1 ]. The difficulties arise since, in the absence of confinement, the 
long-range dipolar tails created by the saturated gluons can extend to arbitrary large distances, 
and thus contribute to scattering even at very large impact parameters (b 3> Rh(to)), where 
physically there should be no contribution at all. This is illustrated by eq. ( 3.34[ ): We have 
previously argued that, for b± within the grey area, the dominant contribution comes from short- 
range scattering (the first term in the r.h.s. of eq. ( 3.34| ) ) . But if one extrapolates this formula 
at very large b 3> Rh(t), then, clearly, the long-range contribution, which has only a power-law 
fall-off with b±, will eventually dominate over the short-range contribution, which decreases 
exponentially. When this happens, however, the impact parameters are so large (b — R(r, Q 2 ) > 
l/uin) that the long-range scattering is controlled by the exchange of very soft (k± m^) 
quanta, which in a full theory would be suppressed by the confinement. That is, in a more 
complete theory which would include the physics of confinement, the long-range contribution 



n The discussion in this paragraph has been inserted as a partial response to criticism by Kovner and Wiedemann 
p2[ , written in response to the original version of this paper. Please see the note added at the end of the paper 
for further discussion. 
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to eq. (3.34) would be suppressed at very large b± by an additional factor e 2rn *( b R ( T >Q )) ? 



so that the short-range contribution will always dominate, for all impact parameters. But in 
the present, perturbative, setting, the only way to avoid unphysical long-range contributions is 
to start the quantum evolution directly in the grey area, as we did before (rather than try to 
construct this grey area via perturbative evolution from earlier rapidities t -C t, at which the 
points b± of interest were far outside the initial grey area: b S> Rh(j~o))- 

Note finally that what is truly remarkable, and also essential for our conclusion on the 
Froissart bound, is not the suppression of the long-range non- perturbative contributions by the 
confinement — this is only to be expected in the full theory, and can be also enforced in 
the present calculation by appropriately chosing the boundary conditions — , but rather the 
suppression of the long-range perturbative contribution within the grey area (where perturbation 
theory applies, so its predictions must be taken at face value). We mean here, of course, the 
fact that the long-range contribution to eq. ( |3.34 ) falls off like and not like l/b\, as it 



would have been the case if the saturated gluons were uncorrelated. The mechanism for this 
suppression is purely perturbative, and related to saturation: The saturated gluons are globally 
colour neutral, so the monopole fields of the individual gluons are replaced at large distances 
3> 1/Qs by the more rapidly decreasing dipolar field of the whole distribution. To understand 
the relevance of this suppression for the Froissart bound, consider what would happen if the 
saturated sources were statistically independent, i.e., if eq. ( 3.19| ) was replaced by 



(p a v (z ± )p b v ,(u ± )) T = X5 ab 5(r ] - V ')5^(z ± -u ± ). 



Then the exponent in eq. ( |3.26| ), and also the second term in the r.h.s. of eq. ( |3.34| ), would 
change into 

&s £ drjR 2 (r),Q 2 ) , 



which would generate a black disk increasing exponentially with r (cf. eq. ( 3.32 )-( 3.33 )). Thus, 
the long-range contribution would dominate already within the grey area, and the Froissart 
bound would be violated. We thus conclude that colour correlations at saturation are essential 
to ensure unitarity. 



4 Black disk evolution and the Froissart bound 



In this section, we shall exploit the factorization property ( 3.1 6| ) together with the known solution 



■A/t ) to the homogeneous BFKL equation in order to compute the scattering amplitude in the 
grey area, and thus study the evolution of the black disk with increasing energy. After briefly 
recalling the BFKL solution, in Sect. 4.1, we shall then compute the radius of the black disk 
R(t,Q 2 ) and derive the Froissart bound (in Sect. 4.2). Then, in Sect. 4.3, we shall study the 
impact parameter dependence of the saturation scale Ql{T,b±), and deduce a physical picture 
for the expansion of the black disk, to be exposed in Sect. 4.4. 

4.1 Scattering amplitude in the BFKL approximation 



Eq. ( 3.16| ) for the scattering amplitude in the grey area involves the solution M T (r±) to the 



homogeneous BFKL equation, i.e., the BFKL equation without impact parameter dependence. 
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This solution is well known, and we shall briefly recall here the relevant formulae, at the level of 
accuracy of the present calculation. (See Refs. |37], ^] for a similar approach and more details.) 
The solution can be expressed as a Mellin transform with respect to the transverse coordinate: 

K(r± = l/Q) = f — ( ^) e^Wm-^W-iKi-X)} , (4.1) 



where ip(\) is the di-gamma function, and A is an arbitrary reference scale, of order Aqcd- 
The contour C in the inverse Mellin transform is taken on the left of all the singularities of the 
integrand in the half plane Re A > 0. Note that, since M T (r±) is a function of r"j_, we find it 
convenient to use the momentum variable Q 2 = l/r\ to characterize the transverse resolution 
of the dipole. From now on, we shall again use the notation M T {Q 2 ) = M T (r± = l/Q), which 
was already introduced in Sect. 2 (cf. eq. Q2.9| )). 

We are interested here in a regime where the energy is very high, 6l s t 3> 1, and the dipole is 
small: Q 2 3> A 2 . In these conditions, it is appropriate to evaluate the integral ( |4.1| ) in the saddle 
point approximation. Higher is the energy, better is justified this approximation, and closer is 
the saddle point Ao — which is a function of (In Q 2 / A 2 ) ja s r — of the so-called "genuine BFKL" 
saddle-point at Ao = 1/2. (This is the saddle point which governs the asymptotic behaviour of 
the solution to the BFKL equation at very large energy.) In fact, for 

^lng«l, (42) 
which is the most interesting regime here, the saddle point is easily estimated as: 



with (5 = 28£(3). In fact, the recent analysis in Ref. |29[ shows that eq. (|4.3| ) remains a 
good approximation for the saddle point even for comparatively low energies, such that a s r ~ 
ln(Q 2 /A 2 ). This saddle point gives the standard BFKL solution, which, after multiplication 
with the profile function (cf. eq. ( 3.1 6| ) ) , provides the scattering amplitude in the grey area in 
the present approximation: 

/ A 2 „u« s t 




where u = 4 In 2 is the customary BFKL exponent. The factor \J2tt(56l s t in the denominator 
comes from integrating over the Gaussian fluctuations around the saddle point. When expo- 
nentiated, this gives a contribution cx ln(a s r) which is subleading at large energy and will be 
ignored in what follows. It is then convenient to rewrite eq. (|4,4j) as follows: 



1 Q 2 if Q 2 \ 2 \ 
N T {Q Z , b±) ~ exp { -2m n b + ua s T - - In - — ^ I In -^j 1 \, (4.5) 

where we have also used S(b±) e - 2m * b ^ as appropriate for sufficiently large b (b 3> i?o, cf. the 
discussion after eq. fl2.6|) ). This is the most interesting case here, since we consider the high 
energy regime in which the black disk is already quite large. 
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Eq. fl4.5| ) is valid for those values of the parameters r, Q 2 and b± for which our previous 
approximations are justified, namely, such that the conditions A 2 <C Q 2 (t, b) <C Q 2 are satisfied. 
As it was anticipated in Sect. 2, and will be verified below in this section, these conditions 
are realized within a corona at R(t,Q 2 ) <C b± <C Rh(t), which, with increasing energy, moves 
further and further away from the center of the hadron. 

When decreasing b± towards R(t,Q 2 ) at fixed r, or, equivalently, increasing r at fixed b±, 
the scattering amplitude increases towards one, and the BFKL approximation (|4.5| ) ceases to be 
valid. (The dipole resolution Q 2 is always fixed in these considerations.) But it is nevertheless 
legitimate to use eq. (|4.5|) in order to estimate the boundary of its range of validity, that is, 
the black disk radius R(t,Q 2 ), or the saturation scale Q 2 (r,b±). Indeed, the non-linear effects 



become important when the BFKL solution (4.5) becomes of order one. This condition can 



be written either as an equation for R(t,Q 2 ) for given r and Q 2 , namely, eq. ( |2.9D , or as an 
equation for Q 2 (t, b±) for given r and b±, namely, eq. (|2 . 1 1 ) . (One could, of course, similarly 



introduce and compute also a critical rapidity t(Q 2 , b±) at which blackness is reached for given 
Q 2 and b±, but this is less interesting for our subsequent discussion. See, however, Sect. 4.4.) 

4.2 The black disk radius 

In this subsection, we shall use eqs. (|2.9|) and (|4.5|) to compute the radius of the black disk and 
study some limiting cases. Eq. ( ^9| ) amounts to the condition that the exponent in eq. 
vanishes^, which immediately implies: 

2m n R(r, Q 2 ) = ua s r - i In g - ^ |!) . (4.6) 

The right hand side is positive as long as Q 2 < A 2 e cctsT = Q 2 (t, 6 = 0), with 
P , 1 



c=-^ + ~^P(P + 8u) = 4.84.... (4.7) 

As anticipated by our notations, 

Q 2 {r,b = 0) = A 2 e c " sT (4.8) 

is the saturation scale at the center of the hadron (this will be verified via a direct computation 
in the next subsection). This is as expected: for Q 2 > Q 2 (r,b = 0), the hadron looks grey 
everywhere, so i?(r, Q 2 ) = 0. 

The other extreme situation is when Q 2 ~ A 2 , so that the black disk extends up to the edge 
of the hadron (cf. eq. ( 2.12] )). Equation (|4.6| ) yields then: 



«-W-^r, (4.9) 

which should be seen only as a crude estimate: for such a small Q 2 , our approximations are not 
justified any longer. 



J At the level of the present approximation, one can take n = 1 in eqs. (2.9) and (2.11) without loss of accuracy. 



27 



But the physically interesting case is when Q 2 S> A 2 , but the energy is so large that the 



condition (4.2) is satisfied. Then one can neglect the term quadratic in lnQ 2 /A 2 in eq. 



(since this term vanishes when r — > oo), and deduce that: 



^'^^-^J- (410) 

The term linear in lnQ 2 /A 2 , although subleading at large r (since independent of r), has been 
nevertheless kept in the above equation since, first, we expect this term to give the dominant 
Q 2 -dependence of the cross-section at high energy, and, second, it measures the separation 
between the black disk and the edge of the hadron in the high energy regime. Specifically: 

R H (T)-R(T,Q 2 )K-±-\n^, (4.11) 

which is fixed (i.e., independent of r), but large when Q 2 ^> A 2 . This is important since 
the points b± at which our approximations are justified should lie deeply within this corona: 
R(t,Q 2 ) < i» < Rjj(t). Thus, as anticipated in Sect. 2, the fact that Q 2 S> A 2 ensures the 
existence of a large grey area in which our approximations apply. 

Equation ( 4.10| ) is our main result in this paper. It shows that, at very high energy, the 



radius of the black disk increases only linearly with r, i.e., logarithmically with the energy. This 
is the result anticipated in eq. ( |2.16| ). The corresponding cross-section is given by eq. ( [2.171 ) 
and saturates the Froissart bound, that is, it grows like In 2 s, with a proportionality coefficient 
which is universal (i.e., the same for any hadronic target), and which reflects the combined role 
of perturbative and non-perturbative physics in controlling the asymptotic behaviour at high 
energy. 

In the remaining part of this paper, we shall further explore this result and gain a different 
perspective over it by computing also the saturation scale and studying the geometric scaling 
properties. 

4.3 Saturation scale with the impact parameter dependence 

In this subsection, we shall compute the saturation scale for an inhomogeneous hadron and 
study its variation with the energy and the impact parameter. Previous studies of this kind 
were restricted to a homogeneous hadron [21, 2S, 3C], but, as we shall see, the dependence upon 
the impact parameter introduces some interesting new features. 

By inspection of eq. (|4.5|), it is clear that the saturation condition ( |2.11| ) amounts to the 
following, second order algebraic equation for the quantity p s = (l/a s r) In Q 2 /A 2 : 

p 2 + Pp s -2(3u; = -2P^. (4.12) 

The solution to this equation and the corresponding saturation scale read: 




Ps(r,b) = + + — 1-— - ), (4.13) 
2 2 y p \ tua s T J 

Q 2 (r,b) = A 2 e « sTps( - T > b \ (4.14) 
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Figure 3: The saturation scale Q s {b)/Q s {b 
Woods-Saxon profile function of eqs. (|2,6|) -(| 
distance is measured in units of 1/2711,7,-. 



0) from eqs. (4.13)-(4.14) for a s T = 3 and the 
with Rq = 3/2m n . On the abscisa, the radial 



Note that, in general, the impact parameter dependence in the saturation scale ( [4.14 ) is not 
factorizable. Below, however, we shall recover factorization in some specific limits. 

It can be easily checked that the above equations ( 4.13| )— ( 4.14 ) and eq. ( [4.6D are consistent 
with each other, in the sense that Q 2 {r,b = R{t,Q 2 )) = Q 2 , as it should (cf. eq. (^10|)). 
In particular, one can use eqs. ( |4.13| )-( 4.14| ) to rederive the results in eqs. (f4.7|)-(f4.8|) for the 




for the hadron 



saturation scale Q s {r,b = 0) at the center of the hadron, as well as eq. 
radius. 

A pictorial representation of the 6-dependence of the saturation scale, as emerging from 
eqs. ( 4.13[ )-(4.14), is given in Fig. |3[ As compared to eq. (4.13), in this graphical representation 
we have replaced 2m- w b — ► — In S(b), with S(b) given by a Woods-Saxon profile, cf. eqs. ( |2.G| ) 
(|2.7| ); this is more realistic than the exponential at short distances, b ^ Rq, where it has a much 
slower decrease, but it shows the expected fall-off S(b) ~ e -2m 7r (b-R ) a ^ l ar g 6 r distances. As 
manifest on this figure, Q 2 (t, b) is itself very similar to an exponential for all distances 6 ^ Rq. 
This can be understood via a further study of eq. ( 4.13| ), which will also reveal that, in fact, there 
is a change in the slope of the exponential with increasing b : To a very good approximation, 
the plot in Fig. | can be seen as the superposition of two exponentials, one at small b, the other 
one at large 6, which have different exponential slopes. 

To see this, note that the function 



A(r,6) 



1 



2171^5 



1 



(4.15) 



Loa s r Rh{t) 

which enters the square root in eq. ( [4.13|) , is positive semi-definite for b within the hadron radius 
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(b < Rh(t)), and monotonically decreasing with b from A(r, b = 0) = 1 to A(r, 6 = Rh{t)) = 0. 
This suggest two different approximations according to whether A is close to one (for b sufficiently 
small) or close to zero (for b sufficiently close to Rh(t)). (Note that the factor 8to/P multiplying 
A(t, b) in eq. ( 4.13 ) is a number of order one, 8u/f3 ~ 0.67, so it does not interfere with our 
order-of-magnitude estimates.) 

(I) If A is close to one, which happens when b is much smaller than the hadron radius: 



2m n b 



«1, or b < Rh(t), 



(4.16) 



one can evaluate the square root in eq. fl4.13| ) in an expansion in powers of 1 — A (this is 
equivalent to an expansion of p s (r,b) in powers of b around p s (T,b = 0) = c, cf. eq. (|4.7|) ). To 
linear order in this expansion, one obtains: 



p s (r,b) ~ c 



2m n b 



X, 



Ij3 + Suj 



4/? 



X s a s T 

which gives (with 7 = 1/A S ~ 1.55): 

Q 2 s {t, b) ~ A 2 e c ^ T e ~ 2 ^ b = Q 2 (r, b = 0) [5(6)] 7 



0.644..., 



(4.17) 



(4.18) 



This is in a factorized form, although, as compared to the corresponding factorized structure 
of the scattering amplitude ( |3.16| ), it features some "anomalous dimension" 7 for the profile 
function. The value Q 2 (t, b = 0) at the center of the hadron is the same as the saturation scale 
for a homogeneous hadron previously found in Refs. |21, 2£, 3C]. Also, the constant A s which 



appears in eq. ( 4.17 ) is the value of the saddle point Ao in the Mellin representation (4.1) for 
Q 2 = Q 2 s (t, b = 0) (i.e., eq. Q with ln(Q 2 /A 2 ) = ca s r) |§. 

Equation (|4.1 8] ) shows that the saturation scale decreases exponentially with the distance b 
from the center of the hadron, with a typical decay scale ~ 1/27771^. (Of course, this exponential 
law applies only for values b which are not too close to the center, b <; Rq, cf. Fig. ||.) 
(II) When b is sufficiently close to Rh(t), in the sense that: 



A(r,6) 



Rh(t) 



Rh(t) 



« 1, 



(4.19) 



than one can expand eq. ( |4.13j ) in powers of A (this is an expansion of p s (r, b) around p s (t, b 
Rh) = 0). To lowest order in this expansion, one obtains p s (r,b) ~ 2o;A(r, 6), and therefore: 



Q 2 s (r,b) 



^2 e 2aja s T ^-Am^b 



(4.20) 



Thus, the saturation scale in the tail of the hadron distribution is still in a factorized form, 
but the exponential slopes are different as compared to the corresponding form near the center, 
eq. flOp , both for the increase with r — which is now controlled by the BFKL exponent 
2lo ~ 5.55 > c — and for the decrease with b — where the "anomalous dimension" 7 of 



eq. ( 4. 18|) has been now replaced by 2. 

These changes can be easily understood by reference to eq. 
vanish when Q 2 



: 5|): The exponent there must 



Qj(r,b). If b satisfies the condition ( 4.19 ), then the first two terms in the 



exponent, which are the large terms, almost cancel each other, so the other terms there must 
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be relatively small, in the sense of eq. Q4.2| ). Then, the term quadratic in ln(<5 2 /A 2 ) is much 
smaller than the linear term, and can be neglected. We thus end up with 

f 1 A 2 1 

M T {Q\ b ± ) « ex P j-2m^ + W a s r--ln^l = J_ e ^- 2 ^ b , (4.21) 



which, together with the saturation criterion ( |2.11 ), provides indeed the expression (4.20) for 



the saturation scale. To summarize, when the "diffusion" term in the BFKL solution (|4.4j ) 
becomes negligible, then the energy dependence and the Q 2 -dependence of the solution are fully 
controlled by the "genuine" BFKL saddle-point at Ao = 1/2. 

The change of behaviour from eq. ( |4 . 1 8| ) to eq. ( 4.2C| ) is also visible on a logarithmic plot 



of the saturation scale in eqs. ( f!.13| )-( f4.14| ) as a function of b. In Fig. ^ we have displayed 



the function p s (r, b) of eq. ( 4.13|) as a function of b/Rjf(r), together with its small-distance and 



long-distance approximations, as given by eq. ( f4.17| ), and prior to eq. ( 4.20| ), respectively. As 



explicit on this figure, the transition between the two regimes is rather smooth, and takes place 
at intermediate values b ~ Rh(t)/2. 

4.4 Expansion of the black disk 

Let us finally consider the implications of the previous results on Q s for the expansion of the 
black disk. This is interesting since, as we shall discover in Sect. 5 below, the two domains 
(I) and (II) are characterized by different "geometric scaling" laws for the black disk radius 
R(t,Q 2 ), and thus for the dipole cross-section. 
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Equation ( 4.18|) together with the definition ( p.lO| ) of R(t,Q 2 ) imply (recall that 7 = 1/A S ): 



2m 7T R(T, Q 2 



x ( - , & 



(4.22) 



This is valid as long as R(t,Q 2 ) satisfies the condition ( 4.16| ), in practice, for R(t,Q 2 ) ^ 
Rh(t)/2. For a given Q 2 3> A 2 , this happens only within an intermediate range of energies, to 
be specified shortly. Of course, eq. (|4.22j ) is just an approximate form of the general expression 
(|4.6| ) valid in this intermediate range of energies, but the approximations necessary to derive 
eq. ( |4.22| ) may not be easily recognized at the level of eq. (4.6). These approximations will be 
clarified in Sect. 5 below. 

But at sufficiently high energy, the border of the black disk lies in domain (II), that is, it is 
relatively close to the edge of the hadron, in the sense of eq. ( f4.19| ). To verify this, one can use 
eq. Q4.ll ) to deduce that the ratio: 



R h (t)-R(t, Q 2 
Rh(t) 



1 



In 



Q 2 



2ua,r A 2 



(4.23) 



is decreasing with r, and therefore necessarily satisfies the condition ( 4.1S| ) for sufficiently large 
r. Also, it can be easily checked that the expressions ( 4.2CQ for the saturation scale in domain 
(II) and ( |4.10| ) for the black disk radius at high energy are consistent with each other, via 
eq. ( |2.1C| ) . This is as it should since both expressions have been obtained via the same high energy 
approximation, namely, they follow from the asymptotic form ( |4.21 ) of the BFKL solution. 

The evolution of the black disk with increasing r at fixed Q 2 is pictorially illustrated in Fig. 
||. The black disk appears at the center of the hadron at a critical rapidity f\{Q 2 ) such that 
Q 2 (fi,b = 0) = Q 2 . This condition together with eq. (^4.8|) implies: 



n(Q 2 ) = — 



ca. 



1 Q2 



(4.24) 



For t > fi(Q 2 ), but not much larger, the black disk remains confined to domain (I), as illustrated 
by the smallest disk on the left of Fig. ||. However, the expansion rate of the black disk, which 
is equal to cA s ~ 3.12 in appropiate units (cf. eq. ( |4.22| )), is faster than the corresponding rate 
lo/2 ~ 1.39 for the borderline b{r) = Rh(t)/2 between the two domains. Thus, at some new, 
larger, critical value, which can be easily estimated from eq. ( |4,22| ) as 



HQ 2 



1 



1 Q 2 

ln A7 



(4.25) 



(c - uj/2X s )a s 

the black disk reaches domain (II), and than extends further within this domain (as illustrated 
by the two larger disks in Fig. ||). But within domain (II), the expansion rate of the black 
disk slows down to uj ~ 2.77, cf. eq. ( 4.1C| ), which is the same rate as for the hadron outer 
border Rh(t). Thus, with increasing r, the radial distance between the border of the black disk 
and the edge of the hadron (i.e., the width of the grey area) remains constant (cf. eq. ( (4.11| )), 
so the relative size of the grey area with respect to the hadron size is smaller and smaller (cf. 
eq. ( f4.23| )). At large r, this grey area represents the tail of the hadron wavefunction, in which 
the BFKL equation remains valid even at arbitrarily large energy. 

Given the (quasi)exponential decrease of the saturation scale Q 2 (r,b) with b, it is straight- 
forward to verify that the conditions A 2 <C Q 2 (r, b) <C Q 2 hold at any point b± in the corona 
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higher energy 



(II)/' 



(I) 



T"<T ; 




V <X 

Figure 5: A pictorial representation of the expansion of the black disk with increasing rapidity. 
The dotted line circle of radius b{r) = Rjj(t)/2 separates between domains (I) and (II). 



R(t,Q 2 ) <C b± <C Rh(t), as necessary for the consistency of our approximations. Indeed, these 
conditions are satisfied as soon as the separation between b± and the (inner or outer) edges of 



this corona is of order l/m n or larger. For instance, eq. (4.20) can be rewritten as: 

Q 2 s (T,b) ~ A2 e 4m,(fiH(r)-b)_ ( 42 g) 

Similarly, by using eq. ( PIOD , the expressions ( |4.18| ) and ( |4.20 |) can be recast into the form: 



2 -2 7 m T (fe-_R(r,Q 2 )) 



Qi(r,b) ~ Q 2 e 

where 7 ~ 1.55 if b± is in domain (I), and 7 = 2 for b± in domain (II). 



(4.27) 



5 Geometric scaling at high energy 

"Geometric scaling" refers to the property of the dipole-hadron total cross section <t(t, Q 2 ), 
eq. (|2.5|), to depend upon the two kinematical variables r and Q 2 only via the combination 
Qq(t)/Q 2 (the "scaling variable") where Qq(t) is some suitable momentum scale which increases 
as a power of the energy: Qq(t) cx e Ar . This property is interesting since it can be related to a 
similar property of the virtual photon total cross section cr 7 * p (x, Q 2 ) which is actually seen in 
the HERA data on deep inelastic scattering (for x < 0.01 and Q 2 < 400 GeV 2 ) |§. 

Clearly, this property cannot hold for arbitrary r and Q 2 , since it is known to be violated 
by the solutions to the linear evolution equations (DGLAP or BFKL) at very high Q 2 . On 
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the other hand, in the saturation regime at Q 2 < Q 2 (t), this property is physically motivated, 
since the saturation scale is then the only scale in the problem |2(J (at least for a homogeneous 
hadron), so it is tempting to identify it with the momentum scale Qq(t) introduced above. 

So far, studies of geometric scaling have been performed only for a homogeneous hadron, so 
they have naturally focused on the corresponding property of the scattering amplitude M T {Q 2 ) 
(because a(r,Q 2 ) = 2irR 2 N T (Q 2 ) in this case). Such previous studies — which relied on either 
numerical [39|, or (approximate) analytic p9| , |30| solutions to the BK equation — not 

only confirmed the existence of geometric scaling in the saturation regime, but also showed that 
this property extends up to momenta Q 2 considerably larger than Q 2 (t). In particular, in Ref. 
p9| , the upper limit for "extended scaling" has been estimated as Q 2 max ~ Qs(t)/A 2 , which is 
roughly consistent with the phenomenology |3^] (see also below). 

Our purpose in this section is to extend the analysis of Ref. [^9| by taking into account the 
transverse inhomogeneity in the hadron, within the formalism developed in the previous sections 
of this paper. Note that, since the impact parameter is integrated over in the formula (^5|) for 
<t(t, Q 2 ), the connection between the scaling properties of the cross-section and those of the 
scattering amplitude M T (r±, b±) becomes more subtle now. In particular, it is not a priori clear 
what should be the scale which plays the role of Qq(t) in the scaling variable, or even whether 
geometric scaling exists at all (since the inhomogeneous problem is not a single-scale problem 
any longer). 

To progressively introduce the effects of the inhomogeneity, let us start with a situation where 
the rapidity r is small enough for the hadron to look "grey" everywhere: Q 2 > Q 2 (t, 6 = 0), or 
r < fi(Q 2 ), cf. eq. (|4.24| ). In that case, and within the present approximations, the scattering 
amplitude M T (Q 2 ,b±) is factorized at any b± : M T {Q 2 , b± ) ~ ftf T (Q 2 )S(b±). Thus, the study of 
the scaling properties of the total cross-section: 

a(r,Q 2 ) ~ 2Af T (Q 2 ) J d 2 b ± S(b ± ) (all grey), (5.1) 



reduces to the corresponding study of M T {Q 2 ) in the homogeneous case p9| . It is instructive to 
briefly rederive here the main result in Ref. [^9| (in a very schematic way): 

The homogeneous scattering amplitude Af T (Q 2 ) is given by (|4.5|) with 6 = 0. As such, this 
formula shows no scaling. To see the scaling emerging, we shall replace the arbitrary reference 
scale A 2 in this equation with the saturation scale Q 2 (t) = A 2 e CQsT (which is the same as the 
central saturation scale Q 2 (r, 6 = 0) in the inhomogeneous case; cf. eq. (f4.8[)). We write: 



, Q 2 , Q 2 , Q 2 s (r) . Q 2 , _ 

A2 = Q|(r) ~h?~ = ln QW) + CasT - ( } 



Since, by construction, the saturation scale is such that the exponent in ( [4.5D (with 6 = 0) 
vanishes for Q 2 = Q 2 (t), it is clear that, after the replacement (5^), we are left only with terms 
involving, at least, one power of ln(Q 2 /Q 2 (t)) : 

K(C?) = «P< -A, b - ^ (in^) ) ■ (5.3) 



Here, A s rj 0.64 has been generated as X s = 1/2 + c/f3 (cf. eqs. ((Oj) and ( 4.17| )). So far, (|5 



is just a rewriting of eq. ([4 . 5[) . But this is suggestive of the conditions under which geometric 
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scaling should be expected: This emerges when the (central) saturation scale Q 2 (t) is sufficiently 
close to Q 2 (although still smaller than it) for the quadratic term in eq. (|5,3|) to be negligible 
compared to the linear term. When this happens, eq. (|5.3|) can be approximated as: 



K(Q 2 ) - i^) S (5-4) 

which shows geometric scaling indeed, with Qq(t) = Q 2 (t). This is valid as long as: 
O 2 

1 < ln 775~TT « 2A,/?a s r, (5.5) 

where the lower limit is simply the condition that we are in a regime where eq. (|Tl|) applies: the 
hadron looks "grey" everywhere. Since 2\ s f3 ~ 43.37 is a large number, eq. ( |5.5| ) gives a rather 
large window, which however extends beyond the validity range of the BFKL saddle point 
A more complete analysis shows that eq. (5J5) should be replaced by(^|: 



i^l Q2 ^ l Q * {t) ~ «\ 

1 < In Q 2 ( T ) < m A 2 = ca « T ' ( 5 - 6 ) 

which for the present purposes is rewritten as a range for r: 

ln TT < T < -r ln T?- ( 5 - 7 ) 



2ca s A 2 ca s A 2 ' 

The lower limit in the equation above, which arises from the upper limit in eq. ( |5.6| ), is the 
smallest value of r at which the BFKL solution starts to behave like a scaling function, for 
a given Q 2 . As for the upper limit — in which we recognize the critical value fi(Q 2 ) for the 
emergence of the black disk, eq. ( [4.24 ) — , this is necessary only for the validity of eq. ([O]) . For 



higher values of r, geometric scaling may still hold, but in order to see it, the calculation should 
be modified to account for the formation of the black disk. 

Specifically, for r > f\(Q 2 ), the cross-section (|2.5|) can be then decomposed into a "black" 
contribution plus a "grey" one, which are evaluated as (for R(t,Q 2 ) > l/m n ): 

cr(r, Q 2 ) ~ 2ttR 2 (t,Q 2 ) + 2M T (Q 2 ) J d 2 b ± S(b ± ) Q(b - R(t,Q 2 )) 

. 2,R 2 (r,Q 2 ) + 2 -^^l. (5.8) 



The second line is obtained from the first one after replacing S(b) ~ e m , and noticing 
that Af T (Q 2 )e~ 2rri7TR ( T 'Q ) = 1. Eq. ( |5.8| ) confirms that, when the energy is large enough for 
R(t,Q 2 ) > I/to^, the total cross-section is dominated by the black disk, as anticipated in 
Sect. 2. It further shows that the scaling properties of the cross section are determined by 
the corresponding properties of the radius of the black disk, which can be inferred from the 
discussion in Sect. 4. As in Sects. 4.3-4.4, we are led to distinguish between two regimes: 



13 Note that the strong inequality on the logarithm in eq. J5.5| ) has been replaced in eq. (5.6) by a normal 
inequality, which is equivalent to a strong inequality on the argument of the log. This is the condition Q 2 <C 
Q1(t)/A 2 alluded to before H. 
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I) If ti(Q 2 ) < r < t 2 (Q 2 ), cf. eqs. ( |4.24| )-( p~25|), t he black disk lies entirely inside domain 
(I), and the corresponding radius is given by eq, ( f4,22j ), which is now rewritten as0 : 



R(t, Q 2 



X s , Q 2 s (r,b = 0) 



{I) 



2771^ 



In- 



Q 2 



(5.9) 



This shows scaling, with the same scaling variable Q 2 (t)/Q 2 as in the "all grey" regime (J5T 

II) In the high-energy regime r > f2(Q 2 ), the edge of the black disk is in domain (II), so its 
radius is given by eq. ( 4.1 0| ) . This shows scaling too, but with a different scaling variable: 



R(t, Q 2 



1 ^QM 



2 n 2wa a T 



(II) 4m n "' Q 2 ' 
This gives the scaling law for the total cross section at very high energy: 



<t(t, Q 2 ) ~ 2ttR 2 (t, Q 2 ) 



7T 



2mi 



In 



Ql(r) 

Q 2 



(5.10) 



(5.11) 



To summarize, for fixed Q 2 and intermediate energies corresponding to the following range 

Q 2 



of rapidities (cf. eqs. (5/7) and ( 4.25 )): 



i . Q 2 

2ca a A 2 



< r < 



In 



(c - uj/2X s )a s A 2 



= HQ 1 



(5.12) 



the dipole-hadron cross-section exhibits geometric scaling, with the scale set by the central satu- 
ration scale Q 2 (t) = Q 2 (r,b = 0). Eq. (5.12) allows for a significant window since, numerically, 
2c ~ 9.68 while c — u/2X s ~ 2.69. On the other hand, at higher energies, corresponding to 
r > T2(Q 2 ), there is scaling again, but the relevant scale is rather the asymptotic scale Q^o(t). 

Of course, the scaling breaks down, strictly speaking, for r ~ f2(Q 2 ), i-e., in the transition 
regime where the black disk crosses from domain (I) to domain (II). But it so happens that, 
numerically, the difference between the exponents c ~ 4.84 and 2uj ~ 5.55 in the corresponding 
scaling variables (Q 2 (t) and, respectively, Q^ir)) is quite small, so the transition from one 
scaling law to the other takes places rather fast, as also manifest in the plot in Fig. |||. 

It should be also stressed that, with incresing r, the geometric scaling becomes less and 



less relevant, since the (5 2 -dependence of the cross-section ( 5.11 ) eventually becomes a sublead- 
ing effect: When the condition ( f4.2j ) is fulfilled, the leading-order term in eq. ( |5.11| ) is simply 
proportional to r 2 , with the scale set by the pion mass (cf. eq. ( 2.1 7| ) ) . 

For completness, let us conclude with a discussion of the scaling properties of the local 
scattering amplitude N T (Q 2 , b±). Within the grey area, the factorization property J\f T (Q 2 , b±) = 
S(b±)M T (Q 2 ) implies that any scaling property of the homogeneous solution M T (Q 2 ) transmits 
automatically to M T {Q 2 , b±), but with a scaling variable involving the local saturation scale 
Q 2 (r,b) : 



M T {Q 2 



QM 
Q 2 



K(Q 2 ,b±) 



Q 2 s (r,b) 

Q 2 



(5.13) 



14 Incidentally, this also shows that the way to deriv e eq . ( |l.22|) fr om t h e gen eral expression (16) is via manip- 
ulations similar to those leading from eq. (4.5) to eq. (5.4), cf. eqs. (5.2 )— (5.3). 
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with Qs(t, b) = (5o('r)e _2m7rb/ ' A . This is precisely the scaling solution (3.11) that we used in our 
arguments in Sect. |3.1| . In eq. ( 5.13j ), the momentum scale Qq(t) and the power A depend upon 
the rapidity r, in the expected way: When r is in the intermediate range ( |5.12| ), Qq(t) is the 
saturation scale at the center (J 2 (t, b = 0), and A = A s 0.64. At higher rapidities, r > f2(Q 2 ), 
one rather has Qq(t) = Q^r) and A = 1/2. 

In both cases, eq. ( |5.13| ) holds only in the grey area at b > R(t,Q 2 ). But inside the black 
disk, geometric scaling holds as well, and almost trivially, since there A/" T (Q 2 , ~ 1 and the 



deviation from one shows scaling too 



We thus conclude that, for any r above the 



lower limit in eq. (|5. 
in the hadron disk 



the scattering amplitude M T (Q 2 , b±) shows geometric scaling everywhere 



6 Conclusions and outlook 

In this paper, we have proposed a simple mechanism combining perturbative gluon saturation 
and non-perturbative boundary conditions which ensures the saturation of the Froissart bound 
in dipole-hadron scattering at high energy. Gluon saturation has been implemented via the 
non-linear evolution equations derived in perturbation theory in Refs. |7|, ||, ^, that we 
have used only for impact parameters in the central region of the hadron, where the gluon 
density is high and perturbation theory is applicable. On the basis of these equations, we have 
shown that, with increasing energy, the non- linear effects ensure not only the unitarization of 
the scattering amplitude at fixed impact parameter, but also the factorization of the impact 
parameter dependence of the scattering amplitude in the outer "grey area" , where the unitarity 
limit has not yet been reached. This factorization, together with the exponential fall-off of the 
non-perturbative initial condition at large distances, leads to a total cross-section which grows 
like In 2 s. The coefficient of this growth is universal, i.e., independent of the hadronic target, 
and has been computed here. 

Our analysis makes explicit the deep connection between unitarization effects and the forma- 
tion of the Colour Glass Condensate. The "black disk" within which the unitarity limit has been 
reached is precisely the region of the hadron where the gluons "seen" by the incoming dipole are 
saturated, whereas the "grey area" corresponds to a region of lower density, in which the BFKL 
evolution still applies. The transition from "grey" to "black", i.e., the expansion of the black 
disk, is described by the non-linear evolution equations for the Colour Glass Condensate. The 
colour correlations at saturation are essential for both factorization and unitarization: They en- 
sure colour neutrality for the saturated gluons, and thus suppress the non-unitary contributions 
due to the long-range Coulomb tails. 

We have given the first computation of the saturation scale for an inhomogeneous hadron 
and identified two factorization regimes at different impact parameters. This has interesting 
consequences for the geometric scaling properties of the total cross-section: Our analysis predicts 
two different scaling laws in different ranges of energy. 

For simplicity, our analysis has been carried out for an external probe which is a qq dipole 
of fixed transverse size. But its extension to deep inelastic scattering should be straightforward. 
To this aim, the dipole-hadron cross-section must be averaged over all transverse sizes of the qq 



pair, with a weight given by the light-cone wavefunction of the virtual photon [|24| , 25, 26] 



It is important to point out two limitations of our approach, which call for further studies. 
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First, it was essential for our approximation scheme that the dipole is small, r± <C 1/Aqqd- in 
particular, this restricts the applicability of our results to deep inelastic scattering at relatively 
large photon virtuality, Q 2 3> A.q CD . It would be most interesting to see if one can phenomeno- 
logically extend this analysis towards the strong coupling regime at Q 2 ~ Aqcd usm g the "soft 
pomeron" , and thus get at least an estimate for the coefficient of the In 2 s growth in that regime. 

Second, we have performed our analysis at fixed coupling. It is an interesting open question 
how our results would be modified by the running of the coupling. Strictly speaking, this running 
is one of the next-to-leading order effects which so far have not been systematically included 
in the non-linear evolution equations (see, however, (^0|). But the experience with the BFKL 
equation, for which the NLO corrections have been recently computed suggests that a 
significant part of these corrections could be indeed taken into account by including the running 
of the coupling. Then, the natural question is, what should be the scale at which the coupling 
must be evaluated. 

In a previous paper [29], where the saturation scale Q s (r) for a homogeneous hadron has been 
computed, we found it natural to evaluate the coupling at the saturation momentum. Indeed, in 
that problem, the rapidity r and the external momentum Q 2 (the resolution of the dipole) were 
increased simultaneously, in such a way to preserve the saturation condition Q 2 (t) = Q 2 . For 
that particular running, we have shown that the saturation momentum changes its parametric 
dependence upon rapidity with respect to the fixed coupling case, from e CQsT to e Kv/T+T ° (with 
constant k and To). This result has been subsequently confirmed in Ref. f30fl . 

On the other hand, in the present paper the physical situation is quite different: The external 
momentum Q 2 is fixed, and with increasig r, we simultaneously increase the impact parameter 
b, in such a way that the condition Q 2 (r,b) = Q 2 (which defines the edge of the black disk) 
remains satisfied. It seems therefore natural to evaluate the coupling at the fixed external scale 
Q 2 , in which case the fixed coupling expressions obtained in this paper would remain valid 
(after the trivial replacement a s — > a s (Q 2 )). A further argument in this sense is provided by the 
discussion in Sect. 3.1, which shows that the dominant scattering involves only nearby colour 
sources, within the area covered by the incoming dipole or slightly further away, so that the 
typical transferred momenta are of order Q 2 . 

A still different possibility for the running, which would be closer in spirit to Ref. [29], would 
be to evaluate the coupling at the local saturation scale Q 2 (r, b) . Although this scale and the 
external scale Q 2 are identified at the edge of the black disk, they are nevertheless different 
at the points b± which lie further away in the grey area (where Q 2 > Q 2 s {r,b)). Thus, this 
choice for the running would probably modify the current formulae, in such a way to provide a 
generalization of the results in Ref. (2t| to the case of an inhomogeneous hadron. 

Of course, the description above is extremely crude and exploratory, and a full analysis of 
the running coupling case is required before any strong conclusions can be drawn. 
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Note added 

After the submission of the original version of this paper, there appeared a preprint by Kovner 
and Wiedemann which disputes the claims of this paper |fjq| . In Sec. 3.2, we have added some 
clarifying comments which we hope will ameliorate this dispute. In this note, we would like to 
discuss this in more detail. 

We are first of all pleased to notice that the authors of Ref. [^2| have to some extent agreed 
with our criticism to their papers |fl^ |, and modified their arguments accordingly. For instance, 
they now agree that the dominant contribution to the scattering within the "grey area" comes 
from local sources (although their exact definition of the "grey area" is somewhat different 
from ours). Also, they admit that the long-range, non-unitary, contribution to the scattering 



amplitude that they previously computed in Ref. [14| applies only for a coloured probe, and not 
for a dipole. 

But in spite of their agreement on these points, the authors of Ref. |4^] dispute our conclusion 
on "Froissart bound from gluon saturation", and instead conclude that the Froissart bound is 
violated by the perturbative evolution equations in Refs. J/], || |9|, pX| ]. In its essence (but using 
the language of our Sect. 3.2), the argument of Ref. [^2| can be formulated as follows: Assume 
that one starts the perturbative evolution at some early rapidity to, at which the point b± of 
interest lies far outside the initial "grey area" (b 3> Rh(t~o)), i.e., in the "white area", according 
to the terminology of Ref. p^] . In this perturbative setting, the dominant contribution to 
the scattering amplitude at b± is the long-range contribution, which decreases only as 



cf. eq. (3.34). With increasing r, the central region expands, and the point b± is eventually 
incorporated in the "grey area" at some "time" r%. For r > t%, the evolution proceeds locally 
at b±, but with an initial condition at r = t\ which decreases as a power of b±, rather than as 



an exponential. Thus, the first term in the r.h.s. of eq. ( 3,34| ) is replaced by e UJOlaT and the 



"blackness" condition M T (r±,b±) ~ 1 yields a black disk increasing exponentially with r. 

Although mathematically correct (at least, in its above formulation; we disagree with some 



of the intermediate steps and detail points in Ref. [42]), this argument is physically irrelevant: 
It simply signals a pathology of perturbation theory when this is used in a region where it is 
not supposed to work, namely, in the "white area", where physics is truly controlled by the 
confinement. In the discussion in Sect. 3.2, we pointed out that, under plausible assumptions 
about the non-perturbative physics in this region, the true evolution respects the exponential 
fall-off of the initial condition. 

In any case, our point in this paper was not to argue that, by integrating the perturbative 
evolution up to arbitrarily large r, one would generate cross-sections which saturate the Froissart 
bound indefinitely. Rather, we have shown that by applying perturbation theory where it is 
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expected to work — namely, in the central region at b < Rh(t), where the condition Q s (t, b) 3> 
Aqcd is satisfied by definition (cf. eq. (|2,12| )) — , one finds that (i) it is possible to compute the 
expansion rate of the black disk, and (ii) this rate is such that the Froissart bound is saturated. 
None of these points is trivial: 

(i) It was not a priori obvious that the expansion of the black disk is controlled by pertur- 
bation theory. We have demonstrated this for the case of a small dipole (Q 2 S> Aqcd)> f° r 
which we have shown that the "grey area" (= the region in which Q 2 (t, b) decreases from Q 2 
down to Aqcd) nas a considerable extent (cf. eq. ( |4.11[ )). This is the region which controls the 
expansion, and perturbation theory is applicable here. 

(ii) It was not a priori clear that the perturbative evolution of the black disk will respect the 
Froissart bound. Indeed, the saturated gluons from the black disk generate long-range forces 
which remain perturbative at impact parameters within the "grey area", and thus cannot be 
discarded by invoking confinement. (This was precisely the objection to Froissart bound raised 
in Ref. JT^j.) We have shown that these long-range forces, although present, are nevertheless 
harmless: Because the saturated gluons are globally colour-neutral (an aspect which has been 
overlooked in Ref. (]j~4]| ), they generate only dipolar forces, whose contribution to the scattering 
amplitude in the grey area remains smaller than the corresponding contribution of the local 
sources (which saturates the Froissart bound). 

The only place where non-perturbative physics has entered our argument was in providing 
the b± -dependence of the initial condition for the short-range contribution. On physical grounds, 
this is fixed by confinement, and thus is certainly of the exponential type (and not of the power- 
law type, as a naive perturbative evolution starting in the "white area" would predict). 

To conclude, Kovner and Wiedemann seem to be addressing a different question: Whether or 
not one can naively apply the perturbative evolution equation in all regions of impact parameter 
space — including over distances many times the pion Compton wavelength — , and still get 
sensible results. They negatively answer this question, and we definitely agree with them on 
this point. On the other hand, they further make the statement that one cannot compute the 
high-energy cross section using the ideas associated with (perturbative) saturation. The body 
of this paper clearly contradicts this claim. 
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Abstract 

We demonstrate that the dipole-hadron cross-section computed from the non-linear evo- 
lution equation for the Colour Glass Condensate saturates the Froissart bound in the case 
of a fixed coupling and for a small dipole (Q 2 ^> ^qcd)- That is, the cross-section increases 
as the logarithm squared of the energy, with a proportionality coefficient involving the pion 
mass and the BFKL intercept (a s -/V c /7r)4 In 2. The pion mass enters via the non-perturbative 
initial conditions at low energy. The BFKL equation emerges as a limit of the non-linear 
evolution equation valid in the tail of the hadron wavefunction. We provide a physical pic- 
ture for the transverse expansion of the hadron with increasing energy, and emphasize the 
importance of the colour correlations among the saturated gluons in suppressing non-unitary 
contributions due to long-range Coulomb tails. We present the first calculation of the satu- 
ration scale including the impact parameter dependence. We show that the cross-section at 
high energy exhibits geometric scaling with a different scaling variable as compared to the 
intermediate energy regime. 
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1 Introduction 



One of the striking features of the physics of strong interaction is that at high energies, cross 
sections are slowly, but monotonously, increasing with s (= the total center-of-mass energy 
squared). For instance, the data for the pp and pp total cross sections at high energy can be 
reasonably well fitted by both s s (with 6 ~ 0.09) and In 2 s, although the In 2 s growth appears 
to be favored by recent investigations Q. 

At a theoretical level, it was proven many years ago that in the limit s — > oo, the hadronic 
total cross sections must rise no faster than (ir/m 2 ) In 2 s |2[ ||, |[]. This is the Froissart bound, 
which is a consequence of very general principles, such as unitarity and analiticity, but does 
not rely on any detailed dynamical information. So, in reality, this bound may very well be 
not saturated, although the measured cross-sections seem to do so (at least, in so far as they 
are consistent with a squared log dependence upon the energy). In fact, a simple, qualitative 
mechanism realising this In 2 s growth has been proposed by Heisenberg ||] already before Frois- 
sart bound has been rigorously proven (a similar argument is briefly discussed by Froissart [0]). 
However, we are not aware of any field-theoretical implementation of this, or other, argument. 
(See also Ref. || for a recent review and more references, and Ref. for a recent analysis of 
the data.) 

The goal of this paper is to see how the Froissart bound can be consistent with modern 
pictures of high energy strong interactions. To keep the discussion as simple as possible while 
still encompassing the interesting physics, we shall consider the scattering of a colour dipole off a 
hadronic target at very high energy. The "colour dipole" may be thought of as a quark-antiquark 
pair in a colourless state, like a quarkonium, or a fluctuation of the virtual photon in deep 
inelastic scattering. More generally, arbitrary hadronic probes can be considered as collections 
of "colour dipoles" at least in the approximation in which the number of colours N c is large (so 
that a gluon excitation can be effectively replaced by a qq pair). For our approximations to be 
justified, we shall assume that the dipole is "small", in the sense that it has a small transverse 
size, rj_ <C 1/Aqcd, or a large transverse resolution: Q 2 = 1/r 2 3> A-qcd- 

On the particular example of the dipole-hadron scattering, we shall develop an argument 
involving both perturbative and non-perturbative features of QCD, like non-linear quantum 
evolution, parton saturation, and confinement, which will lead us to the conclusion that the 
total cross-section respects, and even saturates, the Froissart bound (at least, in the case of a 
fixed coupling, that we shall exclusively consider in this paper). In its essence, our argument 
may be viewed as a modern version of Heisenberg's original mechanism. But our main point is 
to place this mechanism within the context of our present theoretical understanding, and clarify 
to which extent perturbation theory may play a role in determining this result. In brief, we shall 
demonstrate that this mechanism is naturally realized by the solutions to non-linear evolution 



equations derived within perturbation theory, in Refs. [7], ||, ||, 10(1 , but with non-perturbative 



initial conditions at low energy. In particular, our analysis will confirm, clarify and extend the 
conclusions reached in an early study [11] based on the GLR equation [12, 13], which is a simplified 
version of the non-linear evolution equations that we shall use in this paper. 

The fact that perturbation theory is the appropriate tool to describe the high energy be- 
haviour of hadronic cross-sections is by itself non-trivial, and deserves some comment. The first, 
and most certainly true, objection to the use of weak coupling methods comes from the fact 
that the Froissart bound involves the pion mass, and this surely must arise from confinement. 
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And, indeed, this is how the pion mass enters also our calculations: via the non-perturbative 
initial condition that we shall assume, and which specifies the impact parameter dependence of 
the dipole-hadron scattering amplitude at low energy. But starting with this initial condition, 
we shall then study its evolution with increasing energy within perturbation theory, and show 
that the resulting cross-section saturates the Froissart bound at high energy. 

The adequacy of perturbation theory even for such a limited purpose — the study of the 
quantum evolution of the cross-section with s — is still non-trivial ||, and has been in fact 
disputed in recent papers pM . It is first of all clear that ordinary perturbation theory, in the form 
of the BFKL equation |l5| — which resums the dominant radiative corrections at high energy, 
but neglects the non-linear effects associated with high parton densities — , fails to describe the 
asymptotic behaviour at high energy: The BFKL equation predicts a power-law growth for the 
total cross section: a ~ s WQs , with u> = 4 In 2 and a s = a s N c /ir, which clearly violates the 
Froissart bound. Besides, with increasing energy, the solution to the BFKL equation "diffuses" 
towards smaller and smaller transverse momenta, thus making the applicability of perturbation 
theory questionable. 

But the objections in Ref. [14] actually apply to the more recently derived non-linear evo- 
lution equations |@, |8| [9|, [ji]], in which non-linear effects cure the obvious pathologies of the 
BFKL equation. (Alternative derivations of some of these equations, at least in specific limits, 
have been given in Refs. ]i~2[ , |i~3| , 16, [17]]. See also Refs. [18, 19] for recent reviews and more 
references.) Specifically, the non-linear effects associated with the high density of gluons in the 
hadron light-cone wavefunction (the "Colour Glass Condensate" [|l(J) lead to gluon saturation 
[12, 2C, 21, 22], with important consequences for the high-energy scattering: First, this intro- 
duces a hard intrinsic momentum scale, the "saturation scale" , which is a measure of the density 
of the saturated gluons in the impact parameter space, and grows like a power of the energy. 
This limits the infrared diffusion [^] and thus provides a better justification for using weak 
coupling methods at high energy. Second, this ensures the unitarization of the dipole-hadron 



scattering amplitude at fixed impact parameter [gj, [25|, gg, [27], That is, with increasing 
energy, the hadron eventually turns "black" (i.e., the dipole is completely absorbed, or the scat- 
tering amplitude reaches the unitarity limit) at any given point b± in the impact parameter 
space. 

However, by itself, the unitarization at given b± is not enough to guarantee the Froissart 
bound for the total cross-section, which involves an integration over all the impact parameters. 
Indeed, as a quantum mechanical object, the hadron has not a sharp edge, but rather a dif- 
fuse tail, so, with increasing energy, the "black disk" can extend to larger and larger impact 
parameters. The radial expansion of the black disk is controlled by the scattering in the sur- 
rounding "grey area", where the gluon density is relatively low and the BFKL equation still 
applies. Given the problems of the BFKL equation alluded to before, it is not a priori clear 
whether this expansion is slow enough to ensure that the Froissart bound is respected. This 
would require a black disk radius which increases at most logarithmically with s. In Heisenberg's 
original argument (but using the current terminology), such a logarithmic increase was ensured 
by a compensation between the power-law increase of the scattering amplitude with s and its 
exponential decrease with b = \b±\. Such an exponential fall-off at large impact parameters is, of 
course, a true property of full QCD, and also the crucial assumption about our initial condition, 
but it is not clear whether this property can be preserved by the perturbative quantum evolution, 
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which involves massless gluons and therefore long range interactions. 

In fact, it was the main point of Ref. to argue that the Coulomb tails associated with the 
saturated gluons should replace the exponential fall-off with b = \b±\ of the initial distribution 
by just a power-law fall-off, which would be then too slow to ensure the Froissart bound: the 
corresponding black disk radius would increase as a power of s. 

However, as we shall explain in this paper, the argument in Ref. [0] is irrelevant for the problem 
at hand, and also incorrect in its original formulation]]]. Part of the confusion in Ref. |Q comes 
from non-recognizing that the saturated gluons are actually colour neutral over a relatively short 



distance (of the order of the inverse saturation scale) [22, jig] , and thus cannot produce Coulomb 
tails at large distances. Rather, the colour field produced by the saturated gluons is merely a dipolar 
field, whose fall-off with b is sufficiently fast to respect the Froissart bound for the scattering of an 



external dipole. (In Ref. |14], this was masked by the fact that the authors were truly computing 



the scattering of a coloured external probe, although their verbal arguments were formally developed 
for a "dipole".) As we shall see, for the relevant impact parameters, this long-range dipole-dipole 
scattering is subleading at high energy as compared to the short-range scattering off the local sources. 

More precisely, we shall demonstrate that the region which controls the evolution of the cross 
section is the "grey area" outside, but close to the black disk. In this area, we expect perturbation 
theory to apply, since the local saturation scale is much larger than Aqcd- For an incoming dipole 
in this area, the dominant interactions are those with the non-saturated colour sources within 
a "saturation disk" (i.e., a disk with radius equal to the inverse saturation scale) around the 
impact parameter b±. This is a consequence of two physical facts: a) the non- linear effects limit 
the contribution of the distant colour sources, and b) being colourless, the dipole couples only 
to the local electric field (as opposed to the long-range gauge potentials), so it is less influenced 
by colour sources which are far away. 

Since, moreover, the local saturation length is much shorter than the typical scale ~ 1/Aqcd 
for transverse inhomogeneity in the hadron (this is where the condition that the dipole is "small" , 
i.e., Q 2 S> Aq CD , is essential), it follows that the quantum evolution proceeds quasi-locally in 
the impact parameter space. This in turn implies that, within the grey area, the ^"dependence of 
the scattering amplitude factorizes out, and is therefore determined by the initial condition at low 
energy. On general physical grounds, we shall assume this initial condition to have an exponential 
fall-off as e - 2m ^ b at large distances (indeed, pion pairs must control the long distance tail of the 
hadron wavefunction; see, e.g., pM and Refs. therein). Then an argument similar to the original 



one by Heisenberg can be used to conclude that the total cross-section increases like In 2 s. 

The coefficient in front of In 2 s in our final result is also interesting, as it reflects the subtle 
interplay between the perturbative and non-perturbative physics contributing to this result. 
Specifically, we shall find that, for any hadronic target, 

vr {(va s ) 2 2 n n 

a w ^ — in s as s — > oo, f 1.1 ) 

2 m~ 

where the pion mass in the denominator enters via the exponential fall-off of initial condition, 
while the factor ua s = 4(ln2)a s N c /ir in the numerator is recognized as the "BFKL intercept". 
This latter comes up because, in deriving this result, we will have to consider the solution to 



1 See however the new preprint J42| where a modified version of this argument has been presented; we shall 
comment on this new argument in the Note added at the end of this paper. 
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the BFKL equation at large energy for fixed Q 2 . This is of course the limit for which the 
BFKL equation has been originally proposed |]l5[ , but not also the limit used in more recent 
applications of this equation within the context of saturation (e.g., in studies of the saturation 



scale [ 2X|, g^, |3(J, or of the "geometric scaling" J29fl) for a homogeneous hadron. 

The difference with Refs. pl| , |29| , |30[ | occurs because we consider here a different physical 
problem: Rather than studying the quantum evolution at a fixed impact parameter — which 
would then limit the applicability of the BFKL equation to not so high energies, such that the 
local saturation scale remains below Q 2 — , we rather follow the expansion of the black disk 
with increasing s, and use the BFKL equation only in the outer grey area at sufficiently large 
b, where the gluon density remains small even when the energy is large. In other terms, by 
increasing the energy at fixed Q 2 and simultaneously moving towards larger impact parameters, 
one always finds a corona where the local saturation scale is much smaller than Q, but much 
larger than Kqcd- For points b± in this corona, we can prove the factorization of the scattering 
amplitude into a b± -dependent "profile function" which is determined by the initial condition 
at low energy, and an energy- and Q 2 -dependent factor which can be computed by solving the 
homogeneous (i.e., no b± -dependence) BFKL equation. In the high energy limit at fixed Q 2 , 
this calculation yields the cross-section in eq. (|l . l| ) . 

In addition to eq. (LI), we shall find a variety of new features within our approach. For 
example, we shall compute the impact parameter dependence of the saturation scale, and dis- 
cover an entirely nontrivial structure for the radial distribution of matter inside the hadron, and 
its evolution with increasing energy. Related to that, we shall find that the geometric scaling 
arguments which have been used to characterize deep inelastic scattering at high energies be- 
come modified by the appearence of two different scales (associated both with gluon saturation) 
which govern the geometric scaling of the total cross-section in different ranges of the energy. 
As a result of our analysis, we shall derive an intuitive picture for the expansion of the hadron 
in the transverse plane. 

The outline of the paper is as follows: 

In the second section, we qualitatively and semi-quantitatively discuss how the Froissart 
bound becomes saturated within the framework of our knowledge of gluon saturation, non- 
linear evolution, and its linearized, BFKL, approximation. This discussion introduces the main 
arguments to be demonstrated by the technical developments in the nextcoming sections. 

In the third section, we study the properties of a non-linear evolution equation for the 
scattering amplitude, the Balitsky-Kovchegov (BK) equation J?], We argue that the dominant 
contribution to the scattering in the grey area comes from virtual dipoles whose size is much 
smaller then the saturation length. This leads us to conclude that the corresponding scattering 
amplitude factorizes in the way alluded to before. In the rest of the third section, we explore 
the physical significance of this factorization by using the efective theory for the Color Glass 
Condensate [|n], |H| . We show that an essential ingredient for factorization and Froissart bound 
is the colour-neutrality of the saturated gluons within the black disk. 

In the fourth section, we compute the radius of the black disk using the solution to the 
BFKL equation. Then, we extend our results to compute the impact parameter dependence of 
the saturation momentum. We identify two range of values of the impact parameter at which 
also the saturation scale factorizes, i.e., it is the product of an exponentially decreasing function 
of b± times a factor increasing like a power of s. These two ranges correspond to the two limits 
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of the BFKL solution alluded to before: For b± sufficiently close to the center of the hadron, 
the increase with the energy is the same as for the saturation scale of a homogeneous hadron, 
previously studied in Refs. pi], |2^, |30f1 . This comes from a solution to the BFKL equation in 
the intermediate regime where InQ 2 ~ a s Ins. On the other hand, near the edge of the hadron, 
the increase with s is rather controlled by the high energy solution at InQ 2 <C a s lns. 

In the fifth section, we discuss the implication of the different behaviours found at small and 
large impact parameters for geometrical scaling in deep inelastic collisions. 

In the last section, we summarize our results and present our conclusions. 



2 Saturation scale and the Froissart bound 

The dipole-hadron collision will be considered in a special frame, the "dipole frame" [24, [lS|| , 
in which the physical interpretation of our results becomes most transparent: This is the frame 
in which the effects of the quantum evolution are put solely in the wavefunction of the hadron, 
which carries most of the total energy. This being said, it should be stressed that our final 
results are independent of this choice of the frame — although their interpretation may look 
different in other frames — since at a mathematical level they are based on boost invariant 
equations, namely, the non- linear evolution equation for the scattering amplitude R |8L [K| and 
its linearized, BFKL [GiJ, approximation. In fact, the relevant non-linear equation — namely, 
the Balitsky-Kovchegov (BK) equation to be presented in Sect. 3 — has been independently 
derived in the hadron rest frame [0, ||, where the evolution refers to the incoming dipole wave- 
function, and in the infinite momentum frame (which for the present purposes is equivalent to 
the "dipole frame" alluded to before) , from the evolution of the "Colour Glass Condensate" [[h| . 

So, let us consider an incoming dipole of transverse size r±, with Q 2 = l/r 2 S> Aq CD , which 
scatters off the hadron target at large invariant energy squared s, or rapidity gap r = ln(s/Q 2 ). 
In the dipole frame, most of the total energy is carried by the hadron, which moves nearly at 
the speed of light in the positive z direction. Moreover, any further increase in the total energy 
is achieved by boosting the hadron alone. Thus, the dipole rapidity Tdi po ie is constant, and 
chosen such as a s Tdi po ie "C 1, so that we can neglect higher Fock space components in dipole 
wavefunction: the dipole is just a quark-antiquark pair, without additional gluons. On the other 
hand, the hadron wavefunction has a large density of small-x gluons, which increases rapidly 
with r. Here, x is the longitudinal momentum fraction of the gluons which participate in the 
collision, and is related to the rapidity r as r = ln(l/x). 

In this special frame, the unitarization effects in the dipole-hadron collision can be assimilated 
to the saturation effects in the target wavefunction. This is so since it is the same momentum 
scale, namely the saturation scale Q 2 , which sets the border for both types of effects. 

A priori, Q s is an intrinsic scale of the hadron, proportional to the gluon density in the 



transverse plane at saturation |12L 13 



Q 2 s (r,b±) ~ ^_xG(x,Q 2 ,6 ± ). (2.1) 
In this equation, 
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is the number of gluons with longitudinal momentum fraction x and transverse size Axj_ ~ 1/Q 
per unit rapidity and per unit transverse area, at transverse location b±. The more standard 
gluon distribution xG(x, Q 2 ) is obtained by integrating xG(x, Q 2 , b±) over all the points b± in 
the transverse plane. Note that the distribution in b± is indeed a meaningful quantity since 
we consider gluons with relatively large transverse momenta, Q 2 ^> Aq CD , which are therefore 
localized over distances Ax± much smaller than the typical scale for transverse inhomogeneity 
in the hadron, namely, Ab± ~ 1/Aqcd- 

The saturation scale separates between two physical regimes: At high transverse momenta 
k± S> Q s (t, b±), we are in the standard, perturbative regime: The gluon density is low, but it in- 
creases very fast, (quasi)exponentially with r, according to linear evolution equations like BFKL 
[15] or DGLAP [pi] ]. At low momenta k± ^ Qs{t, b±), the non-linear effects are strong even if 
the coupling is weak, and lead to saturation: The gluon phase-space density is parametrically 
large, dN/drd 2 k±d 2 b± ~ l/a s , but increases only linearly with r [^, ||^, |l8| . 

Prom eq. ( |2.lD , one expects the saturation scale to increase rapidly with r, so like the gluon 
distribution at high momenta. For a hadron which is homogeneous in the transverse plane (no 



dependence upon b±), the r-dependence of Q s is by now well understood |21j, |29|, |30|, |32j, |33|, (23(1 , 
and will be reviewed in Sect. 4 below. Understanding the r- and ^-dependences of the saturation 
scale in the general inhomogeneous case is intimately related to the problem of the high energy 
behaviour of the total cross-section, so this will be a main focus for us in this paper. 

To appreciate the relevance of the saturation scale for the dipole scattering, note that a small 
dipole, i.e., a dipole with transverse size r± <C 1/Q S (t, b± ) , where b± is the impact parameter, 
couples to the electric field created by the colour sources in the target. Thus, its scattering 
amplitude is proportional to the correlator of two electric fields, which is the same as the gluon 
distribution xG(x, Q 2 ,b±) with Q 2 ~ l/r\ @ : 

TT 2 Ct C 

N T (r±,b±) ~ r\ N2 S _ 1 xG(x, l/rj_,b±) (single scattering). (2.3) 

This equation, together with the solution xG(x, l/rj_, b±) to the BFKL equation, predicts an 
exponential increase of the scattering amplitude with r. If extrapolated at high energy, this 
behaviour would violate the unitarity requirement M T < 1. However, eq. (2.3) assumes single 
scattering, so it is valid only as long as the gluon density is low enough for the condition 
A/" r (rj_, 6j_) <C 1 to be satisfied. At high energies, where the gluon density is large, multiple 
scattering becomes important, and leads to unitarization. Assuming that the successive collisions 
are independent, one obtains [24] : 

A/" r (rj_, b±) ~ 1 — exp j — r l ^^ S xG(x, l/rj_, 6j_)| (multiple scattering), (2.4) 

which clearly respects unitarity. 

So far, this only demonstrates the role of multiple scattering in restoring unitarity. The deep 
connection to saturation follows after observing that the condition for multiple scattering to be 
important — that is, that the exponent in eq. (|2.4|) is of order one — is the same as the condition 



(2.1) for gluon saturation in the hadron wavefunction. This is natural since the dipole is a direct 
probe of the gluon distribution in the hadron, so the non-linear effects in the dipole-hadron 
scattering and in the gluon distribution become important at the same scale. But this also 
shows that, in the non-linear regime at Q 2 ^ Qs(t, b±) one cannot assume independent multiple 
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scatterings, as in eq. fl2.4|): Rather, the dipole scatters coherently off the saturated giuons, with 



a scattering amplitude which, unlike eq. (^4), cannot be related to the gluon distribution (a 
2-point function) alone, but involves also higher n-point functions. This amplitude satisfies a 
non- linear evolution equation to be discussed in Sect. 3. 

But it is nevertheless true that, as suggested by eq. (|2.4[), the scattering amplitude becomes 
of order one in the saturation regime at Q 2 5; Qlij, b±). (This has been verified via both analytic 



p2| , 22, |29| and numerical |32|, 33, |23| investigations of the non-linear evolution equation.) Thus, 
when increasing the energy at fixed Q 2 = 1/r 2 ,, the unitarity limit M T (r±, b±) = 1 is eventually 
reached at any given impact parameter b± : for the incoming dipole, the hadron looks locally 
"black". 

However, the unitarization of the local scattering amplitude is not enough to guarantee the 
Froissart bound for the total cross-section. The latter is obtained by integrating the scattering 
amplitude over all impact parameters: 

a(r,r ± ) = 2 / d 2 b ± M T (r ± ,b ± ). (2.5) 



It is easy to see that difficulties with the Froissart bound can arise only because the hadron does 
not have a sharp edge. Indeed, if in transverse projection the hadron was a disk of finite radius 
i?0) then for sufficiently large energy it would become black at all the points within that disk, 
and the total cross-section would saturate at the geometrical value 27tRq. 

But in reality a hadron is a quantum bound state of the strong interactions, so its wave- 
function has necessarily an exponential tail, with the scale set by the lowest mass gap in QCD, 
that is, the pion mass. Specifically, in the rest frame of the hadron, the distribution of matter 
is typically of the Woods-Saxon type: 

P(r) = ( 2 - 6 ) 



1 + exp (^a) 



where Rq is the typical radial size of the hadron under consideration (this increases as A 1 / 3 for 
a nucleus with atomic number A), while the thickness a = 1/2771,^ is universal (i.e., the same 
for all hadrons). This latter involves twice the pion mass because of isospin conservation: At 
high energy, one probes the giuons in the hadron wavefunction, and giuons have zero isospin, so 
they couple to the external probe via the exchange of (at least) two pions. It is therefore 2m,,- 
which controls the exponential fall off of the scattering amplitude, or of p(r), at large distances: 
p(r) oc e -i™-(r-Ro) f or b-R > l/2m w . 

Since, moreover, small-x giuons have large longitudinal wavelengths, a high energy scattering 
is sensitive only to the distribution integrated over z, that is, to the transverse profile function : 

2 7 ,2 



fdzp^b 2 ± + z 



Jdzp(z) 

(b = \b±\) which is normalized at the center of the hadron: S(b = 0) = 1. Note that, independent 
of the detailed form of p(r) in the central domain at r < Rq, the function S(b) decreases 
exponentially0, 5(6) ~ e -^(b-R )^ for 5 — Rq > l/2m n . 

2 But power law corrections to this exponential decrease may be numerically important |28|. 
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Based on these considerations, we shall assume that the scattering amplitude for the (rel- 
atively low energy) dipole-hadron scattering in the target rest frame — which is our initial 
condition for the quantum evolution with r — has the following factorized structure: 

K (r±,b ± ) = M T0 (r ± )S(b), (2.8) 

with To = Tdipoie and the profile function S(b) introduced above. At low energy and high Q 2 , the 
factorization of the ^"dependence is natural, since consistent with the DGLAP equation (see, 



e.g., (T^]). For what follows, the crucial feature of the initial amplitude (2.8) is its exponential 
fall off ~ e -^(b-Ro) at i arge distances b > R . 

Starting with this initial condition, we increase the energy by boosting the hadron to higher 
and higher rapidities. Clearly, the gluon distribution xG(x, l/rj_, b±) at any b± will increase with 
r, and the BFKL approximation suggests that this increase should be exponential. Thus, even 
points b± which were originally far away in the tail of the hadron wavefunction (J)—Rq ^> l/2m n ), 
and did not contribute to scattering at the initial rapidity To, will eventually give a significant 
contribution, and even become black when the energy is high enough. That is, with increasing 
r, the black disk may extend to arbitrarily large impact parameters. 

At this point, it is useful to introduce some more terminology. By the "black disk" we mean 
the locus of the points b± in the transverse plane at which the unitarity limit A^-(rj_,fej_) = 1 
has been reached in a dipole-hadron collision at rapidity r and transverse resolution Q 2 = 
Equivalently, the condition Q 2 < Q 2 (t, b±) is satisfied — i.e., the gluons with transverse 
momenta ~ Q are saturated — at all the points in the black disk. Given the shape of the initial 
matter distribution Q2.6| ) — which is isotropic and decreases from the center of the hadron 
towards its edge — it is clear that the "black disk" is truly a disk, with center at b = 
and a radius R(t,Q 2 ) which increases with r and decreases with Q 2 . The black disk radius 
is determined by any of the two following conditions (for more clarity, we shall often rewrite 
Af T (Q 2 ,bj_) =fif T (r± = 1/Q,b±) in what follows): 

Af T (Q 2 , b±) = k for b = R(t,Q 2 ), (2.9) 

or 

Q 2 s (r,b±) = Q 2 for b = R(t,Q 2 ), (2.10) 

which are equivalent since, in turn, the saturation scale is defined by: 

M T (Q 2 ,b ± ) = K for Q 2 = Q 2 (r,b ± ). (2.11) 

In these equations, k is a number smaller than one, but not much smaller (e.g., k = 1/2), whose 
precise value is a matter of conventionP]. For qualitative arguments, and also for quantitative 
estimates at the level of the approximations to be developed below, one can take k = 1. 

We shall also need below the "edge of the hadron" at rapidity r, by which we mean the radial 
distance Rh{j) at which the saturation scale becomes of order Aqcd (this would correspond to 
the black disk seen by a large dipole with resolution Q 2 ~ Aq CD , e.g., a pion) : 

Q 2 s (r,b ± ) = A 2 QCD for b = R H (r). (2.12) 



3 The separation between the saturation regime at small Q 2 and the low-density regime at high Q 2 being not 
a sharp one, ther e is some ambiguity in defining the borderline Q 2 = Q 2 {t, b±). This is fixed by choosing the 



number k in eq. ( 2.11 ) 



9 



With these definitions at hand, we now return to the discussion of the Froissart bound. At 
sufficiently high energy, the total cross-section is dominated by the contribution of the black 
disk (this will be verified in Sect. 5) : a ~ obd with 

°bd(t,Q 2 ) = 2 J d 2 bj_ Q(R(t, Q 2 ) - b)M T {Q 2 , b ± ) « 2ttR 2 (t,Q 2 ). (2.13) 

Thus, the question about the Froissart bound becomes a question about the expansion of the 
black disk with r : To respect this bound, R(t, Q 2 ) must grow at most linearly with r. 

One can easily construct a "naive" argument giving such a linear increase (this is similar in 
spirit to the old argument by Heisenberg g): Starting with an initial distribution like ( |2.6| ), 
assume that, with increasing r, the gluon density increases in the same way at all the points b± 
(outside the black disk), so that the b± -dependence of the scattering amplitude factorizes out, 
and is fixed by the initial condition: 

K(Q 2 ,b ± ) « S(b ± )K(Q 2 ). (2.14) 

Let us furthermore assume that the function M T (Q 2 ) at large r is given by standard perturbation 
theory, that is, by the solution to the BFKL equation at high energy : N T (Q 2 ) oc e ujasT , 
where ui = 4 In 2 and a s = N c a s /n. Under such (admittedly crude) assumptions, the scattering 
amplitude at large r and b 3> Rq is given by: 

K(Q 2 ,b ± ) « \ ^ e^ T e- 2m - b , (2.15) 

V Q 

where we have also included the leading Q 2 -dependence of the asymptotic BFKL solution at 
high energy []l5[ . (A 2 is some arbitrary reference scale, of order Aq CD .) This expression together 



with the saturation condition (2.9) imply: 



and the resulting cross-section saturates the Froissart bound indeed: 
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Our main objective in this paper will be to show that the above, seemingly naive, argument 
is essentially correct, and the results in eqs. ( 2.16| ) and ( 2.1 7| ) are truly the predictions of the 



non- linear evolution equation for N T (Q 2 , b±) at sufficiently large r. This is non-trivial since the 
naive argument might go wrong for, at least, two reasons: 

i) At very high energies, the non-linear effects become important, and the use of the BFKL 
equation becomes questionable. For instance, the unitarization of the local scattering amplitude 
J\f T (Q 2 ,6j_) is precisely the result of such non-linear effects, which are taken into account by 
replacing the BFKL equation with the BK equation. 

ii) Although non-linear, the quantum evolution described by the BK equation remains per- 
turbative, so it involves massless gluons and long-range effects which could not only invalidate the 
factorization property ( p,14[) , but also replace the exponential fall-off of the initial distribution 
by just a power-law fall-off (an eventuality in which the Froissart bound would be violated). 
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Nevertheless, as we explain now (and will demonstrate in Sect. 3 below), none of these two 
objections apply to the problem of interest. Indeed: 

i) However large is r, there exists an outer corona at R(t, Q 2 ) < b < Rh{t) where the hadron 
looks still "grey", i.e., where M T (Q 2 , b±) -C 1 and the BFKL approximation applies. It is this 
"grey area" which controls the expansion of the black disk, and therefore the evolution of the 
total cross-section at high energy. 

ii) The quantum evolution within the grey area is quasi- local in b± , because of the non- linear 
effects which limit the range of the relevant interactions to Ab± <C 1/Q S (t, b± ) . 

To be more specific, note that, in order to study the expansion of the black disk with r, 
one needs to consider the evolution of the scattering amplitude N T (Q 2 ,b±) at points b± which 
lie outside the black disk, but relatively close to it. Indeed, when r — > r + dr with a s dr ~ 1 
(which is the typical increment in the high energy regime of interest: a s r 3> 1), the black disk 
expands by incorporating the points b± within the range R<b < R + dR with R = R(t, Q 2 ) 



and dR ~ cf. eq. (2.16). Such points are sufficiently far away from the black disk for 

the local saturation scale to be small compared to Q 2 — that is, they are in the "grey area" 
— , but also sufficiently far away from the edge of the hadron (cf. eq. fl2.12p ) for Q s (r,b±) 
to be a "hard" scale. That is, the following conditions are satisfied for any b± of interest: 
^■QCD ^ Ql{ T ^) *^ Q 2 - Both inequalities are important for our argument, as we explain now: 
The fact that Q 2 s (r,b) < Q 2 ensures that the dominant contribution to the evolution 
dAf r (Q 2 , b±) of the scattering amplitude comes from nearby colour sources, i.e., from the sources 
which are located within a saturation disk around b± (see also Fig. ||): 

|2l - M<< d^r (2 ' 18) 

and therefore lie themselves inside the grey area. This is so because sources which lie further 
away are shielded by the non-linear effects. Besides, being a colour singlet, the dipole is not 
sensitive to the long-range gauge potentials. 

The fact that 1/Q s (t, b±) <C 1/Aqcd implies that the transverse inhomogeneity in the 
hadron can be neglected when computing the contribution of such nearby sources to dN T (Q 2 , b±). 
That is, all the relevant sources act as being effectively at the same impact parameter, equal to 
b±. This explains the factorized expression ( 2.14| ) for the scattering amplitude. 



Since, moreover, M T {Q 2 , z±) <C 1 for any z± satisfying ( 2.1g| ), it follows that the function 
M T (Q 2 ) in eq. ( 2.14 ) can be computed by solving the linearized (and homogeneous) version of 



the BK equation, namely the BFKL equation without b± dependence. 

Clearly, it was essential for the previous arguments that the dipole is "perturbative" : Q 2 
Aq CD . This ensures that the separation Rh(t) — R(t, Q 2 ) between the black disk and the hadron 
edge (i.e., the width of the "grey area") is sufficiently large for the condition Q 2 (r,b) 3> Aq CD 
to apply at all the points b± of interest. Besides, one can argue that Q 2 is the scale at which 
the QCD coupling should be evaluated (see the discussion in the Conclusions). 

A factorization assumption similar to eq. ( |2.14| ) has been already used in the literature, in 
particular, in relation with the Froissart bound []ll|, |35|| , and also as an Ansatz in the search for 
approximate analytical [|32| or numerical |^] solutions to the BK equation. But in previous work, 
this assumption was always based on experience with the (homogeneous) DGLAP equation, and 
not duly justified in the small-x regime. 
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Figure 1: A pictorial representation of the dipole-hadron scattering in transverse projection 
(only half of the hadron disk is shown). The 6-dependence of the saturation scale illustrated by 
the lower plot is the one to be found in Sect. 4.3. 
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That such a factorization is highly non-trivial in the presence of long-range gauge interac- 
tions is also emphasized by a recent controversy about this point, put forward in Ref. [ [L"f| . 
Specifically, in Ref. [jl4| it has been shown that, as far as the scattering of a coloured probe 
off the hadron is concerned, the long-range fields created by the saturated gluons provide a 
non-unitarizing contribution to the respective cross-section. On the basis of this example, the 



authors of Ref. [14] have concluded that the non-linear BK equation provides "saturation with- 
out unitarization" . In Sect. 3.2 below, we shall carefully and critically examine the arguments in 
Ref. |jl4|], and demonstrate that, for the physically interesting case where the external probe is a 
(colourless) dipole, there is no problem with unitarity at all. The long-range interactions between 
the incoming dipole and the saturated gluons give only a small contribution to the scattering 
amplitude in the grey area, because the saturated gluons form themselves a dipole (i.e., they are 
globally colour neutral), and the dipole-dipole interaction falls off sufficiently fast with the sep- 
aration between the two dipoles. The dominant contribution comes rather from the short-range 
scattering within the grey area (cf. eq. ( 2.18[ )), for which the factorization assumption ( |2.14| ) is 
indeed justified. 

3 Quantum evolution and black disk radius 

In the effective theory for the Colour Glass [pj|], the dipole-hadron scattering is described as 
scattering of the qq pair off a stochastic classical colour field which represents the small-x gluons 
in the hadron wavefunction. At high energy, one can use the eikonal approximation to obtain: 

K(r±,b ± ) = l-S T (x ± ,y ± ), S T (x ± , y ± ) = (ti(V^(x ± )V(y ± ))) T , (3.1) 

with rj_ = x±—y± the size of the dipole and b± = (x± + y±)/2 the impact parameter (the quark 
is at xj_, and the antiquark at y±). The S'-matrix element S T ( x± , y± ) involves the Wilson lines 
(path ordered exponentials along the straightline trajectories of the quark and the antiquark) 
and V built with the colour field of the target hadron. For instance, 

V\x ± ) = Pexp ji 5 jT dna a v (x ± )t a Y (3.2) 

where a v (x±) is the stochastic "Coulomb field" created by color sources (mostly gluons) at 
rapidities r' < r, and has longitudinal support at (space-time) rapidity^ n < r. Thus, the 
integral over n in eq. fl3.2| ) is in fact an integral over the longitudinal extent of the hadron (in 
units of space-time rapidity) seen by the external probe in a scattering at rapidity r. That is, 
the actual width of the hadron depends upon the energy of the collision. This is so since, with 
increasing energy, gluon modes with larger and larger longitudinal wavelengths participate in 



the collision, so that the hadron looks effectively thicker and thicker [ 10 1 



The brackets in the definition (3.1) of the S'-matrix element refer to the average over all the 



configurations of the classical field with some appropriate probability distribution W T [a] : 

(tr(V+(s_L)V(yj_))) T = / [da] tr(V^x ± )V(y ± ))W T [a\. (3.3) 



4 The space-time rapidity is defined as rj = ln(a; _ P + ), where x~ — (t — z)/y2 is the light-cone longitudinal 
coordinate, and P + is the light-cone momentum of the hadron. Light-cone vector notations are defined in the 
standard way, that is, v = (1/V2)(v° ± v 3 ). With the present conventions, the hadron is a right mover, while 
the dipole is a left mover. 
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This probability distribution is not known directly, but its variation corresponding to integrating 



out gluons in the rapidity window (r, r + dr) can be computed [34, 10]. This leads to a functional 
evolution equation for W T [a;] whose precise form is not needed here (see Ref. |H| for details). 
Suffices it to say that, via equations like (|3.3|) , the functional equation for W T [a] can be translated 
into an hierarchy of ordinary evolution equations for the n-point functions of the Wilson lines . 
This procedure yields the same equations as obtained by Balitsky within a different approach, 
which focuses directly on the evolution of Wilson line operators 0. (The fact that the infinite 
hierarchy of coupled equations by Balitsky can be reformulated as a single functional equation 
has been first recognized by Weigert ||.) In the limit where the number of colours N c is large, 
a closed equation can be written for the 2-point function fl3,l| ) (with a s = N c a s /ir) : 

d c , s _ f d 2 z± (x± - y±) 2 

-S T {x±,y±) = -a s 



dr J 2tt (x± - z±) 2 (y± - z±) 2 

x(S T (x±,y±) - S T (x±, z±)S T (z±,y±)y (3.4) 

The same equation has been derived independently by Kovchegov |8| within the Mueller's dipole 
model [25]. We shall refer to eq. ( |3.4| ) as the Balitsky-Kovchegov (BK) equation. 

3.1 Scattering in the grey area 

In this subsection, we shall study the scattering amplitude in the grey area, and prove the 
factorization property Q2.14| ). For more clarity, we shall formulate our arguments at the level of 
the BK equation. But one should keep in mind that our final conclusions are not specific to the 
large N c limit: the same results would have been obtained starting with the general non- linear 
evolution equations in Refs. |7|, PJH]]. 

Specifically, we shall use eq. (|3.4| ) to demonstrate that the dominant contribution to dS T /dr 
in the grey area comes from short-range scattering, i.e. from points z± such that 

1*1 -b±\ « 1 « -t— — • (3.5) 
Q s {T,b ± ) Aqcd 

As explained in Sect. 2, the impact parameters of interest are such that the following in- 
equalities are satisfied (with Q 2 = 1/rjJ: Q 2 3> Q 2 (r,b) 3> Aq Cd . That is, the dipole is small 
not only as compared to the typical scale for non-perturbative physics and transverse inhomo- 
geneity in the hadron, namely I /Aqcd, but also as compared to the shorter scale 1/Q s (r,b±), 
which is the local saturation length. This implies that the dipole is only weakly interacting with 
the hadron: S T (x±,y±) ~ 1, or J\f T (r± , b± ) <C 1. But this does not mean that we are a priori 
allowed to linearize eq. ( |3.4j ) with respect to J\f T . Indeed, the r.h.s. of this equation involves an 
integral over all z±, so the virtual dipoles with transverse coordinates (x±, z±) or (z±,y±) can be 
arbitrarily large. In fact, we shall see below that the dominant contribution comes nevertheless 
from z± which is relatively close to in the sense of eq. (|3.5|) , but the upper limit 1/Q s (t, b±) 
in this equation is a consequence of the non-linear effects. 

To see this, it is convenient to divide the integral over z± in eq. ( p.4[ ) into two domains 
("short-range" and "long-range"): 

(A) \z± - b±\ < l/Q s (T,b±) , (B) 1/Q s (r,b±) <C \z±-b±\. (3.6) 
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It is straightforward to compute the contribution of domain (B) to the r.h.s. of eq. ( |3.4[ ): In this 
range, — z±\ ~ \z± — ~ \z± — 6j_|, so the virtual dipoles are both relatively large, and 
therefore strongly absorbed. Thus, to estimate their contribution, one can set S r (x±,z±) 
and S T (z±,y±) 



0, and approximate the ("dipole" [24, 25, [T 



(x ± - z ± ) 2 (y ± - z ± ) : 
This gives (with u\ = (z± 



-^-S T (x±,y±) 



(B) 



b±) 2 > l/Q 2 s (r,b ± )) : 
-a s r 2 L S T (x±,y ± ) / 

-Y(»io2(r,6±))«r(xi.,»±; 



kernel in the BK equation as: 

(3.7) 



du 2 . 



(3.8) 



Since S T (x± , y± ) is of order one for the small dipole of interest, we deduce the following order- 
of-magnitude estimate (which we write for dM T /dr, for further convenience) : 



(B) 



a s rj_Q 2 s (T, b±). 



(3.9) 



Note that, even for this "long range" contribution, the integral in eq. ( |3.8[) is dominated by 
points u± which are relatively close to the lower limit 1/Q s (t, b±) ; this is so because the dipole 
kernel ( |3.7| ) is rapidly decreasing at large distances | z± — b± \ S> r± . 

To evaluate the corresponding contribution of domain (A), we note first that, in this domain, 
all the dipoles are small, so the scattering amplitude is small, M T <C 1, for any of them. It is 
therefore appropriate to linearize the r.h.s. of eq. (|3.4|) with respect to N T (below, u± = x± — z±): 



d_ 



K(r±,b ± ) 



-a. 



(A) 



L /«» d 2 U ± 



2ir U]_(rj_ — u±) 



x \ N T (r ± ,b ± ) -Af T (u±,b± 



u± - r ± 



(3.10) 



2 / 



This is recognized as the BFKL equation in the coordinate representation. Since both u± and r± 
are small as compared to 1/Q s (t, b±), and therefore much smaller than 1/Aqcd, it is appropriate 
to neglect the hadron inhomogeneity when evaluating the r.h.s. of this equation. That is, all 
the functions M T in the r.h.s can be evaluated at the same impact parameter, namely b±. 

To obtain an order-of-magnitude estimate for the r.h.s. of eq. ( |3,10 ), we need an estimate 
for the function M T ( r ±,b±) in the regime where r± <C 1/Q s (t, b±). An approximate solution 
valid in this regime will be constructed in Sect. ||. But for the present purposes, we do not need 
all the details of this solution. Rather, it is enough to use the following "scaling approximation" 



K(r±,b±) 



rlQ 2 s (r,b ± ) 



(3.11) 



with A < 1. In Ref. [29j1 , this approximation has been justified for a homogeneous hadron (no 
dependence upon b±). In Sect. 5 below, we shall find that, in the regime of interest, geometric 
scaling remains true also in the presence of inhomogeneity. 
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The highest value A = 1 corresponds to the "double logarithmic regime"[] in which the dipole 
is extremely small, ln(Q 2 /Q 2 (t, b±)) S> 1, or, equivalently, its impact parameter b± is far outside 
the black disk, b ^> R(t,Q 2 ). Here, however, we are mostly interested in points b± which are 
not so far away from the black disk, since we would like to study how the latter expands by 
incorporating points from the grey area. In this regime, i.e., for b > R(t,Q 2 ) but such that 
(b — R(t, Q 2 )) / R(t, Q 2 ) <C 1, the scattering amplitude is given by eq. (|3.11 ) with a power A 
which is strictly smaller than one (see Sect. 5). 

To simplify the evaluation of eq. ( 3.10| ), we shall divide domain (A) in two subdomains, in 
which further approximations are possible: (A.I) one of the two virtual dipoles, say u±, is much 
smaller than the other one^| : u± <C r±; (A.II) both virtual dipoles are larger than the original 
one, although still smaller than one saturation length: r± <C «i ~ \u± — r±\ <C 1/Q s (t, b±). 

In domain (A.I), the first and third "dipoles" in the r.h.s. of eq. ( |3.10| ) cancel each other 
(since r\ ~ (r± — u±) 2 ), and we are left with 

d 2 u± 



d_ 



K(r±,b±) 



2a, 



(A.I) 



1 



2vr u\ 



u 2 ± Q 2 s (r,b±)) ~ a s (r 2 ± Q 2 (r,b ± ) 



(3.12) 



where the factor of 2 takes into account that one can choose any of the two virtual dipoles as 
the small one. 

In domain (A.II), one can neglect M T (r±, b±) <C M T (u±, b±) ~ N T (u± — r_i_, b±), and obtain: 



-^K(r±,b A 



2a, 



(A.II) 



ITT — {U ± Q s {T,b A 



a, 



(r 2 ± Q 2 (r,b ± ) 



A l-(r 2 ± Q 2 s (r,b ± )) 
1-A 



l-A 



(3.13) 



which is of the same order as the (A.I)-contribution ( p.!2| ) when A < 1, but is logarithmically 
enhanced over it, and also over the long-range contribution (|3.9|), when A — > 1 : 



-^K(r±,b A 



a. 



Q 2 s (r,h 



In 



Q 2 s {rM 



when A = 1. 



(3.14) 



(A.II) 

By comparing eqs. ( |3.9| ) and ( p.l2| )-( ^.14 ), it should be clear by now that, for any A < 1, 
the short-range contribution, domain (A), dominates over the long-range one, domain (B). In 
other words, from the analysis of the non- linear BK equation, we have found that, for a "small" 
incoming dipole, the dominant contribution to the quantum evolution o£ N T (r±,b±) comes 
from still "small" virtual dipoles {(x±, z±), (z±, y±)} (see Figure ||). This has two important 
consequences, (i) One can linearize the BK equation with respect to M T , as we did already in 
eq. ( 3.10 ). This gives the BFKL equation, (m) One can ignore the transverse inhomogeneity in 
the BFKL equation. That is, one can replace eq. ( |3.10| ) by 



— Af T (r±,b±) 



a. 



d 2 z ± 



2tt (z± - x±) 2 (z± - y±_Y 



:{N T (z ± - x ± ,b ± ) +M T (z ± -y±,b±) - Af T (r±,b±)}, (3.15) 



3 Strictly speaking, there is no geometric scaling in this regime, but for power counting purposes one can assume 



Indeed, at very high Q , the giuon distribution 



that Af T {r±, b±) oc r ± xG(x, l/r ± , b±) (cf. eq. (2.4)) is linear in r\ 
xG(x, Q 2 ,b±_) is only weakly depen dent upon Q . 

6 That is, in the notations of eq. (^4), the point z± is within the area occupied by the original dipole (x±,y±), 
and much closer to x± than to y±. 
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in which all the amplitudes N T are evaluated at the same impact parameter, namely at b±. 
Note that, as compared to eq. ( 3.13] ), there is no need to insert an upper cutoff ~ 1/Q S in 



the integral in eq. ( 3.15| ). This is so since, to the accuracy of interest, the solution J\f T (r±,b±) 



to eq. ( 3. 15| ) is actually insensitive to such a cutoff. One can understand this on the basis of 



eq. ( p. 13|) : For large r± (with r 2 ± Q 2 s (r, b±) <C 1 though), the solution has the "scaling" behaviour 
in eq. (]3.1l| ) with a power A which is strictly smaller than one (see the discussion in Sects. 4 and 



5). Then the integral in eq. ( 3.13 ) is dominated by points u± which are close to the lower limit 
rj_, i.e., by virtual dipoles which are not much larger than the incoming dipole. The dependence 
upon the upper cutoff ~ 1/Q S is therefore a subleading effect, which can be safely ignored. 



We finally come to the last step in our argument: Since the dependence of eq. ( 3.15| ) upon 



b± is only "parametric", it is clear that the impact parameter dependence of the solution is 
entirely fixed by the initial condition. This, together with eq. ( |2.8| ), implies that M T (r±, b±) has 
the factorized structure: 

K(r±,b ± ) = S(b ± )Af T {r ± ) (3.16) 

where S(b±) is the transverse profile of the initial condition, while M T {r±) satisfies the homo- 
geneous BFKL equation and will be discussed in Sect. ||. A brief inspection of the previous 
arguments reveals that the terms neglected in our approximations are suppressed by either 
powers of r^Q^r, b±) (e.g., the long-range contribution (|3.9|), or the cutoff-dependent term in 
eq. fl3.13|) ), or by powers of Aq CD /Q 2 (t, b±) (the inhomogeneous effects in eq. ( 3.10 )). This 



specifies the accuracy of the factorized approximation in eq. (3.16). 



3.2 More on the saturated gluons 

In the previous subsection, we have seen that colour sources located far away from the impact 
parameter of the dipole, such as | — 6j_ | > 1/Q s (t, b±), do not significantly contribute to 
the scattering amplitude in the grey area. In what follows, we shall examine more carefully a 
particular contribution of this type, namely, that associated with the saturated gluons within 



the black disk: \z±\ < R(t,Q 2 ). Indeed, it has been recently argued 14] that, by itself, this 
contribution would lead to unitarity violations. To clarify this point, we shall compute this 
contribution within the effective theory for the Colour Glass Condensate, where the physical 
interpretation of the result is transparent. The same result will be then reobtained from the 
BK equation. Our analysis will confirm that, for impact parameters b± within the grey area, 
this long-range contribution is indeed subleading, and can be safely neglected at high energies. 



As we shall see, non-unitary contributions of the type discussed in Ref. [14| appear only in 
the physically uninteresting case where the exernal probe carries a non-zero colour charge (as 
opposed to the colourless dipole). 

According to eqs. ( |3. 1[ ) and (|3.2| ), the scattering amplitude at rapidity r depends upon the 
Coulomb field a^(x±) at all the space-time rapidities rj < r. In general, this is related to the 
colour sources in the hadron via the two-dimensional Poisson equation — V^_a^(xj_) = p^(x±), 
with the solution: 

<4(x±) = J cPz±(x ± \-^ T \z ± )p°(z ± ). (3.17) 

In this equation, p^(x±) is the colour charge density (per unit transverse area per unit space- 
time rapidity) of the colour sources at space-time rapidity r]. Here, we are only interested in 
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such colour sources which are saturated. To isolate their contribution, it is important to remark 
that these sources have been generated by the quantum evolution up to a "time" equal to 77, 
so the corresponding saturation scale is Q s {r],b±) (and not Q s (t, b±)). Thus, the integration in 
eq. (|3.17 ) must be restricted to \z±\ < R(f], Q 2 ), with R(rj, Q 2 ) the black disk radius at rapidity 



77 (cf. eq. Q2.10D ), and Q 2 the typical momentum carried by the Fourier modes of a^(x±) (as 
usual, this is fixed by the transverse size of the incoming dipole). 

There is a similar restriction on the values of 77 : For given Q 2 , there is a minimum rapidity 
f\{Q 2 ) below which there is no black disk at all: R(tj,Q 2 ) = for 77 < fi(Q 2 )- This is the 
rapidity at which the black disk first emerges at the center of the hadron, namely, at which 
Q 2 (fi,b = 0) = Q 2 (see Sect. 4.4 below). Thus, in order to count saturated sources only, the 
integral over 77 in Wilson lines like ([T^) must be restricted to the interval fi(Q 2 ) < 77 < t. In 
Fig. ^, the saturated sources (with momentum Q 2 ) occupy the lower right corner, below the 
dashed line which represents the profile of the black disk as a function of 77. 

The external point x± in eq. ( |3.17 ) is at the impact parameter of the quark (or the antiquark) 



in the dipole, so it satisfies \x±\ 3> R(t, Q 2 ) > R(rj, Q 2 ) for any 77 < r. It is therefore appropriate 
to evaluate the field ( |3.17| ) in a multipolar expansion: 

1 r R d l f R 

o%(x±) = (xx\—^r\0x) J d 2 z ±P ^(z ± ) - ^— (x ± \— ^2-|0±>_/ d 2 z ± z l p^z ± ) +••■ 

= (zjJ-Ij-IOj.) Q a + ^rK +..-, (3.18) 

where R = R(rj, Q 2 ), Q a is the total colour charge within the black disk, T> l a is the corresponding 
dipolar moment, etc. 

To compute the scattering amplitude (|3.1|), one has to construct the Wilson lines and V 
with the field ( 3.18j ) and then average over a (or, equivalently, over p) as in eq. (|3.3| ). In what 



follows, it is more convenient to work with the probability distribution for p, i.e., W r [p]. In 
general, this distribution is determined by a complicated functional evolution equation, which 
is very non- linear g jj]}. 



However, as observed in Ref. [22], this equation simplifies drastically in the saturation 
regime, where the corresponding solution W T [/o] is essentially a Gaussian in p. This can be 
understood as follows: The non-linear effects in the quantum evolution enter via Wilson lines 
like eq. ( |3.2| ). At saturation, the field a in the exponential carries typical momenta k± ~ Q s and 
has a large amplitude a ~ 1/g. Thus, the Wilson lines are strongly varying over a transverse 
distance 1/Q s (r, b±). When observed by a probe with transverse resolution Q 2 <C Q 2 (i~,b±), 
these Wilson lines are rapidly oscillating and average to zero. Thus, at saturation, one can 
drop out the Wilson lines, and all the associated non-local and non-linear effects. Then, the 
probability distribution W T [/9] becomes indeed a Gaussian, which by the same argument is local 
in colour and space-time rapidity, and also homogeneous in all the (longitudinal and transverse) 
coordinates. The only remaining correlations are those in the transverse plane, and, importantly, 



these are such as to ensure colour neutrality [18] 



Specifically, the only non-trivial correlation function of the saturated sources is the two-point 
function, which reads @| (see also Sect. 5.4 in Ref. flgl ) : 

(p a v (z ± )p b v ,(u ± )) T = S ab 6(r ] - V ')X(z ± -u ± ), 
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1\ T\ X 

Figure 2: The longitudinal profile of the hadron as it appears in a scattering at given r and Q 2 . 
The longitudinal coordinate is on the horizontal axis, and is measured in units of space-time 
rapidity. The tranverse coordinate is on the vertical axis. A longitudinal layer at rapidity rj is 
delimited for more clarity. The wavy line represents the colour field a v created at point x± by 
the (saturated) source p v at z±. The enclined dashed line represents the limit of the black disk, 
which increases linearly with 77, as we shall see in Sect. 4. 
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A(A:_l) = -k\. (3.19) 



7T 



For given rj and k\ ~ Q 2 , eq. ( |3.19| ) holds for points z± and u± within the black disk of radius 
R(rj,Q 2 ). The crucial property of the 2-point function ( |3.19| ) is that it vanishes as k\ — > 0. 
Physically, this means that, globally, the saturated gluons are colour neutral^, as anticipated: 
(Q a Q a ) = 0, where Q a is the total colour charge (at given rj) in the transverse plane. In fact, 
since the Wilson lines average to zero over distances Ab± ^ 1/Q s (rj,b±), it follows that colour 
neutrality is achieved already over a transverse scale of the order of the saturation length: 

/ d 2 x ± [ d 2 y ± {p a Jx L )p a Jy L )) = 0, (3.20) 
JASi JAS± 

where AS± is, e.g., a disk of radius R > 1/Q s (r],b±) centered at b±. 

This immediately implies that, as soon as the black disk is large enough, the overall charge of 
the saturated gluons vanishes, Q a = 0, so we can ignore the monopole field in eq. ( |3.18| ). Here, 
"large enough" means, e.g., R(t, Q 2 ) 3> 1/Q, which guarantees that the radius of the black disk 
is larger than the saturation length 1/Q s {rj,b±) at any b± within the disk and at any r] in the 
interval f\(Q 2 ) < f] < r (since then Q s (rj,b±) > Q). In these conditions, the dominant field of 
the saturated gluons at large distances is the dipolar field in eq. ( 3.18| ). 

Below, we shall need the two-point function of this field: 

(a a v (x ± )a^(y ± )) T = 6 ab 6( V - rj) 1t ,(x ± - y ± ). (3.21) 
From eq. ( |3.1§| ), we deduce (with colour indices omitted, since trivial): 

with (cf. eq. Q)): 



pR rR 

(V l V 3 ) = / d 2 z ± d 2 u± z % u 3 \{z ± - u± 



j2 r R j2 f Q ° d 2 k ± ikl , Zl _ Ul) d 2 k\ 
dZ ^J dV i WdkIT 
~ 8 ij 2R 2 (ri,Q 2 ) . (3.23) 

where R = R(rj,Q 2 ), Q s = Q s (Vjb±) (with b± = (z± + u±)/2 < R), and formal manipulations 
like integrations by parts or the use of the Fourier representation of the 5-function were permitted 
since RQ S 3> 1. Thus, finally, 



lv(x ± ,y ± ) = -L R 2 (v,Q 2 ), (3.24) 



which, we recall, is valid only as long as x±, y± 3> R(r], Q 2 ). 

We are now in a position to compute the scattering amplitude ( |3.1| ) for the scattering between 
the incoming dipole and the dipolar colour charge distribution within the black disk. To this 



7 Since the distribution of p is a Gaussian, the fact that {Q a Q a ) vanishes is equivalent to Q a = 0, which means 
colour neutrality indeed. 
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aim, we have to average the product tr(V"f (x±)V(y±)) over the Gaussian random variable a v 
with two-point function ( |3.21 ). The result of this calculation is well-known (see, e.g., [|H|): 



S T (x±,y±) = exp 



9 2 C F 



drj 7r,(x_ L ,X_ L ) +7r?(y±,y±) - 2lrj{x±,y±) 



(3.25) 



where C F = t a t a = (N 2 - l)/2N c and the lower limit r in the integral is a shorthand for 
fi(Q 2 ) with Q 2 = The "tilde" symbol on S is to remind that this is not the total S- 

matrix element, but just the particular contribution to it coming from the saturated gluons. An 
immediate calculation using eqs. ( |3T24|) and flOp yields (with a s = 2a s C F /'K ~ a s N c /ir in the 
large N c limit) : 



S T (x±,y±) = exp<-a s 



(x± - y±) 
2xlyl 



d v R 2 ( v ,Q 2 ) 



* exp|-a s ^ J_ drj R 2 ( V ,Q 2 )j, 
where in the second line we have replaced in the denominator x\ 



(3.26) 



b\ (which is appropriate 



since b± S> R(t,Q 2 ) 3> r±). The exponent in eq. ( 3.26) ) vanishes when the dipole shrinks to a 



point, rj_ — > 0. This is the expected dipole cancellation, manifest already on eq. ( |3. 25 ). 

The previous derivation makes the physical interpretation of eq. ( 3.26 ) very clear: The 
exponent there is the square of the potential ~ gt a r_\_jpR(rj) for the interaction between two 
dipoles — the "external dipole" of size r± and the dipole made of the saturated gluons (at a 
given space-time rapidity 77), with size R{rj) — separated by a large distance b. There exists one 



layer of saturated gluons at any 77 within the interval f\{Q 2 ) < r/ < r, so eq. ( 3.26 ) involves an 
integral over this interval. 

Of course, eq. ( |3.26 ) can be also obtained directly from the BK equation, although, in that 
context, its physical interpretation in terms of dipole-dipole scattering may not be so obvious. 
In fact, this is just a particular piece of what we have called "the contribution (B)" in the 
previous subsection, i.e., the contribution of the points z± satisfying \z± — b±\ » l/Q s (r, b± ) . 
If z± now refers to the saturated gluons within the black disk, then it is further restricted by 
\z±\ < R(t,Q 2 ), which implies \z± — b±\ rs b for the same reasons as above. Then, a simple 
calculation similar to eq. ( |3.8[ ) immediately yields 



Q-S T {x±,y±) 



—a. 



b\ 



R 2 (t, Q 2 )S T (x ± ,y±) 



(3.27) 



which after integration over r is indeed equivalent to eq. ( |3.26D f\ 

For comparison with eq. ( |3.26j ), it is interesting to compute also the 5-matrix element for 
a coloured external probe, e.g., a quark, which scatters off the saturated gluons in the eikonal 
approximation. A calculation entirely similar to that leading to eq. ( |3.26 ) yields (b± is the 



used in Refs. |23, |18| to derive eq 
currently interested in 



The mismatch by a factor of two betw een eqs. (3.26) and (3.27) is inherent to the mean field approximation 
{22] , [lit to derive eq. (3.2E), and is completely irrelevant for the kind of estimates that we are 
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transverse location of the quark) : 



^( tr ^( 6 J-)> T = exp 1-^^^7,(6^,6 



1 



exp i -a s ^ J_ d V R 2 (r], Q 1 ) \ (3.28) 



where the second line follows after using eq. ( |3.24 ) 



To summarize, the amplitude for the scattering off the saturated gluons decreases like 1/6* 
for an external dipole, but only as l/6j_ for a coloured probe. This difference turns out to be 
essential: because of it, this long-range scattering plays only a marginal role for the dipole, while 
it leads to unitarity violations in the case of the coloured probe (although the very question of 
unitarization makes little physical sense for a "probe" which is not a colour singlet). 

To see this, assume the long-range contributions shown above to be the only contributions, 
or, in any case, those which give the dominant contribution to the cross-section. Then, one 
can rely on the previous formulae to estimate the rate of expansion of the black disk. Namely, 
assume that, for the purpose of getting an order-of-magnitude estimate, one can extrapolate 
eqs. ( gj2g ) and ( g2p up to energies where the black disk approaches the incidence point b± of 
the external probe. Then, we expect the exponents in these equations to become of order one 
for b ~ R(t, Q 2 ). For the external dipole, this condition implies: 



r 2 



R ^ Q 2 ) J. dr l R (V>Q ) ~ 1- (3-29) 

This gives (recall that Q 2 = 1/tj_): 

Q 4 R\t, Q 2 ) = a s J T _ dT] Q 2 R 2 (v, Q 2 ), (3.30) 

or, after taking a derivative w.r.t. r , 

2^-(q 2 R 2 (t,Q 2 )) = a 8 , (3.31) 

whose solution R 2 (t,Q 2 ) increases linearly with r. 

By contrast, for a coloured probe, the same condition yields: 

a s J T _ dT] R 2 (r], Q 2 ) = R 2 (t, Q 2 ) (3.32) 

or after taking a derivative w.r.t. r: 

^-R 2 (r,Q 2 ) = a s R 2 (r,Q 2 ), (3.33) 
ar 



which gives an exponential increase with r, as found in Ref. [14]. 

We thus see that the violation of unitarity by long-range Coulomb scattering reported by the 
authors of Ref. H] is related to their use of an external probe which carries a non-zero colour 
charge. This case is physically ill defined, and therefore uninteresting (note, indeed, that (trV*) is 
not a gauge-invariant quantity); in particular, its relevance for the problem of gluon saturation in the 
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target wavefunction remains unclear to us (since the relation between "blackness" and saturation 
holds only for dipole probes; cf. the discussion prior to eq. (|2.9| )). 

On the other hand, for the physically interesting case of an external dipole, the contribution 
( 3.31] ) to the expansion of the black disk not only is consistent with unitarity — if this was the only 
contribution, the cross-section oc R 2 (t, Q 2 ) would increase linearly with r — , but at large r, is even 
negligible as compared to the corresponding contribution of the short-range scattering (which gives 
a cross-section increasing like r 2 , cf. eqs. (|2.16| )- (|2.17| )). 

These considerations are conveniently summarized in the following, schematic, approximation to 
the scattering amplitude in the grey area, which follows from the previous analysis in this section: 
M T (r±,b±) is the sum of two contributions, a short-range contribution, cf. eq. (2.15), and a long- 
range one, cf. eq. fl3.26|) : 

M T (r ± ,b ± ) « W^e^e" 2 ^ 6 + a, A- [ d n R 2 ( V ,Q 2 ) , (3.34) 

with the short-range contribution determined by the solution to the homogeneous BFKL equation 
( |3.15| ) together with the assumed exponential fall-off of the initial condition (see Sects. 4 and 5 
below for more details), and the long-range contribution obtained by keeping only the lowest-order 
term in eq. ( P-26| ) (which is enough since we are in a regime where S T (x±,y±) ~ 1). At the 
initial rapidity To, the long-range contribution vanishes, while the short-range contribution reduces 
t0 e -2m.(6-R(r ,Q 2 )) wjth r( Tq ,Q 2 ) given by eq. pig 1 ), as it shouldfl 

For given r and r_|_ = l/Q, eq. ( 3.34| ) applies at impact parameters in the grey area, R(t,Q 2 ) < 
b < Rh(t), but it can be extrapolated to estimate the boundaries of this area, according to eqs. ( |2.9|) - 
(|2.12j ). It is easy to check that, at high energy, both these boundaries are determined by the short- 
range contribution, i.e., the first term in the r.h.s. of eq. ( |3.34j ). Thus, this contribution dominates 
the scattering amplitude at any b± in the grey area. By comparison, the long-range contribution is 
suppressed by one power of l/a s T. 

For instance, the black disk radius is obtained by requiring M T (r±,b±) ~ 1 for b± ~ R(t,Q 2 ) 
(cf. eq. (^9|)). If one assumes the short-range contribution to dominate in this regime, one obtains 
the estimate ( |2.16 ) for the black disk radius: R(t,Q 2 ) ~ (w/2m T )a s T. By using this result, one 
can evaluate the corresponding long-range contribution, and thus check that this is comparatively 
small, as it should for consistency with the original assumption: 



< 1. (3.35) 



^R^ml dr > R2 M 2 



r I mt 



«(r,Q 2 )~ a s T/ mjr CH S T 

By contrast, if one starts by assuming that the long-range contribution dominates, then one is 
running into a contradiction, since in this case R(t, Q 2 ) oc ^cx s t, cf. eq. (|3.31j ), and the short-range 
contribution increases exponentially along the "trajectory" b = R(t,Q 2 ). 

A similar conclusion holds for b = Rh{t) [since Rh{t) = R{t,Q 2 = Aq C£) ), cf. eq. ( |2.12|) 1, 
and therefore for any point b± within the grey area. (The hadron radius Rh{t) will be evaluated in 
Sect. 4.2 below.) Thus, the short-range contribution is indeed the dominant one in the grey area. 
This contribution preserves the exponential fall-off of the initial condition, and therefore saturates 
the Froissart bound, as explained in Sect. 2. 

9 That is, the initial scattering amplitude N T ,A r l, bx) is equal to one within the black disk < R(to,Q 2 )), 
and decreases exponentially outside it. Eq. ( 3.34 ) applies, of course, only at impact parameters outside the black 
disk. 
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It is essential for the consistency of the previous arguments — which combine perturbative 
quantum evolution with non-perturbative initial conditions — that perturbation theory has been 
applied only in the regime where it is expected to be valid, namely, in the central region at b < Rh{t), 
where the gluon density is high and the local saturation scale Q s (r,b) is much larger than Aqcd- 
It has been enough to consider this region for the present purposes since this includes both the 
black disk and the grey area which controls its expansion. Within this region, the perturbative 
evolution equations of Refs. 0, ||, ^, |l0| can be trusted, and the additional approximations that we 
have performed on these equations are under control as well. When supplemented with appropriate 
boundary conditions — which are truly non-perturbative, since reflecting the physics of confinement 
— , these equations allow one to compute the rate for the expansion of the black disk, and, more 
generally, to follow this expansion as long as the black disk remains confined within the region of 
applicability of perturbation theory (which includes, at least, the central area at the initial rapidity 
t : b < R h (t )). 

But if one attempts to follow this expansion up to much higher rapidities, where R(t,Q 2 ) 3> 
Rh{ t o), then a strict application of the perturbative evolution with initial conditions at r = tq 
would run into difficulties^]. The difficulties arise since, in the absence of confinement, the long- 
range dipolar tails created by the saturated gluons can extend to arbitrary large distances, and thus 
contribute to scattering even at very large impact parameters (b ^> Rh{tq)), where physically there 
should be no contribution at all. This is illustrated by eq. ( P-34| ): We have previously argued that, 
for b± within the grey area, the dominant contribution comes from short-range scattering (the first 
term in the r.h.s. of eq. ( p.34| )). But if one extrapolates this formula at very large b 3> Rh{t), 
then, clearly, the long-range contribution, which has only a power-law fall-off with b±, will eventually 
dominate over the short-range contribution, which decreases exponentially. When this happens, 
however, the impact parameters are so large (b — R(t,Q 2 ) ^> that the long-range scattering 

is controlled by the exchange of very soft (k± ^ m n ) quanta, which in a full theory would be 
suppressed by the confinement. That is, in a more complete theory which would include the physics 
of confinement, the long-range contribution to eq. ( 3.34 ) would be suppressed at very large b± by 
an additional factor e - 2 "M6--R(T,Q )) ^ SQ tnat tne sn ort-range contribution will always dominate, for 
all impact parameters. But in the present, perturbative, setting, the only way to avoid unphysical 
long-range contributions is to start the quantum evolution directly in the grey area, as we did before 
(rather than try to construct this grey area via perturbative evolution from earlier rapidities tq <C r, 
at which the points b± of interest were far outside the initial grey area: b ^> Rh(tq)). 

Note finally that what is truly remarkable, and also essential for our conclusion on the Froissart 
bound, is not the suppression of the long-range non-perturbative contributions by the confinement 
— this is only to be expected in the full theory, and can be also enforced in the present calculation 
by appropriately chosing the boundary conditions — , but rather the suppression of the long-range 
perturbative contribution within the grey area (where perturbation theory applies, so its predictions 
must be taken at face value). We mean here, of course, the fact that the long-range contribution 
to eq. ( |3.34| ) falls off like l/b 4 ^, and not like \jb\, as it would have been the case if the saturated 
gluons were uncorrelated. The mechanism for this suppression is purely perturbative, and related to 
saturation: The saturated gluons are globally colour neutral, so the monopole fields of the individual 

10 The discussion in this paragraph has been inserted as a partial response to criticism by Kovner and Wiedemann 
[ ^2[ , written in response to the original version of this paper. Please see the note added at the end of the paper 
for further discussion. 
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gluons are replaced at large distances S> 1/Q S by the more rapidly decreasing dipolar field of the 
whole distribution. To understand the relevance of this suppression for the Froissart bound, consider 
what would happen if the saturated sources were statistically independent, i.e., if eq. ( |3.19| ) was 
replaced by 

(a#*i)^(«l)>t = Xd ab 5( V - V ')5^(z ± -u ± ). 

Then the exponent in eq. ( |3.26| ), and also the second term in the r.h.s. of eq. ( [3.34 ), would change 
into 2 

a S ^r [_ drjR 2 ( V ,Q 2 ), 

which would generate a black disk increasing exponentially with r (cf. eq. ( p.32| )-( |3.33] )). Thus, the 
long-range contribution would dominate already within the grey area, and the Froissart bound would 
be violated. We thus conclude that colour correlations at saturation are essential to ensure unitarity. 

4 Black disk evolution and the Froissart bound 



In this section, we shall exploit the factorization property ( 3.161 ) together with the known solution 



M T (t± ) to the homogeneous BFKL equation in order to compute the scattering amplitude in the 
grey area, and thus study the evolution of the black disk with increasing energy. After briefly 
recalling the BFKL solution, in Sect. 4.1, we shall then compute the radius of the black disk 
R(t,Q 2 ) and derive the Froissart bound (in Sect. 4.2). Then, in Sect. 4.3, we shall study the 
impact parameter dependence of the saturation scale Q 2 (t, b±), and deduce a physical picture 
for the expansion of the black disk, to be exposed in Sect. 4.4. 

4.1 Scattering amplitude in the BFKL approximation 

Eq. ( 3.16] ) for the scattering amplitude in the grey area involves the solution M T (r±) to the 



homogeneous BFKL equation, i.e., the BFKL equation without impact parameter dependence. 
This solution is well known, and we shall briefly recall here the relevant formulae, at the level of 
accuracy of the present calculation. (See Refs. |37], |{| for a similar approach and more details.) 
The solution can be expressed as a Mellin transform with respect to the transverse coordinate: 

M T (r ± = l/Q) = [ ^ ( e « s r W (i)-#)^(i-A)} j (41) 



where ip(\) is the di-gamma function, and A is an arbitrary reference scale, of order Aqcd- 
The contour C in the inverse Mellin transform is taken on the left of all the singularities of the 
integrand in the half plane Re A > 0. Note that, since M T (r±) is a function of r^_, we find it 
convenient to use the momentum variable Q 2 = l/r -2 to characterize the transverse resolution 
of the dipole. From now on, we shall again use the notation M T (Q 2 ) = M T (r± = l/Q), which 
was already introduced in Sect. 2 (cf. eq. ( |2.9| )). 

We are interested here in a regime where the energy is very high, 6l s t 3> 1, and the dipole is 
small: Q 2 3> A 2 . In these conditions, it is appropriate to evaluate the integral ( |4.1| ) in the saddle 
point approximation. Higher is the energy, better is justified this approximation, and closer is 
the saddle point Aq — which is a function of (lnQ 2 /A 2 )/a s r — of the so-called "genuine BFKL" 
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saddle-point at Ao = 1/2. (This is the saddle point which governs the asymptotic behaviour of 
the solution to the BFKL equation at very large energy.) In fact, for 

(4.2) 

which is the most interesting regime here, the saddle point is easily estimated as: 



with (5 = 28£(3). In fact, the recent analysis in Ref. [29| shows that eq. (|4.3| ) remains a 
good approximation for the saddle point even for comparatively low energies, such that a s r ~ 
ln(Q 2 / A2 )- This saddle point gives the standard BFKL solution, which, after multiplication 
with the profile function (cf. eq. ( 3.16| )), provides the scattering amplitude in the grey area in 
the present approximation: 

MAQ 2 ,b ± ) ~ S(b ± )J-» y= = exp<;-— ^(ln^l }, (4.4) 




Q 2 ^2irf3a s T 

where u = 4 In 2 is the customary BFKL exponent. The factor ^/27r/3a s r in the denominator 
comes from integrating over the Gaussian fluctuations around the saddle point. When expo- 
nentiated, this gives a contribution oc ln(a s r) which is subleading at large energy and will be 
ignored in what follows. It is then convenient to rewrite eq. (|4,4j) as follows: 



1 . Q 2 if, Q 2 " 2 ' 



where we have also used S(b±) « e- 2 ™^ as appropriate for sufficiently large b (b 3> i?o, cf. the 
discussion after eq. (2.6)). This is the most interesting case here, since we consider the high 



energy regime in which the black disk is already quite large. 

Eq. flL5|) is valid for those values of the parameters r, Q 2 and b± for which our previous 
approximations are justified, namely, such that the conditions A 2 <C Q 2 (t, b) <C Q 2 are satisfied. 
As it was anticipated in Sect. 2, and will be verified below in this section, these conditions 
are realized within a corona at R(t,Q 2 ) <C b± <C Rh(t), which, with increasing energy, moves 
further and further away from the center of the hadron. 



When decreasing b± towards R(r, Q ) at fixed r, or, equivalently, increasing r at fixed b 



the scattering amplitude increases towards one, and the BFKL approximation (|4,5| ) ceases to be 
valid. (The dipole resolution Q 2 is always fixed in these considerations.) But it is nevertheless 
legitimate to use eq. (|4.5|) in order to estimate the boundary of its range of validity, that is, 
the black disk radius R(t,Q 2 ), or the saturation scale Q 2 (t, b±). Indeed, the non-linear effects 



become important when the BFKL solution (4.5) becomes of order one. This condition can 



be written either as an equation for R(t,Q ) for given r and Q , namely, eq. ( |2.9[ ), or as an 
equation for Q 2 (t, b±) for given r and b±, namely, eq. (|2.11 ). (One could, of course, similarly 



introduce and compute also a critical rapidity t(Q 2 , b±) at which blackness is reached for given 
Q 2 and b±, but this is less interesting for our subsequent discussion. See, however, Sect. 4.4.) 
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4.2 The black disk radius 



In this subsection, we shall use eqs. (|2,9|) and (|4.5|) to compute the radius of the black disk and 
study some limiting cases. Eq. ( |2.9| ) amounts to the condition that the exponent in eq. 
vanishes^, which immediately implies: 



1 . Q 2 if, Q 2X 2 



2m^(r, Q 2 ) = ua s r -.^- — ^J . (4.6) 
The right hand side is positive as long as Q 2 < A 2 e cclsT = Q 2 (t, 6 = 0), with 

P , 1 



c=-^ + -y/P(J3 + 8u;) =AM.... ( 47 ) 

As anticipated by our notations, 

Ql(r,b= 0) = A 2 e c " sT (4.8) 

is the saturation scale at the center of the hadron (this will be verified via a direct computation 
in the next subsection). This is as expected: for Q 2 > Q 2 (r,b = 0), the hadron looks grey 
everywhere, so R(t, Q 2 ) = 0. 

The other extreme situation is when Q 2 ~ A 2 , so that the black disk extends up to the edge 
of the hadron (cf. eq. ( 2.12j )). Equation (|4.6| ) yields then: 



Rh{t) « ^r, (4.9) 

which should be seen only as a crude estimate: for such a small Q 2 , our approximations are not 
justified any longer. 

But the physically interesting case is when Q 2 3> A 2 , but the energy is so large that the 



condition (^2) is satisfied. Then one can neglect the term quadratic in InQ /A in eq. 



(since this term vanishes when r — > oo), and deduce that: 

The term linear in InQ 2 /A 2 , although subleading at large r (since independent of r), has been 
nevertheless kept in the above equation since, first, we expect this term to give the dominant 
Q 2 -dependence of the cross-section at high energy, and, second, it measures the separation 
between the black disk and the edge of the hadron in the high energy regime. Specifically: 

R H (r) - R(t,Q 2 ) « ^ln|J, (4.11) 

which is fixed (i.e., independent of r), but large when Q 2 3> A 2 . This is important since 
the points b± at which our approximations are justified should lie deeply within this corona: 
R(t,Q 2 ) « i) « Rh(t). Thus, as anticipated in Sect. 2, the fact that Q 2 S> A 2 ensures the 
existence of a large grey area in which our approximations apply. 



At the level of the present approximation, one can take k — 1 in eqs. (|2.9D and (2.11) without loss of accuracy. 
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Equation ( 4.1G| ) is our main result in this paper. It shows that, at very high energy, the 



radius of the black disk increases only linearly with r, i.e., logarithmically with the energy. This 
is the result anticipated in eq. ( |2.16| ). The corresponding cross-section is given by eq. ( |2.i7| ) 
and saturates the Froissart bound, that is, it grows like In 2 s, with a proportionality coefficient 
which is universal (i.e., the same for any hadronic target), and which reflects the combined role 
of perturbative and non-perturbative physics in controlling the asymptotic behaviour at high 
energy. 

In the remaining part of this paper, we shall further explore this result and gain a different 
perspective over it by computing also the saturation scale and studying the geometric scaling 
properties. 

4.3 Saturation scale with the impact parameter dependence 

In this subsection, we shall compute the saturation scale for an inhomogeneous hadron and 
study its variation with the energy and the impact parameter. Previous studies of this kind 
were restricted to a homogeneous hadron El M , but, as we shall see, the dependence upon 



the impact parameter introduces some interesting new features. 

By inspection of eq. (|4.5|), it is clear that the saturation condition ( |2.11| ) amounts to the 
following, second order algebraic equation for the quantity p s = (l/a s r) In Q 2 /A 2 : 

P l + (3 Ps -2(3u J = -2P^. (4.12) 

a s T 

The solution to this equation and the corresponding saturation scale read: 




Ps(r,b) = + + — (1-^- , (4.13) 

Q 2 s (r,b) = A 2 e aaTp ^ T ' b) . (4.14) 

Note that, in general, the impact parameter dependence in the saturation scale ( [4.14| ) is not 
factorizable. Below, however, we shall recover factorization in some specific limits. 

It can be easily checked that the above equations ( 4.13| )-( 4.14 ) and eq. ( [4.6| ) are consistent 



with each other, in the sense that Q 2 (r,6 = R(t,Q 2 )) = Q 2 , as it should (cf. eq. (pTOP). 
In particular, one can use eqs. fl4.13| )-( 4.14|) to rederive the results in eqs. (|4.7|) - (|4.8[) for the 



saturation scale Q 2 (r,b = 0) at the center of the hadron, as well as eq. ( |4.9| ) for the hadron 
radius. 

A pictorial representation of the 6-dependence of the saturation scale, as emerging from 
eqs. ( 4.13|) — ( 4.14 ), is given in Fig. ||[ As compared to eq. ( 4.13|) , in this graphical representation 



we have replaced 2m 7r b — > — In S(b), with S(b) given by a Woods-Saxon profile, cf. eqs. fl2.G| ) 
(|2.7| ); this is more realistic than the exponential at short distances, b ^ Rq, where it has a much 
slower decrease, but it shows the expected fall-off S(b) ~ e -"irn w (b-R Q ) ^ larger distances. As 
manifest on this figure, Q 2 (t, b) is itself very similar to an exponential for all distances b Rq. 
This can be understood via a further study of eq. (|4.13| ), which will also reveal that, in fact, there 
is a change in the slope of the exponential with increasing b : To a very good approximation, 
the plot in Fig. |3]can be seen as the superposition of two exponentials, one at small b, the other 
one at large b, which have different exponential slopes. 
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Figure 3: The saturation scale Q s {b)/Q s {b 
Woods-Saxon profile function of eqs. (|2,6|) -(| 
distance is measured in units of 1/2711,7,-. 



0) from eqs. (4.13)— (4.14) for a s T = 3 and the 
with Rq = 3/2m n . On the abscisa, the radial 



To see this, note that the function 

A(r,6) = 1 - = 1 - — -- , 4.15 

Loa s T Rh{t) 

which enters the square root in eq. ( |4. 13|) , is positive semi-definite for b within the hadron radius 
[b < Rh(t)), and monotonically decreasing with b from A(r, b = 0) = 1 to A(r, b = Rh(t)) = 0. 
This suggest two different approximations according to whether A is close to one (for b sufficiently 
small) or close to zero (for b sufficiently close to Rh(t)). (Note that the factor 8co/f3 multiplying 
A(t, b) in eq. (4.13) is a number of order one, 8w//3 ~ 0.67, so it does not interfere with our 
order-of-magnitude estimates.) 

(I) If A is close to one, which happens when b is much smaller than the hadron radius: 

— < 1, or b < Rh(t), (4.16) 

U)Ot s T 

one can evaluate the square root in eq. ( [4.13| ) in an expansion in powers of 1 — A (this is 
equivalent to an expansion of p s (r,b) in powers of b around p s (r, b = 0) = c, cf. eq. (|4.7|) ). To 
linear order in this expansion, one obtains: 



ft(T ,») * c- X, = . /£+»i = 0.644..., (4.17) 

X s a s T y 4/? 

which gives (with 7 = 1/A S ~ 1.55): 

Q 2 s {t, b) ~ A 2 e c ^ T e ~ 2 ^ b = Q 2 (r, b = 0) [5(6)] 7 . (4.18) 
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This is in a factorized form, although, as compared to the corresponding factorized structure 
of the scattering amplitude ( |3.16| ), it features some "anomalous dimension" 7 for the profile 
function. The value Q 2 (t, b = 0) at the center of the hadron is the same as the saturation scale 
for a homogeneous hadron previously found in Refs. 21, 2£, 3C]. Also, the constant A s which 
appears in eq. ( |4.17D is the value of the saddle point Ao in the Mellin representation ( f4.1| ) for 
Q 2 = Q 2 s (r,b = 0) (i.e., eq. (£J) with ln(Q 2 /A 2 ) = ca s r) @. 

Equation ( |4.1§| ) shows that the saturation scale decreases exponentially with the distance b 
from the center of the hadron, with a typical decay scale ~ 1/27771^. (Of course, this exponential 
law applies only for values b which are not too close to the center, 6 <; R , cf. Fig. ||.) 
(II) When b is sufficiently close to Rjj(t), in the sense that: 

Rh(t) 



AM) = 



Rh(t) 



« 1, 



(4.19) 



than one can expand eq. ( |4.13j ) in powers of A (this is an expansion of p s (r, b) around p s (t, b 
Rh) = 0). To lowest order in this expansion, one obtains p s (r, b) ~ 2o;A(r, b), and therefore: 



Q 2 s (r,b) 



jy2 e 2u>a s 



(4.20) 



Thus, the saturation scale in the tail of the hadron distribution is still in a factorized form, 
but the exponential slopes are different as compared to the corresponding form near the center, 
eq. (PD , both for the increase with r — which is now controlled by the BFKL exponent 
2lo ~ 5.55 > c — and for the decrease with b — where the "anomalous dimension" 7 of 
eq. ( 4.18 ) has been now replaced by 2. 

These changes can be easily understood by reference to eq. ( |4.5[ ): The exponent there must 
vanish when Q 2 = Q 2 (r,b). If b satisfies the condition ( 4.19|) , then the first two terms in the 
exponent, which are the large terms, almost cancel each other, so the other terms there must 
be relatively small, in the sense of eq. ( |4.2j ). Then, the term quadratic in ln(Q 2 /A 2 ) is much 
smaller than the linear term, and can be neglected. We thus end up with 



N T (Q 2 , b±) exp < —2m n b + uja s T 





u;a s T— 2m 7r ft 



(4.21) 



which, together with the saturation criterion (|2.11 ), provides indeed the expression ( 4.20 ) for 
the saturation scale. To summarize, when the "diffusion" term in the BFKL solution Q4.4j ) 
becomes negligible, then the energy dependence and the (^-dependence of the solution are fully 
controlled by the "genuine" BFKL saddle-point at Ao = 1/2. 

The change of behaviour from eq. (f4.18| ) to eq. ( |4.2C| ) is also visible on a logarithmic plot 
of the saturation scale in eqs. ( [4.13|) -( [4,14| ) as a function of b. In Fig. ^, we have displayed 
the function /9 s (r, b) of eq. ( 4.13[ ) as a function of b/Rnir), together with its small-distance and 
long-distance approximations, as given by eq. ( |4.17|) , and prior to eq. ( 4.20| ), respectively. As 
explicit on this figure, the transition between the two regimes is rather smooth, and takes place 
at intermediate values b ~ Rh(t)/2. 



4.4 Expansion of the black disk 

Let us finally consider the implications of the previous results on Q s for the expansion of the 
black disk. This is interesting since, as we shall discover in Sect. 5 below, the two domains 
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(I) and (II) are characterized by different "geometric scaling" laws for the black disk radius 
R(t, Q 2 ), and thus for the dipole cross-section. 

Equation ( 4.1 g| ) together with the definition ( p. 10 ) of R(t,Q 2 ) imply (recall that 7 = 1/A S ): 

~ A s (W-ln|J). (4.22) 



2m 7r i?(r, Q 2 



This is valid as long as R(t,Q 2 ) satisfies the condition ( 4. 16| ) , in practice, for R(t,Q 2 ) ^ 
Rh(t)/2. For a given Q 2 3> A 2 , this happens only within an intermediate range of energies, to 
be specified shortly. Of course, eq. ( 4.22] ) is just an approximate form of the general expression 
( |4.6| ) valid in this intermediate range of energies, but the approximations necessary to derive 
eq. ( [4.22| ) may not be easily recognized at the level of eq. (4.6). These approximations will be 
clarified in Sect. 5 below. 

But at sufficiently high energy, the border of the black disk lies in domain (II), that is, it is 
relatively close to the edge of the hadron, in the sense of eq. ( [4.19| ). To verify this, one can use 
eq. ( |4 .11 ) to deduce that the ratio: 



R h (t)-R(t,Q 2 
Rh{t) 



1 , Q 2 

in —pr 

2uja.T A 2 



(4.23) 



is decreasing with t, and therefore necessarily satisfies the condition ( 4.19| ) for sufficiently large 
r. Also, it can be easily checked that the expressions ( |4.2C| ) for the saturation scale in domain 
(II) and ( |4.10| ) for the black disk radius at high energy are consistent with each other, via 
eq. (|2TC| ). This is as it should since both expressions have been obtained via the same high energy 
approximation, namely, they follow from the asymptotic form ( |4.21 ) of the BFKL solution. 
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higher energy 




X"<X' 




\6(t)V / 



V < X 



Figure 5: A pictorial representation of the expansion of the black disk with increasing rapidity. 
The dotted line circle of radius b(r) = Rh{t)/2 separates between domains (I) and (II). 
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The evolution of the black disk with increasing r at fixed Q 2 is pictorially illustrated in Fig. 
The black disk appears at the center of the hadron at a critical rapidity fi(Q 2 ) such that 



Qj(n,b = 0) = Q . This condition together with eq. (4.8) implies: 



HQ 2 ) = T^ln^. (4.24) 



ca s A 2 



For t > fi(Q 2 ), but not much larger, the black disk remains confined to domain (I), as illustrated 
by the smallest disk on the left of Fig. ||. However, the expansion rate of the black disk, which 
is equal to c\ s m 3.12 in appropiate units (cf. eq. fl4.22j )), is faster than the corresponding rate 
co/2 ~ 1.39 for the borderline 5(r) = Rh(t)/2 between the two domains. Thus, at some new, 
larger, critical value, which can be easily estimated from eq. ( [4,22| ) as 

r 2 (Q 2 ) = )— — In % , (4.25) 

the black disk reaches domain (II), and than extends further within this domain (as illustrated 
by the two larger disks in Fig. ||). But within domain (II), the expansion rate of the black 
disk slows down to uj ~ 2.77, cf. eq. ( 4.1CQ , which is the same rate as for the hadron outer 



border Rh{t). Thus, with increasing r, the radial distance between the border of the black disk 
and the edge of the hadron (i.e., the width of the grey area) remains constant (cf. eq. ( [4.11[) ), 
so the relative size of the grey area with respect to the hadron size is smaller and smaller (cf. 
eq. ( |4.23| )). At large r, this grey area represents the tail of the hadron wavefunction, in which 
the BFKL equation remains valid even at arbitrarily large energy. 

Given the (quasi)exponential decrease of the saturation scale Q 2 (r,b) with b, it is straight- 
forward to verify that the conditions A 2 <C Q 2 (t, b) <C Q 2 hold at any point b± in the corona 
R(t, Q 2 ) <C b± <C Rh(t), as necessary for the consistency of our approximations. Indeed, these 
conditions are satisfied as soon as the separation between b± and the (inner or outer) edges of 



this corona is of order 1/ra, or larger. For instance, eq. (4.20) can be rewritten as: 

Q 2 (r,b) ~ A 2 e WRH(r)-&) _ ( 42g ) 

Similarly, by using eq. ( p.10 ), the expressions ( |4.18| ) and ( [4.20 ) can be recast into the form: 



Q s (r,b) ^ Q e~ ^ mn( , (4.27) 

where 7 ~ 1.55 if b± is in domain (I), and 7 = 2 for b± in domain (II). 



5 Geometric scaling at high energy 

"Geometric scaling" refers to the property of the dipole-hadron total cross section a(r,Q 2 ), 
eq. (p.5|), to depend upon the two kinematical variables r and Q 2 only via the combination 
Qq(t)/Q 2 (the "scaling variable") where Qq{t) is some suitable momentum scale which increases 
as a power of the energy: Qq(t) tx e Ar . This property is interesting since it can be related to a 
similar property of the virtual photon total cross section cr 7 * p (x, Q 2 ) which is actually seen in 
the HERA data on deep inelastic scattering (for x < 0.01 and Q 2 < 400 GeV 2 ) @. 

Clearly, this property cannot hold for arbitrary r and Q 2 , since it is known to be violated 
by the solutions to the linear evolution equations (DGLAP or BFKL) at very high Q 2 . On 
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the other hand, in the saturation regime at Q 2 < Q 2 (t), this property is physically motivated, 
since the saturation scale is then the only scale in the problem |2(J (at least for a homogeneous 
hadron), so it is tempting to identify it with the momentum scale Qq(t) introduced above. 

So far, studies of geometric scaling have been performed only for a homogeneous hadron, so 
they have naturally focused on the corresponding property of the scattering amplitude M T {Q 2 ) 
(because a(r,Q 2 ) = 2irR 2 N T (Q 2 ) in this case). Such previous studies — which relied on either 
numerical [39|, or (approximate) analytic p9| , |30| solutions to the BK equation — not 

only confirmed the existence of geometric scaling in the saturation regime, but also showed that 
this property extends up to momenta Q 2 considerably larger than Q 2 (t). In particular, in Ref. 
p9| , the upper limit for "extended scaling" has been estimated as Q 2 max ~ Qs(t)/A 2 , which is 
roughly consistent with the phenomenology |3^] (see also below). 

Our purpose in this section is to extend the analysis of Ref. [^9| by taking into account the 
transverse inhomogeneity in the hadron, within the formalism developed in the previous sections 
of this paper. Note that, since the impact parameter is integrated over in the formula (^5|) for 
<t(t, Q 2 ), the connection between the scaling properties of the cross-section and those of the 
scattering amplitude M T (r±, b±) becomes more subtle now. In particular, it is not a priori clear 
what should be the scale which plays the role of Qq(t) in the scaling variable, or even whether 
geometric scaling exists at all (since the inhomogeneous problem is not a single-scale problem 
any longer). 

To progressively introduce the effects of the inhomogeneity, let us start with a situation where 
the rapidity r is small enough for the hadron to look "grey" everywhere: Q 2 > Q 2 (t, 6 = 0), or 
r < fi(Q 2 ), cf. eq. (|4.24| ). In that case, and within the present approximations, the scattering 
amplitude M T (Q 2 ,b±) is factorized at any b± : M T {Q 2 , b± ) ~ ftf T (Q 2 )S(b±). Thus, the study of 
the scaling properties of the total cross-section: 

a(r,Q 2 ) ~ 2Af T (Q 2 ) J d 2 b ± S(b ± ) (all grey), (5.1) 



reduces to the corresponding study of M T {Q 2 ) in the homogeneous case p9| . It is instructive to 
briefly rederive here the main result in Ref. [^9| (in a very schematic way): 

The homogeneous scattering amplitude Af T (Q 2 ) is given by (|4.5|) with 6 = 0. As such, this 
formula shows no scaling. To see the scaling emerging, we shall replace the arbitrary reference 
scale A 2 in this equation with the saturation scale Q 2 (t) = A 2 e CQsT (which is the same as the 
central saturation scale Q 2 (r, 6 = 0) in the inhomogeneous case; cf. eq. (f4.8[)). We write: 



, Q 2 , Q 2 , Q 2 s (r) . Q 2 , _ 

A2 = Q|(r) ~h?~ = ln QW) + CasT - ( } 



Since, by construction, the saturation scale is such that the exponent in ( [4.5D (with 6 = 0) 
vanishes for Q 2 = Q 2 (t), it is clear that, after the replacement (5^), we are left only with terms 
involving, at least, one power of ln(Q 2 /Q 2 (t)) : 

K(C?) = «P< -A, b - ^ (in^) ) ■ (5.3) 



Here, A s rj 0.64 has been generated as X s = 1/2 + c/f3 (cf. eqs. ((Oj) and ( 4.17| )). So far, (|5 



is just a rewriting of eq. ([4 . 5[) . But this is suggestive of the conditions under which geometric 
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scaling should be expected: This emerges when the (central) saturation scale Q 2 (t) is sufficiently 
close to Q 2 (although still smaller than it) for the quadratic term in eq. (|5,3|) to be negligible 
compared to the linear term. When this happens, eq. (|5.3|) can be approximated as: 



K(Q 2 ) - i^) S (5-4) 

which shows geometric scaling indeed, with Qq(t) = Q 2 (t). This is valid as long as: 
O 2 

1 < ln 775~TT « 2A,/?a s r, (5.5) 

where the lower limit is simply the condition that we are in a regime where eq. (|Tl|) applies: the 
hadron looks "grey" everywhere. Since 2\ s f3 ~ 43.37 is a large number, eq. ( |5.5| ) gives a rather 
large window, which however extends beyond the validity range of the BFKL saddle point 
A more complete analysis shows that eq. (5J5) should be replaced by^: 



i^l Q2 ^ l Q * {t) ~ «\ 

1 < In Q 2 ( T ) < m A 2 = ca « T ' ( 5 - 6 ) 

which for the present purposes is rewritten as a range for r: 

ln TT < T < -r ln T?- ( 5 - 7 ) 



2ca s A 2 ca s A 2 ' 

The lower limit in the equation above, which arises from the upper limit in eq. ( |5.6| ), is the 
smallest value of r at which the BFKL solution starts to behave like a scaling function, for 
a given Q 2 . As for the upper limit — in which we recognize the critical value fi(Q 2 ) for the 
emergence of the black disk, eq. ( [4.24 ) — , this is necessary only for the validity of eq. ([O]) . For 



higher values of r, geometric scaling may still hold, but in order to see it, the calculation should 
be modified to account for the formation of the black disk. 

Specifically, for r > f\(Q 2 ), the cross-section (|2.5|) can be then decomposed into a "black" 
contribution plus a "grey" one, which are evaluated as (for R(t,Q 2 ) > l/m n ): 

cr(r, Q 2 ) ~ 2ttR 2 (t,Q 2 ) + 2M T (Q 2 ) J d 2 b ± S(b ± ) Q(b - R(t,Q 2 )) 

. 2,R 2 (r,Q 2 ) + 2 -^^l. (5.8) 



The second line is obtained from the first one after replacing S(b) ~ e m , and noticing 
that Af T (Q 2 )e~ 2rri7TR ( T 'Q ) = 1. Eq. ( |5.8| ) confirms that, when the energy is large enough for 
R(t,Q 2 ) > I/to^, the total cross-section is dominated by the black disk, as anticipated in 
Sect. 2. It further shows that the scaling properties of the cross section are determined by 
the corresponding properties of the radius of the black disk, which can be inferred from the 
discussion in Sect. 4. As in Sects. 4.3-4.4, we are led to distinguish between two regimes: 



12 Note that the strong inequality on the logarithm in eq. J5.5| ) has been replaced in eq. (5.6) by a normal 
inequality, which is equivalent to a strong inequality on the argument of the log. This is the condition Q 2 <C 
Ql(r)/A 2 alluded to before H. 
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I) If ti(Q 2 ) < r < t 2 (Q 2 ), cf. eqs. ( |4.24| )-( p~25|), t he black disk lies entirely inside domain 
(I), and the corresponding radius is given by eq, ( f4,22j ), which is now rewritten asP| : 



R(t, Q 2 



X s , Q 2 s (r,b = 0) 



{I) 



2771^ 



In- 



Q 2 



(5.9) 



This shows scaling, with the same scaling variable Q 2 (t)/Q 2 as in the "all grey" regime (J5T 

II) In the high-energy regime r > f2(Q 2 ), the edge of the black disk is in domain (II), so its 
radius is given by eq. ( 4.1 0| ) . This shows scaling too, but with a different scaling variable: 



R(t, Q 2 



1 ^QM 



2 n 2wa a T 



(II) 4m n "' Q 2 ' 
This gives the scaling law for the total cross section at very high energy: 



<t(t, Q 2 ) ~ 2ttR 2 (t, Q 2 ) 



7T 



2mi 



In 



Ql(r) 

Q 2 



(5.10) 



(5.11) 



To summarize, for fixed Q 2 and intermediate energies corresponding to the following range 

Q 2 



of rapidities (cf. eqs. (5/7) and ( 4.25 )): 



i . Q 2 

2ca a A 2 



< r < 



In 



(c - uj/2X s )a s A 2 



= HQ 1 



(5.12) 



the dipole-hadron cross-section exhibits geometric scaling, with the scale set by the central satu- 
ration scale Q 2 (t) = Q 2 (r,b = 0). Eq. (5.12) allows for a significant window since, numerically, 
2c ~ 9.68 while c — u/2X s ~ 2.69. On the other hand, at higher energies, corresponding to 
r > T2(Q 2 ), there is scaling again, but the relevant scale is rather the asymptotic scale Q^o(t). 

Of course, the scaling breaks down, strictly speaking, for r ~ f2(Q 2 ), i-e., in the transition 
regime where the black disk crosses from domain (I) to domain (II). But it so happens that, 
numerically, the difference between the exponents c ~ 4.84 and 2uj ~ 5.55 in the corresponding 
scaling variables (Q 2 (t) and, respectively, Q^ir)) is quite small, so the transition from one 
scaling law to the other takes places rather fast, as also manifest in the plot in Fig. |||. 

It should be also stressed that, with incresing r, the geometric scaling becomes less and 



less relevant, since the (5 2 -dependence of the cross-section ( 5.11 ) eventually becomes a sublead- 
ing effect: When the condition ( f4.2j ) is fulfilled, the leading-order term in eq. ( |5.11| ) is simply 
proportional to r 2 , with the scale set by the pion mass (cf. eq. ( 2.1 7| ) ) . 

For completness, let us conclude with a discussion of the scaling properties of the local 
scattering amplitude N T (Q 2 , b±). Within the grey area, the factorization property J\f T (Q 2 , b±) = 
S(b±)M T (Q 2 ) implies that any scaling property of the homogeneous solution M T (Q 2 ) transmits 
automatically to M T {Q 2 , b±), but with a scaling variable involving the local saturation scale 
Q 2 (r,b) : 



M T {Q 2 



QM 
Q 2 



K(Q 2 ,b±) 



Q 2 s (r,b) 

Q 2 



(5.13) 



Incidentally, this also shows that the way to derive ec 



ulations similar to those leading from eq. (4.5) to eq. (5.4), cf. eqs. ( 5.2 )— ( 5.3 ) 



(|4,22|) fr om t h e gen eral expression (4.6) is via manip- 
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with Qs(t, b) = (5o('r)e _2m7rb/ ' A . This is precisely the scaling solution (3.11) that we used in our 
arguments in Sect. |3.1| . In eq. ( 5.13j ), the momentum scale Qq(t) and the power A depend upon 
the rapidity r, in the expected way: When r is in the intermediate range ( |5.12| ), Qq(t) is the 
saturation scale at the center (J 2 (t, b = 0), and A = A s 0.64. At higher rapidities, r > f2(Q 2 ), 
one rather has Qq(t) = Q^r) and A = 1/2. 

In both cases, eq. ( |5.13| ) holds only in the grey area at b > R(t,Q 2 ). But inside the black 
disk, geometric scaling holds as well, and almost trivially, since there A/" T (Q 2 , ~ 1 and the 



deviation from one shows scaling too 



We thus conclude that, for any r above the 



lower limit in eq. (|5. 
in the hadron disk 



the scattering amplitude M T (Q 2 , b±) shows geometric scaling everywhere 



6 Conclusions and outlook 

In this paper, we have proposed a simple mechanism combining perturbative gluon saturation and 
non-perturbative boundary conditions which ensures the saturation of the Froissart bound in dipole- 
hadron scattering at high energy. Gluon saturation has been implemented via the non-linear evolution 
equations derived in perturbation theory in Refs. |7|, & ||, [Top , that we have used only for impact 
parameters in the central region of the hadron, where the gluon density is high and perturbation theory 
is applicable. On the basis of these equations, we have shown that, with increasing energy, the non- 
linear effects ensure not only the unitarization of the scattering amplitude at fixed impact parameter, 
but also the factorization of the impact parameter dependence of the scattering amplitude in the 
outer "grey area", where the unitarity limit has not yet been reached. This factorization, together 
with the exponential fall-off of the non-perturbative initial condition at large distances, leads to a 
total cross-section which grows like In 2 s. The coefficient of this growth is universal, i.e., independent 
of the hadronic target, and has been computed here. 

Our analysis makes explicit the deep connection between unitarization effects and the forma- 
tion of the Colour Glass Condensate. The "black disk" within which the unitarity limit has been 
reached is precisely the region of the hadron where the gluons "seen" by the incoming dipole are 
saturated, whereas the "grey area" corresponds to a region of lower density, in which the BFKL 
evolution still applies. The transition from "grey" to "black", i.e., the expansion of the black 
disk, is described by the non-linear evolution equations for the Colour Glass Condensate. The 
colour correlations at saturation are essential for both factorization and unitarization: They en- 
sure colour neutrality for the saturated gluons, and thus suppress the non-unitary contributions 
due to the long-range Coulomb tails. 

We have given the first computation of the saturation scale for an inhomogeneous hadron 
and identified two factorization regimes at different impact parameters. This has interesting 
consequences for the geometric scaling properties of the total cross-section: Our analysis predicts 
two different scaling laws in different ranges of energy. 

For simplicity, our analysis has been carried out for an external probe which is a qq dipole 
of fixed transverse size. But its extension to deep inelastic scattering should be straightforward. 
To this aim, the dipole-hadron cross-section must be averaged over all transverse sizes of the qq 



pair, with a weight given by the light-cone wavefunction of the virtual photon [24, 25, 26]. 

It is important to point out two limitations of our approach, which call for further studies. 
First, it was essential for our approximation scheme that the dipole is small, r± -C 1/Aqcd- In 
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particular, this restricts the applicability of our results to deep inelastic scattering at relatively 
large photon virtuality, Q 2 3> Aq CD . It would be most interesting to see if one can phenomeno- 
logically extend this analysis towards the strong coupling regime at Q 2 ~ Aqcd using the " soft 
pomeron" , and thus get at least an estimate for the coefficient of the In 2 s growth in that regime. 

Second, we have performed our analysis at fixed coupling. It is an interesting open question 
how our results would be modified by the running of the coupling. Strictly speaking, this running 
is one of the next-to-leading order effects which so far have not been systematically included 
in the non-linear evolution equations (see, however, ^0|)- But the experience with the BFKL 
equation, for which the NLO corrections have been recently computed El]], suggests that a 
significant part of these corrections could be indeed taken into account by including the running 
of the coupling. Then, the natural question is, what should be the scale at which the coupling 
must be evaluated. 

In a previous paper f29fl , where the saturation scale Q s (t) for a homogeneous hadron has been 
computed, we found it natural to evaluate the coupling at the saturation momentum. Indeed, in 
that problem, the rapidity t and the external momentum Q 2 (the resolution of the dipole) were 
increased simultaneously, in such a way to preserve the saturation condition Q 2 (t) = Q 2 . For 
that particular running, we have shown that the saturation momentum changes its parametric 
dependence upon rapidity with respect to the fixed coupling case, from e casT to e K ^+^ (with 
constant k and to). This result has been subsequently confirmed in Ref. ppf . 

On the other hand, in the present paper the physical situation is quite different: The external 
momentum Q 2 is fixed, and with increasig r, we simultaneously increase the impact parameter 
b, in such a way that the condition Q 2 (r,b) = Q 2 (which defines the edge of the black disk) 
remains satisfied. It seems therefore natural to evaluate the coupling at the fixed external scale 
Q 2 , in which case the fixed coupling expressions obtained in this paper would remain valid 
(after the trivial replacement a s — > a s (Q 2 )). A further argument in this sense is provided by the 
discussion in Sect. 3.1, which shows that the dominant scattering involves only nearby colour 
sources, within the area covered by the incoming dipole or slightly further away, so that the 
typical transferred momenta are of order Q 2 . 

A still different possibility for the running, which would be closer in spirit to Ref. [29], would 
be to evaluate the coupling at the local saturation scale Q 2 (r,b). Although this scale and the 
external scale Q 2 are identified at the edge of the black disk, they are nevertheless different 
at the points b± which lie further away in the grey area (where Q 2 > Q 2 s (r,b)). Thus, this 
choice for the running would probably modify the current formulae, in such a way to provide a 
generalization of the results in Ref. (2^] to the case of an inhomogeneous hadron. 

Of course, the description above is extremely crude and exploratory, and a full analysis of 
the running coupling case is required before any strong conclusions can be drawn. 
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Note added 

After the submission of the original version of this paper, there appeared a preprint by Kovner 
and Wiedemann which disputes the claims of this paper |^ |. In Sec. 3.2, we have added some 
clarifying comments which we hope will ameliorate this dispute. In this note, we would like to 
discuss this in more detail. 



We are first of all pleased to notice that the authors of Ref. [42| have to some extent agreed 
with our criticism to their papers [jll], and modified their arguments accordingly. For instance, 
they now agree that the dominant contribution to the scattering within the "grey area" comes 
from local sources (although their exact definition of the "grey area" is somewhat different 
from ours). Also, they admit that the long-range, non-unitary, contribution to the scattering 
amplitude that they previously computed in Ref. []l4[ applies only for a coloured probe, and not 
for a dipole. 

But in spite of their agreement on these points, the authors of Ref. |4^] dispute our conclusion 
on "Froissart bound from gluon saturation", and instead conclude that the Froissart bound is 
violated by the perturbative evolution equations in Refs. ||], || pLQ] . In its essence (but using 
the language of our Sect. 3.2), the argument of Ref. [ |42| can be formulated as follows: Assume 
that one starts the perturbative evolution at some early rapidity to, at which the point b± of 
interest lies far outside the initial "grey area" (b 3> Rh(to)), i.e., in the "white area", according 
to the terminology of Ref. In this perturbative setting, the dominant contribution to 

the scattering amplitude at b± is the long-range contribution, which decreases only as 



cf. eq. (p4p . With increasing r, the central region expands, and the point b± is eventually 



incorporated in the "grey area" at some "time" t%. For r > t%, the evolution proceeds locally 
at b±, but with an initial condition at r = t\ which decreases as a power of b±, rather than as 
an exponential. Thus, the first term in the r.h.s. of eq. ( 3.34| ) is replaced by e UJOlsT /b 4 L , and the 



"blackness" condition M T ( r ±,b±) ~ 1 yields a black disk increasing exponentially with r. 

Although mathematically correct (at least, in its above formulation; we disagree with some 



of the intermediate steps and detail points in Ref. [42]), this argument is physically irrelevant: 
It simply signals a pathology of perturbation theory when this is used in a region where it is 
not supposed to work, namely, in the "white area", where physics is truly controlled by the 
confinement. In the discussion in Sect. 3.2, we pointed out that, under plausible assumptions 
about the non-perturbative physics in this region, the true evolution respects the exponential 
fall-off of the initial condition. 

In any case, our point in this paper was not to argue that, by integrating the perturbative 
evolution up to arbitrarily large r, one would generate cross-sections which saturate the Froissart 
bound indefinitely. Rather, we have shown that by applying perturbation theory where it is 
expected to work — namely, in the central region at b < Rh(t), where the condition Q s {t, b) 3> 
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Aqcd is satisfied by definition (cf. eq. ( |2,12| )) — , one finds that (i) it is possible to compute the 
expansion rate of the black disk, and (ii) this rate is such that the Proissart bound is saturated. 
None of these points is trivial: 

(z) It was not a priori obvious that the expansion of the black disk is controlled by pertur- 
bation theory. We have demonstrated this for the case of a small dipole (Q 2 3> Aq CD ), for 
which we have shown that the "grey area" (= the region in which Q 2 (r, b) decreases from Q 2 
down to Aq CD ) has a considerable extent (cf. eq. ( |4.11[ )). This is the region which controls the 
expansion, and perturbation theory is applicable here. 

(ii) It was not a priori clear that the perturbative evolution of the black disk will respect the 
Froissart bound. Indeed, the saturated gluons from the black disk generate long-range forces 
which remain perturbative at impact parameters within the "grey area", and thus cannot be 
discarded by invoking confinement. (This was precisely the objection to Froissart bound raised 
in Ref. [14].) We have shown that these long-range forces, although present, are nevertheless 
harmless: Because the saturated gluons are globally colour-neutral (an aspect which has been 
overlooked in Ref. [jl4|), they generate only dipolar forces, whose contribution to the scattering 
amplitude in the grey area remains smaller than the corresponding contribution of the local 
sources (which saturates the Froissart bound). 

The only place where non-perturbative physics has entered our argument was in providing 
the b± -dependence of the initial condition for the short-range contribution. On physical grounds, 
this is fixed by confinement, and thus is certainly of the exponential type (and not of the power- 
law type, as a naive perturbative evolution starting in the "white area" would predict). 

To conclude, Kovner and Wiedemann seem to be addressing a different question: Whether or 
not one can naively apply the perturbative evolution equation in all regions of impact parameter 
space — including over distances many times the pion Compton wavelength — , and still get 
sensible results. They negatively answer this question, and we definitely agree with them on 
this point. On the other hand, they further make the statement that one cannot compute the 
high-energy cross section using the ideas associated with (perturbative) saturation. The body 
of this paper clearly contradicts this claim. 
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